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1 Abstract 

> 

Measurements of various features of the fluxes of atmospheric and solar neutrinos 
have provided evidence for neutrino oscillations and therefore for neutrino masses 
and mixing. We review the phenomenology of neutrino oscillations in vacuum and 
in matter. We present the existing evidence from solar and atmospheric neutrinos 
as well as the results from laboratory searches, including the final status of the 
LSND experiment. We describe the theoretical inputs that are used to interpret 
the experimental results in terms of neutrino oscillations. We derive the allowed 
ranges for the mass and mixing parameters in three frameworks: First, each set 
of observations is analyzed separately in a two-neutrino framework; Second, the 
data from solar and atmospheric neutrinos are analyzed in a three active neutrino 
framework; Third, the LSND results are added, and the status of accommodating 
all three signals in the framework of three active and one sterile light neutrinos is 
presented. We review the theoretical implications of these results: the existence of 
new physics, the estimate of the scale of this new physics and the lessons for grand 
unified theories, for supersymmetric models with i?-parity violation, for models of 
extra dimensions and singlet fermions in the bulk, and for flavor models. 
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I. INTRODUCTION 

In 1930 Wolfgang Pauli postulated the existence of the neutrino in order to reconcile data 
on the radioactive decay of nuclei with energy conservation. In radioactive decays, nuclei 
of atoms mutate into different nuclei when neutrons are transformed into slightly lighter 
protons with the emission of electrons: 

neutron — > proton + electron + antineutrino . (1) 

Without the neutrino, energy conservation requires that the electron and proton share the 
neutron's energy. Each electron is therefore produced with a fixed energy while experiments 
indicated conclusively that the electrons were not mono-energetic but observed with a range 
of energies. This energy range corresponded exactly to the many ways the three particles 
in the final state of the reaction above can share energy while satisfying its conservation 
law. The postulated neutrino had no electric charge and, for all practical purposes, does not 
interact with matter; it just serves as an agent to balance energy and momentum in above 
reaction. In fact, Pauli pointed out that for the neutrino to do the job, it had to weigh less 
than one percent of the proton mass, thus establishing the first limit on neutrino mass. 

Observing neutrinos is straightforward — in principle. Pauli had to wait a quarter of a 
century before Fred Reines and Clyde Cowan Jr. observed neutrinos produced by a nuclear 
reactor. In the presence of protons, neutrinos occasionally initiate the inverse reaction of 
radioactive decay: 

v + p -> n + e + . (2) 

Experimentally, one exposes a material rich in protons to a neutrino beam and simply looks 
for the coincident appearance of an electron and a neutron. In the alternative possibility 
where the incident neutrino carries muon flavour, a muon appears in the final state instead. 

By the 1960's neutrino beams passed from a futuristic dream to one of the most important 
tools of particle physics. The technique contributed in important ways to the discovery of 
quarks, the constituents of protons and neutrons. The technique is conceptually simple, 
though technologically very challenging. A very intense beam of accelerated protons is shot 
into a beam-dump which typically consists of a kilometer-long mound of earth or a 100 meter- 
long block of stainless steel. Particle physics does the rest. Protons interact with nuclei in 
the dump and produce tens of pions in each collision. Charged pions decay into a muon and 
a neutrino. The material in the dump will eventually absorb the muons, photons and any 
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other charged particles, so that only neutrinos exit at the opposite end, forming an intense 
and controlled beam. 

Neutrinos are also produced in natural sources. Starting at the 1960's, neutrinos pro- 
duced in the sun and in the atmosphere have been observed. As we will see in this review, 
these observations play an important role in understanding the detailed features of the 
neutrinos. In 1987, neutrinos from a supernova in the Large Magellanic Cloud were also 
detected. 

The properties of the neutrino and in particular the question of its mass have intrigued 
physicists' minds ever since it was proposed (Kayser, Gibrat-Debu and Perrier, 1989; Ra- 
mond, 1999). In the laboratory neutrino masses have been searched for in two types of 
experiments (Boehm and Vogel, 1987): (i) direct kinematic searches of neutrino mass, of 
which the most sensitive is the study of tritium beta decay, and (ii) neutrinoless double-/? 
decay experiments. Experiments achieved higher and higher precision, reaching upper lim- 
its for the electron-neutrino mass of 1CT 9 the proton mass, rather than the 1CT 2 originally 
obtained by Pauli. This raised the question of whether neutrinos are truly massless like 
photons. 

Can one go further below the eV scale (that is 10~ 9 the proton mass) in the search for 
neutrino masses? This is a very difficult task in direct measurements. In 1957, however, 
Bruno Pontecorvo realized that the existence of neutrino masses implies the possibility of 
neutrino oscillations. This phenomenon is similar to what happens in the quark sector, where 
neutral kaons oscillate. Flavor oscillations of neutrinos have been searched for using either 
neutrino beams from reactors or accelerators, or natural neutrinos generated at astrophysical 
sources (the Sun giving the largest flux) or in the atmosphere (as the byproducts of cosmic 
ray collisions). The longer the distance that the neutrinos travel from their production point 
to the detector, the smaller masses that can be signaled by their oscillation. Indeed, the 
solar neutrinos allow us to search for masses that are as small as 1(T 5 eV, that is 1(T 14 of 
the proton mass! 

In fact, in recent years, experiments studying natural neutrino fluxes have provided 
us with the strongest evidence of neutrino masses and mixing. Experiments that mea- 
sure the flux of atmospheric neutrinos have found results that suggest the disappearance of 
muon-neutrinos when propagating over distances of order hundreds (or more) kilometers. 
Experiments that measure the flux of solar neutrinos (Bahcall, 1989) found results that 
suggest the disappearance of electron-neutrinos while propagating within the Sun or be- 
tween the Sun and the Earth. The disappearance of both atmospheric z/^'s and solar z/ e 's is 
most easily explained in terms of neutrino oscillations. As concerns experiments performed 
with laboratory beams, most have given no evidence of oscillations. One exception is the 
LSND experiment which has observed the appearance of electron anti-neutrinos in a muon 
anti- neutrino beam. This signal has not been confirmed so far by any other experiment. 

What can we learn from measurements of neutrino masses about our theories of particle 
physics? The Standard Model (SM) of particle physics is a mathematical description of 
the strong, weak and electromagnetic interactions. Since it was conceived in the 1960's 
by Glashow, Salam and Weinberg, it has successfully passed numerous experimental tests. 
In the absence of any direct evidence for their mass, neutrinos were introduced in the SM 
as truly massless fermions for which no gauge invariant renormalizable mass term can be 
constructed. Consequently, in the SM there is neither mixing nor CP violation in the lepton 
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sector. Therefore, experimental evidence for neutrino masses or mixing or leptonic CP 
violation provides an unambiguous signal for new physics (NP). 

The SM prediction of massless neutrinos is an accidental fact: unlike photons, no pro- 
found principle protects them from having a mass. On the contrary, modern elementary- 
particle theories anticipate ways in which they have small, but definitely non-vanishing 
masses. In fact, there are good theoretical reasons to expect that neutrinos are massive but 
much lighter than all the charged fermions of the SM. Specifically, it is very likely that neu- 
trino masses are inversely proportional to the scale of NP. Consequently, if neutrino masses 
are measured, we can estimate the relevant new scale. 

All dimensionless flavor-blind parameters of the SM, that is, the three gauge couplings 
and the quartic Higgs coupling, are of order one. In contrast, most of the flavor parameters 
— quark and charged lepton masses (except for the mass of the top quark), and the three 
CKM mixing angles — are small and hierarchical. This situation constitutes the flavor puzzle: 
Within the SM, the hierarchical structure can be accommodated but is not explained. If 
neutrinos have Majorana (Dirac) masses, then there are nine (seven) flavor parameters 
beyond the thirteen of the SM. If some of these extra parameters are measured, we will be 
able to test and refine theories that try to solve the flavor puzzle. 

In grand unified theories (GUTs), lepton masses are often related to quark masses. Mea- 
surements of neutrino parameters provide further tests of these theories. 

The values of neutrino parameters that explain the anomalies observed in atmospheric 
and solar neutrino fluxes can be used to address the many theoretical questions described 
above: they imply NP, they suggest an energy scale at which this NP takes place, they 
provide stringent tests of flavor models, and they probe GUTs. 

In this review we first present the low energy formalism for adding neutrino masses to 
the SM and the induced leptonic mixing (Sec. II) and then describe the phenomenology 
associated with neutrino oscillations in vacuum and in matter (Sec. III). In Sees. IV and V 
we discuss the evidence from, respectively, solar and atmospheric neutrinos. We review 
the theoretical modeling that is involved in interpreting the experimental results in terms 
of neutrino oscillations. We briefly describe the techniques used in the derivation of the 
allowed ranges for the neutrino flavor parameters which, in these sections, is performed in the 
framework of mixing between two neutrinos. Sec. VI is devoted to the results from searches 
at laboratory experiments including the final status of the LSND experiment. The two 
more robust pieces of evidence, from solar and atmospheric neutrinos, can be accommodated 
assuming masses and mixing of three standard neutrinos: in Sec. VILA we derive the allowed 
ranges of parameters in this case. The present phenomenological status of the possibility of 
mixing between the three standard neutrinos with a light sterile one, needed to accommodate 
also the LSND result, is discussed in Sec. VII. B. In Sees. VIII and IX we explain the various 
lessons for theory which can be drawn on the basis of these interpretations, focusing on 
models that explain the atmospheric and solar neutrino data through mixing among three 
active neutrinos. In particular, we discuss flavor models, GUTs, supersymmetric models 
with i?-parity violation, and models of extra dimensions. Our conclusions are summarized 
in Sec. X. 



5 



II. THE STANDARD MODEL AND NEUTRINO MASSES 



One of the most beautiful aspects of modern theories of particles physics is the rela- 
tion between forces mediated by spin-1 particles and local (gauge) symmetries. Within the 
Standard Model, the strong, weak and electromagnetic interactions are related to, respec- 
tively, SU(3), SU(2) and U(l) gauge groups. Many features of the various interactions are 
then explained by the symmetry to which they are related. In particular, the way that the 
various fermions are affected by the different types of interactions is determined by their 
representations (or simply their charges in the case of Abelian gauge symmetries) under the 
corresponding symmetry groups. 

Neutrinos are fermions that have neither strong nor electromagnetic interactions. In 
group theory language, they are singlets of SU(3)c x U(1)em- Active neutrinos have weak 
interactions, that is, they are not singlets of SU(2) L . Sterile neutrinos have none of the SM 
gauge interactions and they are singlets of the SM gauge group. 

The SM has three active neutrinos. They reside in lepton doublets, 

Lt=(^^, t = e,n, t. (3) 

Here e, ji and r are the charged lepton mass eigenstates. The three neutrino interaction 
eigenstates, v e , and v T , are defined as the «S£/(2)L-partners of these mass eigenstates. In 
other words, the charged current interaction terms for leptons read 

- £ cc = ^E%7%"^ + + h.c. (4) 
In addition, the SM neutrinos have neutral current (NC) interactions, 

Eqs. (4) and (5) give all the neutrino interactions within the SM. 

The measurement of the decay width of the Z° boson into neutrinos makes the existence 
of three, and only three, light (that is, m u <mz/2) active neutrinos an experimental fact. 
When expressed in units of the SM prediction for a single neutrino generation, one gets 
(Groom et al„ 2000) 

N v = 2.994 ± 0.012 (Standard Model fits to LEP data), 

N v = 3.00 ± 0.06 (Direct measurement of invisible Z width). (6) 



A. The Standard Model Implies m u = 

The SM is based on the gauge group 

G S M = SU(3) C x SU(2) L x U(1) Y , (7) 

with three fermion generations, where a single generation consists of five different represen- 
tations of the gauge group, 
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Q L (3, 2) +1/6 , U R (3, l) +2/3 , £>*(3, l)_ 1/3 , L L (1, 2)_ 1/2 , (8) 

Our notation here means that, for example, a left-handed lepton field L L is a singlet (1) of 
the S77(3)c group, a doublet (2) of the SU(2) L group, and carries hypercharge —1/2 under 
the U(1) Y group. 

The vacuum expectations value (VEV) of the single Higgs doublet 0(1, 2) + i/ 2 breaks the 
symmetry, 

(0) = ( Vj G SM -> 5C/(3) C x C/(1)em. (9) 

Since the fermions reside in chiral representations of the gauge group, there can be no 
bare mass terms. Fermions masses arise from the Yukawa interactions, 

"^Yukawa = Y^^D Rj + Y%Q^j>U Rj + Y^ i( f>E Rj + h.c, (10) 

(where <fi = iT2<f)*) after spontaneous symmetry breaking. The Yukawa interactions of 
Eq. (10) lead to charged fermion masses but leave the neutrinos massless. One could think 
that neutrino masses would arise from loop corrections if these corrections induced effective 
terms 

^4><t>L Li L Lj . (11) 

This however cannot happen, as can be easily understood by examining the accidental sym- 
metries of the Standard Model. (It often happens that, as a consequence of the symmetries 
that define a model and of its particle content, all renormalizable Lagrangian terms obey ad- 
ditional symmetries, that are not a priori imposed on the model. These are called accidental 
symmetries.) Within the SM, with the gauge symmetry of Eq. (7) and the particle content 
of Eq. (8), the following accidental global symmetry arises (at the perturbative level): 

Gf° bal = U{1) B x U(l) e x U(l), x U(1) T . (12) 

Here U(1)b is the baryon number symmetry, and U{l) e ^ T are the three lepton flavor sym- 
metries, with total lepton number given by L = L e + + L T . Terms of the form (11) violate 
Gg 1 ^ 1 and therefore cannot be induced by loop corrections. Furthermore, the U(1)b-l sub- 
group of G| 1 I ^ [ bal is non-anomalous. Terms of the form (11) have B — L = —2 and therefore 
cannot be induced even by nonperturbative corrections. 

It follows that the SM predicts that neutrinos are precisely massless. Consequently, there 
is neither mixing nor CP violation in the leptonic sector. 

B. Extensions of the Standard Model Allow m v / 

There are many good reasons to think that the SM is not a complete picture of Nature. 
For example, the fine-tuning problem of the Higgs mass can be solved by Supersymmetry; 
gauge coupling unification and the variety of gauge representations may find an explanation 
in GUTs; baryogenesis can be initiated by decays of heavy singlet fermions (leptogenesis); 
and the existence of gravity suggests that string theories are relevant to Nature. If any of 
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these (or many other proposed) extensions is indeed realized in Nature, the SM must be 
thought of as an effective low energy theory. That means that it is a valid approximation 
up to the scale A NP which characterizes the NP. 

By thinking of the SM as an effective low energy theory, we still retain the gauge group 
(7), the fermionic spectrum (8), and the pattern of spontaneous symmetry breaking (9) as 
valid ingredients to describe Nature at energies E <C Anp. The SM predictions are, however, 
modified by small effects that are proportional to powers of E/A-^p. In other words, the 
difference between the SM as a complete description of Nature and as a low energy effective 
theory is that in the latter case we must consider also nonrenormalizable terms. 

There is no reason for generic NP to respect the accidental symmetries of the SM (12). 
Indeed, there is a single set of dimension-five terms that is made of SM fields and is consistent 
with the gauge symmetry, and this set violates (12). It is given by 

Z ' J ool. (13) 



Anp 

(While these terms have the same form as (11) we are thinking now not about SM radiative 
corrections as their source but about some heavy fields, related to NP, which can induce 
such terms by tree or loop diagrams.) In particular, (13) violates L (and B — L) by two 
units and leads, upon spontaneous symmetry breaking, to neutrino masses: 

W(i = Z ~f^- (14) 
i Anp 

A few comments are in order, regarding Eq. (14): 

1) Since Eq. (14) would arise in a generic extension of the SM, we learn that neutrino 
masses are very likely to appear if there is NP. 

2) If neutrino masses arise effectively from nonrenormalizable terms, we gain an under- 
standing not only for the existence of neutrino masses but also for their smallness. The 
scale of neutrino masses is suppressed, compared to the scale of charged fermion masses, by 
v/Anp- 

3) The terms (14) break not only total lepton number but also the lepton flavor symmetry 
U(l) e x Z7(1) M x U(1) T . Therefore we should expect lepton mixing and CP violation. 

The best known scenario that leads to (13) is the see-saw mechanism (Ramond, 1979; 
Gell-Mann el al, 1979; Yanagida, 1979). Here one assumes the existence of heavy sterile 
neutrinos iVj. Such fermions have, in general, bare mass terms and Yukawa interactions: 

-C N = l -M Nij WfN 3 + Y£L^j>Nj + h.c. (15) 
The resulting mass matrix (see Sec. II. C for details) in the basis I ] has the following 



form: 



N 3 



( o Y"—\ 

If the eigenvalues of are all well above the electroweak breaking scale v, then the 
diagonalization of M v leads to three light mass eigenstates with a mass matrix of the form 
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(14). In particular, the scale Anp is identified with the mass scale of the heavy sterile 
neutrinos, that is the typical scale of the eigenvalues of M^. 

Two well-known examples of extensions of the SM that lead to a see-saw mechanism for 
neutrino masses are SO(10) GUTs (Ramond, 1979; Gell-Mann et ai, 1979; Yanagida, 1979) 
and left-right symmetry (Mohapatra and Senjanovic, 1980). 

C. Dirac and Majorana Neutrino Mass Terms 

If the only modification that we make to the SM is to assume that it is a low energy 
effective theory, that is, allowing for non-renormalizable terms that are consistent with the 
gauge symmetry and the fermionic content of the SM, then the only way that neutrinos can 
gain masses is through terms of the form (13). These are Majorana mass terms which, in 
particular, violate lepton number by two units. 

One can, however, open up other possibilities by adding new fields. The most relevant 
extension is that where an arbitrary number m of sterile neutrinos za^I, l)o is added to the 
three standard generations of Eq. (8). Now there are, in general, two types of mass terms 
that arise from renormalizable terms: 



Here v c indicates a charge conjugated field, v c = Cu T and C is the charge conjugation 
matrix. 

The first term is a Dirac mass term. It has the following properties: 

(i) Since it transforms as the doublet representation of S77(2)l, it is generated after spon- 
taneous electroweak symmetry breaking from the Yukawa interactions Y^L^v s ^ similarly 
to the charged fermion masses discussed in Sec. II. A. 

(ii) Since it has a neutrino field and an antineutrino field, it conserves total lepton number 
[though it breaks the lepton flavor number symmetries of Eq. (12)]. 

(iii) Md is a complex 3 x m matrix. 

The second term in Eq. (17) is a Majorana mass term. It is different from the Dirac 
mass terms in many important aspects: 

(i) It is a singlet of the SM gauge group. Therefore, it can appear as a bare mass term. 
[Had we written a similar term for the active neutrinos, it would transform as a triplet of 
S77(2)l. In the absence of a Higgs triplet, it cannot be generated by renormalizable Yukawa 
interactions. Such terms are generated for active neutrinos from the non-renormalizable 
Yukawa interactions of Eq. (13).] 

(ii) Since it involves two neutrino fields, it breaks lepton number conservation by two 
units. More generally, such a term is allowed only if the neutrinos carry no additive conserved 
charge. This is the reason that such terms are not allowed for any charged fermions which, 
by definition, carry C/(1)em charges. 

(iii) Mn is a symmetric matrix (as follows from simple Dirac algebra) of dimension mxm. 
It is convenient to define a (3 + m)-dimensional neutrino vector z7, 



M Dij v Li v sj + -M Nij v c si v sj + h.c. 



(17) 




(18) 
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That allows us to rewrite Eq. (17) in a unified way: 




(19) 



where 



M D 
M N 



) 



(20) 



The matrix M v is complex and symmetric. It can be diagonalized by a unitary matrix 
of dimension (3 + m). The resulting mass eigenstates, i>k, obey the Majorana condition, 



There are three interesting cases, differing in the hierarchy of scales between M N and 
M D : 

(i) One can assume that the scale of the mass eigenvalues of Mn is much higher than the 
scale of electroweak symmetry breaking ((f)). This is the natural situation in various exten- 
sions of the SM that are characterized by a high energy scale. We are back in the framework 
of the see-saw mechanism discussed in the previous section. If we simply integrate out the 
sterile neutrinos, we get a low energy effective theory with three light, active neutrinos of 
the Majorana-type. 

(ii) One can assume that the scale of some eigenvalues of M N is not higher than the 
electroweak scale. Now the SM is not even a good low energy effective theory: there are 
more than three light neutrinos, and they are mixtures of doublet and singlet fields. These 
light fields are all of the Majorana-type. 

(iii) One can assume that Mjv = 0. This is equivalent to imposing lepton number 
symmetry on this model. Again, the SM is not even a good low energy theory: both the 
fermionic content and the assumed symmetries are different. (Recall that within the SM 
lepton number is an accidental symmetry.) Now only the first term in Eq. (17) is allowed, 
which is a Dirac mass term. It is generated by the Higgs mechanism in the same way that 
charged fermions masses are generated. If indeed it is the only neutrino mass term present 
and m = 3, we can identify the three sterile neutrinos with the right handed component 
of a four-component spinor neutrino field (actually with its charge conjugate). In this way, 
the six massive Majorana neutrinos combine to form three massive neutrino Dirac states, 
equivalently to the charged fermions. In this particular case the 6x6 diagonalizing matrix 
is block diagonal and it can be written in terms of a 3 x 3 unitary matrix. 

From the phenomenological point of view, it will make little difference for our purposes 
whether the light neutrinos are of the Majorana- or Dirac-type. In particular, the analysis 
of neutrino oscillations is the same in both cases. Only in the discussion of neutrinoless 
double beta decay will the question of Majorana versus Dirac neutrinos be crucial. 

From the theoretical model building point of view, however, the two cases are very dif- 
ferent. In particular, the see-saw mechanism provides a natural explanation for the lightness 
of Majorana neutrinos, while for Dirac neutrinos there is no such generic mechanism. (The 
situation can be different in models of extra dimensions; see Section IX. B.) 
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D. Lepton Mixing 



We briefly review here our notation for lepton mixing. We denote the neutrino mass 
eigenstates by (ui, u 2 , v 3 , . . . , u n ) where n = 3 + m, and the charged lepton mass eigenstates 
by (e, 11, t). The corresponding interaction eigenstates are denoted by (e 1 ,^ 1 ^ 1 ) and V = 
iyLei ^Lfj,, vlt, Vgi, ■ ■ ■ , v sm ). In the mass basis, leptonic charged current interactions are given 
by 



'cc = ~^( e L r L )YU 



W + 



h.c. 



(21) 



Here U is a 3 x n matrix (Schechter and Valle; 1980a, 1980b). 

Given the charged lepton mass matrix and the neutrino mass matrix M v in some 
interaction basis, 



-C M = (e£ ui t[) M e 
we can find the diagonalizing matrices V e and V v : 



( e R \ l 

fJ'R + ^ c M u u + h.c. , 
\rk l 2 



(22) 



V ei M e MjV e = diag(m^, m%m 2 T ), V v] ' MlM u V v = diag(m?, m|, mj, . . .,m 2 n ). (23) 

Here V e is a unitary 3x3 matrix while V u is a unitary n x n matrix. The 3 x n mixing 
matrix U can be found from these diagonalizing matrices: 



(24) 



A few comments are in order: 

(i) Note that the indices % and k run from 1 to 3, while j runs from 1 to n. In particular, 
only the first three lines of V u play a role in Eq. (24). 

(ii) Pi is a diagonal 3x3 phase matrix, that is conventionally used to reduce by three 
the number of phases in U. 

(iii) P u is a diagonal matrix with additional arbitrary phases (chosen to reduce the 
number of phases in U) only for Dirac states. For Majorana neutrinos, this matrix is simply 
a unit matrix. The reason for that is that if one rotates a Majorana neutrino by a phase, 
this phase will appear in its mass term which will no longer be real. 

We conclude that the number of phases that can be absorbed by redefining the mass 
eigenstates depends on whether the neutrinos are Dirac or Majorana particles. In particular, 
if there are only three Majorana neutrinos, U is a 3 x 3 matrix analogous to the CKM matrix 
for the quarks (Maki, Nakagawa and Sakata, 1962; Kobayashi and Maskawa, 1973) but due 
to the Majorana nature of the neutrinos it depends on six independent parameters: three 
mixing angles and three phases. The two Majorana phases do not affect neutrino oscillations 
(Bilenky, Hosek and Petcov, 1980; Langacker, Petcov, Steigman and Toshev, 1987). [CP 
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conservation implies that the three lepton phases are either zero or tt (Schechter and Valle, 
1981; Wolfenstein 1981).] This is to be compared to the case of three Dirac neutrinos, where 
the number of physical phases is one, similarly to the CKM matrix. Note, however, that the 
two extra Majorana phases affect only lepton number violating processes. Such effects are 
suppressed by m u /E and are very hard to measure. 

If no new interactions for the charged leptons are present we can identify their interaction 
eigenstates with the corresponding mass eigenstates after phase redefinitions. In this case 
the charged current lepton mixing matrix U is simply given by a 3 x n sub-matrix of the 
unitary matrix V v . 



III. NEUTRINO OSCILLATIONS 
A. Neutrino Oscillations in Vacuum 

Neutrino oscillations in vacuum would arise if neutrinos were massive and mixed. In 
other words, the neutrino state that is produced by electroweak interactions is not a mass 
eigenstate. This phenomenon was first pointed out by Pontecorvo in 1957 (Pontecorvo, 
1957) while the possibility of arbitrary mixing between two massive neutrino states was first 
introduced in Maki, Nakagawa and Sakata (1962). 

From Eq. (21) we see that if neutrinos have masses, the weak eigenstates, v a , produced 
in a weak interaction (i.e., an inverse beta reaction or a weak decay) are, in general, linear 
combinations of the mass eigenstates z/j 

n 

kHE^h) (25) 
i=i 

where n is the number of light neutrino species and U is the the mixing matrix in Eq.(24). 
(Implicit in our definition of the state \v) is its energy-momentum and space-time depen- 
dence). After traveling a distance L (or, equivalently for relativistic neutrinos, time t), a 
neutrino originally produced with a flavor a evolves as follows: 

M0>= iXk(*)> ■ (26) 
i=i 

It can be detected in the charged-current (CC) interaction v a {t)N' — > £@N with a probability 

n n 

p aP = \(vf>\v a (t))\ 2 HEEWW )!^))! 2 > ( 27 ) 

i=l j=l 

where Ei and m 8 are, respectively, the energy and the mass of the neutrino mass eigenstate 
v%- 

We use the standard approximation that \v) is a plane wave (for a pedagogical discussion 
of the possible quantum mechanical problems in this naive description of neutrino oscillations 
we refer the reader to Lipkin, 1999, and Kim and Pevsner, 1993), \vi(t)) = e~ l Eit \i>i{Q)) . In 
all cases of interest to us, the neutrinos are relativistic: 

E t = ^ + m^p l + ^L. (28) 
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Furthermore, we can assume that pi ~ pj = p ~ E. Then we obtain the following transition 
probability (we only include here the CP conserving piece): 

n— 1 n 

P a p = <W - 4 ]T £ Re^C/^C/^.C/ft] sin 2 x y , (29) 

i=l j=i+l 



where 



Xij^AijL/2, Aij = Am|/(2£), Am?- = m 2 - m 2 , (30) 



and L = t is the distance between the source (that is, the production point of z/ a ) and the 
detector (that is, the detection point of z/3). In deriving Eq. (29) we used the orthogonality 
relation (^(0)1^(0)) = cV,-. It is convenient to use the following units: 

1 Am 2 , LIE 

xa = 1.27 ^ — / . (31) 

13 eV 2 m/MeV V ; 

The transition probability [Eq. (29)] has an oscillatory behavior, with oscillation length 



rose _ ^E ( o 2] 

\v - Am 2 w 



and amplitude that is proportional to elements in the mixing matrix. Thus, in order to 
have oscillations, neutrinos must have different masses (Am^ 7^ 0) and they must mix 

{UaPf* ± 0). 

An experiment is characterized by the typical neutrino energy E and by the source- 
detector distance L. In order to be sensitive to a given value of Am|, the experiment has 
to be set up with E/L « Am| (L ~ Lq s ^). The typical values of LjE for different types of 
neutrino sources and experiments are summarized in Table I. 

If (E/L) 3> Am 2 - (L <C Lq^j), the oscillation does not have time to give an apprecia- 
ble effect because sin 2 ^- <C 1. The case of (E/L) <C Am 2 - (L ^> L^) requires more 
careful consideration. One must take into account that, in general, neutrino beams are not 
monochromatic. Thus, rather than measuring P a p, the experiments are sensitive to the 
average probability 

dE u —a cc (E u )P a(3 (E u )e(E») 
J dE ^fcc(E v )e(E v ) (33) 

n— 1 n 
i=l j=i+l 

where <£> is the neutrino energy spectrum, o C c is the cross section for the process in which 
the neutrino is detected (in general, a CC interaction), and e(E v ) is the detection efficiency. 
For L 3> L™f j: the oscillating phase goes through many cycles before the detection and is 
averaged to (sin 2 Xjj) = 1/2. 

For a two-neutrino case, the mixing matrix depends on a single parameter, 
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TT ( cos 6 sin9\ , nA , 
U =[-sin9cos9) ' ^ 

and there is a single mass-squared difference Am 2 . Then P at3 of Eq. (29) takes the well 
known form 

P a /3 = S a/3 - (25 a(3 - 1) sin 2 29 sin 2 x . (35) 

The physical parameter space is covered with Am 2 > and < 9 < | (or, alternatively, 
< 9 < | and either sign for Am 2 ). 

Changing the sign of the mass difference, Am 2 — > —Am 2 , and changing the octant of the 
mixing angle, 9 — > | — amounts to redefining the mass eigenstates, V\ <-> z/ 2 : P Q/ 3 must be 
invariant under such transformation. Eq. (35) reveals, however, that P a p is actually invariant 
under each of these transformations separately. This situation implies that there is a two- 
fold discrete ambiguity in the interpretation of P a p in terms of two-neutrino mixing: the two 
different sets of physical parameters, (Am 2 , 9) and (Am 2 , | — 9), give the same transition 
probability in vacuum. One cannot tell from a measurement of, say, P efl in vacuum whether 
the larger component of v e resides in the heavier or in the lighter neutrino mass eigenstate. 

Neutrino oscillation experiments measure P a p. It is common practice for the experiments 
to interpret their results in the two-neutrino framework. In other words, the constraints on 
P a/ 3 are translated into allowed or excluded regions in the plane (Am 2 , sin 2 29) by using 
Eq. (35). An example is given in Fig. 1. We now explain some of the typical features of 
these constraints. 

When an experiment is taking data at fixed (L) and (E) , as is the case for most laboratory 
searches, its result can always be accounted for by Am 2 that is large enough to be in the 
region of averaged oscillations, (sin 2 ^) = 1/2. Consequently, no upper bound on Am 2 can 
be achieved by such experiment. So, for negative searches that set an upper bound on the 
oscillation probability, (P a p) < Pl, the excluded region lies always on the upper-right side 
of the [Am 2 , sin 2 26*) plane, limited by the following asymptotic lines: 

• For Am 2 » 1/(L/E), a vertical line at sin 2 29 = 2P L . 

• For Am 2 <C 1/(L/E), the oscillating phase can be expanded and the limiting curve 
takes the form Am 2 sin26 l = A^/P~l/ (L/E), which in a log-log plot gives a straight line 
of slope —1/2. 

If, instead, data are taken at several values of (L) and/or (E), the corresponding region 
may be closed as it is possible to have direct information on the characteristic oscillation 
wavelength. 



B. Neutrinos in Matter: Effective Potentials 

When neutrinos propagate in dense matter, the interactions with the medium affect their 
properties. These effects are either coherent or incoherent. For purely incoherent inelastic 
z/-p scattering, the characteristic cross section is very small: 

a ~ 3f „ 10- 43 cm 2 ( -JL-) 2 . (36) 
7T VI MeV/ v 1 
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The smallness of this cross section is demonstrated by the fact that if a beam of 10 10 
neutrinos with E ~ 1 MeV was aimed at the Earth, only one would be deflected by the 
Earth matter. It may seem then that for neutrinos matter is irrelevant. However, one 
must take into account that Eq. (36) does not contain the contribution from forward elastic 
coherent interactions. In coherent interactions, the medium remains unchanged and it is 
possible to have interference of scattered and unscattered neutrino waves which enhances 
the effect. Coherence further allows one to decouple the evolution equation of the neutrinos 
from the equations of the medium. In this approximation, the effect of the medium is 
described by an effective potential which depends on the density and composition of the 
matter (Wolfenstein, 1978). 

As an example we derive the effective potential for the evolution of v e in a medium 
with electrons, protons and neutrons. The effective low-energy Hamiltonian describing the 
relevant neutrino interactions is given by 

Hw = ^[J^ a (x)4-\x) + |^ a WfW] , (37) 

where the J a 's are the standard fermionic currents, 

J^\x)=V- e (x) la (l- lb )e(x) , (38) 

Ji~\x) = e(x) la (l- l5 )u e (x) , (39) 
J i N \ x ) = ^(x)7 Q (l - 7 5 )z/ e (x) -e(x)[y a (l -75) - 4 sin 2 9 w ^ a \e{x) 

+p(x)[ la (l - #75) - 4 sin 2 9 wla \p{x) - n{x) la {\ - g^lsXx) . (40) 

g^'^ are the axial couplings for neutrons and protons, respectively. For the sake of sim- 
plicity we concentrate on the effect of the charged current interactions. The effective CC 
Hamiltonian due to electrons in the medium is 

H { c e) = ^Jd 5 Pe f(E e ,T) (41) 
x ((e(s,p e )|e(a;)7 a (l - 7 5 )z/ e (:r)z7e(:r)7 Q (l - 7s)e(a;)|e(s,p e ))) 

= % i^(z)7a(l - ls)ve{x) J d 3 p e f(E e ,T)((e(s,p e )\e(x)-f a (l - 7 5 )e(x)|e(s,p e ))) , 

where s is the electron spin and p e its momentum. The energy distribution function of 
the electrons in the medium, f(E e ,T), is assumed to be homogeneous and isotropic and is 
normalized as 

Jd 3 p e f(E e ,T) = l. (42) 

By ^•••^ we denote the averaging over electron spinors and summing over all electrons in the 
medium. Notice that coherence implies that s,p e are the same for initial and final electrons. 
Expanding the electron fields e(x) in plane waves we find 

(e(s,p e )|e(x)7 Q (l - 7 5 )e(x)|e(s,p e )) 
= — (e(s,p e )\u;(p e )al(p e )^ a (l - -f 5 )a s (p e )u s (p e )\e(s,p e )} , (43) 
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where V is a volume normalization factor. The averaging gives 

^-{{e(s,p e )\al(p e )a s (p e )\e{s,p e ))) = N e {p e )^Yl > ( 44 ) 

v " s 

where N e (p e ) is the number density of electrons with momentum p e . We assumed here that 
the medium has equal numbers of spin +1/2 and spin —1/2 electrons, and we used the fact 
that a\(p e )a s (p e ) = M^ s \p e ) is the number operator. We thus obtain: 

((e(s,p e )|e(x)7 Q (l - 7 5 )e(x)|e(s,p e ))) = N e {p e )]^y^){p e )^ a {l - 7 5 )u(«)(p e ) 



2 



s 



2 ^l^-^H^i (45) 

Isotropy implies that / d 3 p e p" e f(E e , T) = 0. Thus, only the p° term contributes upon integra- 
tion, with / d 3 p e f(E e ,T)N e (p e ) = N e (the electron number density). Substituting Eq. (45) 
in Eq. (41) we obtain: 

4 6) = ^^(x)Tb(l - K>e(x) . (46) 

The effective potential for v e induced by its charged current interactions with electrons in 
matter is then given by 



V c = (u e \ j d 3 x Hg\v e ) = ~^fy J d 3 x u\u v = V2G F N e . (47) 



For v e the sign of V is reversed. This potential can also be expressed in terms of the matter 
density p: 

V c = V2G F N e ~ 7.6 Y e 9 3 eV , (48) 

10 i4 g/cm 

where Y e = N ^ N is the relative number density. Three examples that are relevant to 
observations are the following: 

• At the Earth core p ~ 10 g/cm 3 and V c ~ 10 -13 eV; 

• At the solar core p ~ 100 g/cm 3 and V c ~ 10 -12 eV; 

• At a supernova p ~ 10 14 g/cm 3 and Vc ~ eV. 

Following the same procedure we can obtain the effective potentials for v e due to inter- 
actions with different particles in the medium. The results are listed in Table II (Kim and 
Pevsner, 1993). For and u T , Vc = for any of these media while Vn is the same as for v e . 
One can further generalize this analysis to other types of interactions (Bergmann, Grossman 
and Nardi, 1999). 
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C. Evolution Equation in Matter: Effective Mass and Mixing 



There are several derivations in the literature of the evolution equation of a neutrino 
system in matter (see, for instance, Halprin 1986, Mannhein 1988). We follow here the 
discussion in Baltz and Weneser, 1988. Consider a state which is an admixture of two 
neutrino species \u e ) and \vx) or, equivalently, of \ui) and \v2): 

$(x) = $ e (z)k) + *x(x)\u x ) = $i(x)h> + $ 2 (z)h> (49) 

The evolution of $ in a medium is described by a system of coupled Dirac equations: 

rh d i 
E®! = -a x — + $m x + V n $i + V l2 § 2 , 
1 % ox J 

rh d i 

E<$> 2 = ~a x —+ pm 2 + V 22 $2 + V 12 $i , (50) 



where (3 = 70 and a x = 7071 ■ The Vij terms give the effective potential for neutrino 
mass eigenstates. They can be simply derived from the effective potential for interaction 
eigenstates [such as V ee of Eq. (47)]: 



Vi 



(TxH. 



medium 
int 



TJ. V U* 



(51) 



We decompose the neutrino state: $i(x) = Ci(x)fa(x). Here <f>i(x) is the Dirac spinor part 
satisfying: 



(a x {[E - V ti (x)} 2 - ml} 112 + fim, + V n ) fa(x) = £&(a;) . (52) 
So 0j(x) has the form of free particle solutions with local energy £i(x) = E — Vu(x): 



fa(x) 



iS. + m^ 2 



2£ 



X 



v 7 ^ 



L Si+rm 1 . 

where x is the Pauli spinor. We make the following approximations: 



(53) 



(i) The scale over which V changes is much larger than the microscopic wavelength of the 
neutrino: %/V < Tim/E 2 . 

(ii) Expanding to first order in V implies that V\ 2 a x fa ~ fa, V\ 2 a x fa ~ fa and {[E — 
V u (x)] 2 -m^ 2 ^E-Vdx)-§. 

From (i) we find that the Dirac equations take the form: 



EC x fa = -a x —fa + (/3mi + Ki)Ci0i + V 12 C 2 (j) 2 , 



EC 2 fa =- -<\ 



% dx 

h dc 2 



t OX 



fa + (Pm 2 + V 22 )C 2 fa + V X2 C x fa. 



(54) 



Then multiplying by a x and using the equation of motion of fa and (ii), we can drop the 
dependence on the spinor and obtain 
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~ = lE-V^)-§]C 2 -V 12Cl . (55) 
Changing notations Ci jCe (x) — > Vi tU (x), we rewrite Eq. (55) in a matrix form : 



hd_(vi\ E-V u -^ -V 12 
idx\v2) \ -V 12 E-V 22 - 



2E 




(56) 



After removing the diagonal piece that is proportional to E, we can rotate Eq. (56) to the 
flavor basis (H — 1) (Wolfenstein, 1978): 

where we have defined an effective mass matrix in matter: 



2 , ^ + 2£V e - ^cos 2£ ^sin 2tf . 

" I A^ sin2 ^ + 2 £^ + ^cos2# 1 ' 

Here Am 2 = m 2 , — m\. 

We define the instantaneous mass eigenstates in matter, v™, as the eigenstates of M w for 
a fixed value of x (or i). They are related to the interaction eigenstates through a unitary 
rotation, 

U)-w(?)-(!Su^fc)(?)- w 

The eigenvalues of M m , that is, the effective masses in matter are given by (Wolfenstein, 
1978; Mikheyev and Smirnov, 1985): 

2 2 

^(x) = m i+ m 2 + E ( Ve + Vx ) T L^A m 2 cos2 e- A) 2 + (Am 2 sm26) 2 , (60) 

while the mixing angle in matter is given by 

Am 2 sin 26 . . 

Un2e ™= A,ni C os2e-A ' < 61) 

The quantity A is defined by 

A = 2E(V e -V x ). (62) 

In Figs. 2 and 3 we plot, respectively, the effective masses and the mixing angle in matter 
as functions of the potential A, for A > and Am 2 cos 26 > 0. Notice that even massless 
neutrinos acquire non-vanishing effective masses in matter. 

The resonant density (or potential) Ar is defined as the value of A for which the difference 
between the effective masses is minimal: 
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A R = Am 2 cos 29 . (63) 

Notice that once the sign of V e — V x (which depends on the composition of the medium and 
on the state X) is known, this resonance condition can only be achieved for a given sign of 
Am 2 cos 29, i.e. for mixing angles in only one of the two possible octants. We learn that 
the symmetry present in vacuum oscillations is broken by matter potentials. Also if the 
resonant condition is achieved for two neutrinos it cannot be achieved for antineutrinos of 
the same flavor and viceversa. The mixing angle tan# m changes sign at A R . As can be seen 
in Fig. 3, for A > A R we have 9 m > 9. 

We define an oscillation length in matter: 

L osc = L ° Am (64) 



^(Am 2 cos 29 - A) 2 + (Am 2 sin 29) 2 ' 



where the oscillation length in vacuum, Lq sc , was defined in Eq. (32). The oscillation length 
in matter presents a resonant behaviour. At the resonance point the oscillation length is 

T OSC 

L° sc = . (65) 

R sm29 K ' 

The width (in distance) of the resonance, 5r R , corresponding to 5A R = 2 Am 2 sin 2 29, is 
given by 

5A R 2 tan 29 
Sr R = TdAT = I > ( 66 ) 

where we have defined the resonance height: 

1 dA 



h R 



\ dr 



(67) 

R 



For constant A, i.e., for constant matter density, the evolution of the neutrino system is 
described just in terms of the masses and mixing in matter. But for varying A, this is in 
general not the case. 



D. Adiabatic versus Non-adiabatic Transitions 

Taking time derivative of Eq. (59), we find: 

l(r;)=^)(?) + ^™)(?) • (68) 

Using the evolution equation in the flavor basis, Eq. (57), we get 

1 ( % ) = ^E U "^ M " U ^ ( % ) - * Uttf ((>rn) (J?}. (69) 

For constant matter density, 9 m is constant and the second term vanishes. In general, using 
the definition of the effective masses /i«(t) in Eq. (60), and subtracting a diagonal piece 
((j, 2 + fil)/2E x /, we can rewrite the evolution equation as: 
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J_ / -A(t) -AiE9 m {t) \ (v? \ 

4E \ AiE9 m (t) A(t) )\ ) 1 ) 

where we defined A(t) = — /i 2 (£)- 

The evolution equations (70) constitute a system of coupled equations: the instantaneous 
mass eigenstates, z/™ - , mix in the evolution and are not energy eigenstates. The importance 
of this effect is controlled by the relative size of the off-diagonal piece 4E9 m (t) with respect 
to the diagonal one A(t). When A(t) ^> 4E9 m (t), the instantaneous mass eigenstates, z/™, 
behave approximately as energy eigenstates and they do not mix in the evolution. This is 
the adiabatic transition approximation. From the definition of 9 m in Eq. (6f ) we find: 

Am 2 sin 29 ■ . , 

*■ = ~mtT A ■ (71) 

The adiabaticity condition reads then 

2EAAm 2 sin 29 



A(t) > 



A(ty 



(72) 



Since for small mixing angles the maximum of 9 m occurs at the resonance point (as seen 
in Fig. 3), the strongest adiabaticity condition is obtained when Eq. (72) is evaluated at the 
resonance (the generalization of the condition of maximum adiabaticity violation to large 
mixings can be found in Friedland, 2001 and Lisi, Marrone, Montanino, Palazzo and Petcov, 
2001). We define the adiabaticity parameter Q at the resonance as follows: 

_ Am 2 sin 2 29 _ 4 nSr R 
^ Ecos29h R L°i c ' 1 ' 

where we used the definitions of Ar, 5rn, and Hr in Eqs. (63), (66), and (67). Written in 
this form, we see that the adiabaticity condition, Q ^> 1, implies that many oscillations take 
place in the resonant region. Conversely, when Q < 1 the transition is non-adiabatic. 

The survival amplitude of a v e produced in matter at t Q and exiting the matter at t > t 
can be written as follows: 

A(v e -> v e - 1) = £ A(u e (t ) -> Ui(t )) A(ui(to) -> Vj(t)) A(uj(t) -> u e (t)) (74) 



with 

A(u e (t ) -> Ui(t )) = (isi(t )\v e (t )) = U* ei (9 mfi ) 

where f/*j(^ m ,o) is the (ei) element of the mixing matrix in matter at the production point 
and U e j(9) is the (ej) element of the mixing matrix in vacuum. 

In the adiabatic approximation the mass eigenstates do not mix so 

A(ui(to) -> Vj(t)) = % (vi{t)\v t {t )) = Sij exp |i jf* Eiit'jdt'} . (75) 
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Note that is a function of time because the effective mass fii is a function of time, 

EM-p+^p-. (76) 
Thus the transition probability for the adiabatic case is given by 



P[y e -> i/ e ; t) 



X> ei (0R*(0 m ,o)exp 



(77) 



2£V 

For the case of two-neutrino mixing, Eq. (77) takes the form 

P[y e -> i/ e ;t) = cos 2 6> m cos 2 6> + sin 2 9 m sin 2 6> + isin26> m sin 29 cos , (78) 

2 \ 2£/ y 

where 

S(t) = f A(t')dt' = f J (Am 2 cos 29 - A(t')) 2 + (Am 2 sin 29) 2 dt' , (79) 

J to Jt 

which, in general, has to be evaluated numerically. There are some analytical approximations 
for specific forms of A(t'): exponential, linear ...(see, for instance, Kuo and Pantaleone, 
1989). For 5(t) 3> E the last term in Eq. (78) is averaged out and the survival probability 
takes the form 

P(v e -> u e ;t) = ]■ [1 + cos2^ m cos2^] (80) 

In Fig. 4 we plot isocontours of constant survival probability in the parameter plane 
(Am 2 , tan 2 9) for the particular case of the sun density for which A > 0. Notice that, 
unlike sin 2 29, tan 2 9 is a single valued function in the full parameter range < 9 < n/2. 
Therefore it is a more appropriate variable once matter effects are included and the sym- 
metry of the survival probability with respect to the change of octant for the mixing angle 
is lost. As seen in the figure, for 9 < 7r/4, P[y e — > v e ) in matter can be larger or smaller 
than 1/2, in contrast to the case of vacuum oscillations where, in the averaged regime, 
p ee c = 1 - |sin 2 2^ > \. 

In Fig. 4 we also plot the limiting curve for Q — 1. To the left and below of this curve 
the adiabatic approximation breaks down and the isocontours in Fig. 4 deviate from the 
expression in Eq. (80). In this region, the off-diagonal term 9 m cannot be neglected and the 
mixing between instantaneous mass eigenstates is important. In this case we can write 

A(ui(t ) -> Uj(t)) = (^-(*)l^(*fl)X^(*ii)l^(«fl)>(^(*ii)l^(*o)) (81) 

where tn is the point of maximum adiabaticity violation which, for small mixing angles, 
corresponds to the resonant point. The possibility of this level crossing can be described in 
terms of the Landau-Zener probability (Landau, 1932; Zener, 1932): 

Plz = \(^(t R )\u, t (t R ))\ 2 (i^j). (82) 

Introducing this transition probability in Eq. (74) we find that in the non-adiabatic regime 
(after averaging out the oscillatory term), the survival probability can be written as 
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P(i/ e -i/ e ;t) = -[! + (! 



2P LZ )cos29 m cos29] . 



(83) 



The physical interpretation of this expression is straightforward. An electron neutrino pro- 
duced at A > Ar consists of an admixture of v\ with fraction cos 2 9 m and v 2 with fraction 
sin 2 6* m . In particular, for very small mixing angles in vacuum, 9 m ~ n/2 (see Fig. 3) so v e 
is almost a pure ^2 (to) state. When the neutrino state reaches the resonance, v 2 (^1) can 
become v 2 {v\) with probability [1 — Plz] or v\ (i/ 2 ) with probability Plz- So after passing 
the resonance, the v e flux contains a fraction of v\. P el = sin 2 9 m PLz + cos 2 # m (l — Plz), 
and a fraction of z/ 2 ; -Pe2 — cos 2 9 m PLz + sin 2 # m (l — Plz)- At the exit z/ x consists of z/ e 
with fraction cos 2 9 and i/ 2 consists of u e with fraction sin 2 9 so (Parke, 1986; Haxton, 1986; 
Petcov, 1987) 



cos 9P el + sin 2 9P e2 , 



(84) 



which reproduces Eq. (83). 

The Landau-Zener probability can be evaluated in the WKB approximation (this deriva- 
tion follows Kim and Pevsner, 1993): 



P L z = exp 



-21m 



rto 
Jt R 



2E 



dt 



= exp 



-21m 



JA R 



^ 3(A)- tf(A) 



2E\A\ 



dA 



(85) 



where to (or Aq) is a root of n 2 — /i 2 =0 lying in the upper half plane of the complex variable t 
(or A), Aq = Am 2 (cos 29 + % sin 29). Assuming that \A\ is slowly varying near the resonance, 
we get 



I 

Ja r 



^ fj%(A) - ii\(A) 



lm 1 c^ll — ^ILdA 

2E\A\ 2E\A\ R 



lm 



/x 2 (A)-/i 2 (A) 



dA . 



(86) 



Shifting the integral A 

rA 



A — Am 2 cos29, we obtain: 



lm 



f r if 

/ /4(A) - /4(A) dA = lm 
Jar 1 1 Jo 



'; Am 2 sin 28 



dA\j A 2 + (Am 2 sin 29) 2 = ^(Am 2 sin2#) 2 . 

(87) 



So we finally find 



P LZ = exp 



(88) 



where Q is the adiabaticity parameter defined in Eq. (73). In this derivation we assumed that 
\A\ is almost constant. This assumption is strictly valid only for linearly growing densities 
A oc r. For other forms of A, the adiabaticity parameter Q has to be multiplied by some 
factor (for a list of some of these factors see Kuo and Pantaleone, 1989). For instance, for 
an exponential density, A ~ exp(— r), the factor is (1 — tan 2 9). Moreover this approximate 
derivation is not valid for Q <^ 1. The general form of the Landau-Zener probability for an 
exponential density can be written as (Petcov, 1988; Krastev and Petcov, 1988): 



LZ 



exp(— 7sin 2 9) — exp(— 7) 
1 - exp (-7) 



7 = 7T 



Am 2 
E\A/A\ R 



(89) 
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When v e is produced at A 3> A R and 9 is small, 9 m ~ 90°. In this case the survival 
probability is simply given by 




Since Q ~ Am 2 sin 2 29/ E, the isocontours of constant probability in this regime correspond 
to diagonal lines in the (Am 2 , tan 2 9) plane in a log-log plot, as illustrated in Fig. 4. 



E. Propagation in the Sun: MSW Effect 

As an illustration of the matter effects discussed in the previous section we describe now 
the propagation of a v e — ux neutrino system in the matter density of the sun. For the sake 
of concreteness we assume that X is some superposition of fi and r. 

The solar density distribution decreases monotonically (see Fig. 5). For R < O.9P it 
can be approximated by an exponential N e (R) = N e (0) exp (— P/r ), with r = P /1O.54 = 
6.6 x 10 7 m = 3.3 x 10 14 eV" 1 . 

After traversing this density the dominant component of the exiting neutrino state de- 
pends on the value of the mixing angle in vacuum and the relative size of Am 2 cos 29 versus 
Ao = 2 E Gf N e fi (at the neutrino production point): 

(i) Am 2 cos 29 A : matter effects are negligible and the propagation occurs as in 
vacuum. The survival probability at the sunny surface of the Earth is 

P ee ( Am 2 cos 29 > A Q ) = 1 - ^sin 2 29 > ]- . (90) 

(ii) Am 2 cos 29 > A : the neutrino does not pass the resonance but its mixing is affected 
by the matter. This effect is well described by an adiabatic propagation: 

P ee (Am 2 cos 29 > Aq) = ]- [1 + cos 29 m cos 29} . (91) 

Since the resonance is not crossed, cos 29 m has the same sign as cos 29 and the corresponding 
survival probability is also larger than 1/2. 

(hi) Am 2 cos 29 < A : the neutrino can cross the resonance on its way out. In this 
case, as discussed in the previous section, for small mixing angle in vacuum, v e ~ z/™ at the 
production point and remains v™ till the resonance point (for larger mixing but still in the 
first octant v e is a combination of v™ and z/™ with larger z/™ component). It is important in 
this case to find whether the transition is adiabatic. For the solar density, Q ~ 1 corresponds 
to 

(AmVeVW20 
(P/MeV)cos2# V ; 

For Q ^> 1 the transition is adiabatic and the neutrino state remains in the same linear 
combination of mass eigenstates after the resonance determined by 9 rn . As seen in Fig. 3, 
9 m (that is, the u e component of the state) decreases after crossing the resonance and, 
consequently, so does the survival probability P ee . In particular, for small mixing angle, v 2 
at the exit point is almost a pure u x and, consequently, P ee can be very small. Explicitly, 
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P ee (Am 2 cos 29 < A , Q > 1) = - [1 + cos 29 mfi cos 29} (93) 

can be much smaller than 1/2 because cos2^ m an d cos 26? can have opposite signs. Note 
that the smaller the mixing angle in vacuum the larger is the deficit of electron neutrinos 
in the outgoing state. This is the MSW effect (Wolfenstein, 1978; Mikheyev and Smirnov, 
1985). This behaviour is illustrated in Fig. 6 where we plot the electron survival probability 
as a function of Am 2 /E for different values of the mixing angle. 

For smaller values of Am 2 /E (right side of Fig. 6) we approach the regime where Q < 1 
and non-adiabatic effects start playing a role. In this case the state can jump from i>2 into 
v\ (or vice versa) with probability Plz- For small mixing angle, at the surface v\ ~ v e and 
the v e component of the exiting neutrino increases. This can be seen from the expression 
for P ee , 

P ee (Am 2 cos 29 < A , Q < 1) = ]- [1 + (1 - 2P LZ )cos 29 m cos 29] , (94) 

and from Fig. 6. For large mixing angles this expression is still valid but Plz has to be 
evaluated at the point of maximum adiabaticity violation which does not coincide with the 
resonant point as discussed in Friedland, (2001) and Lisi, Marrone, Montanino, Palazzo and 
Petcov, (2001). 



IV. SOLAR NEUTRINOS 

Solar neutrinos are electron neutrinos produced in the thermonuclear reactions which 
generate the solar energy. These reactions occur via two main chains, the pp-chain and the 
CNO cycle, shown in Figs. 7 and 8, respectively. There are five reactions which produce 
v e in the pp chain and three in the CNO cycle. Both chains result in the overall fusion of 
protons into 4 He: 

4p -> 4 He + 2e + + 2v e + 7, (95) 

where the energy released in the reaction, Q = 4m p — m4 He — 2m e — 26 MeV, is mostly 
radiated through the photons and only a small fraction is carried by the neutrinos, (E 2Ue ) = 
0.59 MeV. 

In order to precisely determine the rates of the different reactions in the two chains which 
would give us the final neutrino fluxes and their energy spectrum, a detailed knowledge of 
the sun and its evolution is needed. Solar Models (Bahcall and Ulrich, 1988; Bahcall and 
Pinsonneault, 1992 and 1995; Bahcall, Basu, and Pinsonneault, 1998, 2001; Turck-Chiez 
et al, 1988) describe the properties of the sun and its evolution after entering the main 
sequence. The models are based on a set of observational parameters: the surface luminosity 
(3.844 x 10 26 W), the age (4.5 x 10 9 years), the radius (6.961 x 10 s m) and the mass (1.989 x 
10 30 kg), and on several basic assumptions: spherical symmetry, hydrostatic and thermal 
equilibrium, equation of state of an ideal gas, and present surface abundances of elements 
similar to the primordial composition. Over the past four decades, the solar models have been 
steadily refined as the result of increased observational and experimental information about 
the input parameters (such as nuclear reaction rates and the surface abundances of different 
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elements), more accurate calculations of constituent quantities (such as radiative opacity 
and equation of state), the inclusion of new physical effects (such as element diffusion), and 
the development of faster computers and more precise stellar evolution codes. 

We use as Standard Solar Model (SSM) the most updated version of the model developed 
by Bahcall, Pinsonneault and Basu, (2001) (BP00). (At present, the various solar models 
give essentially the same results when interpreting the bulk of experimental data in terms 
of neutrino parameters. In this sense, the implications for particle physics discussed in this 
review are independent of the choice of solar model.) In Fig. 9 we show the energy spectrum 
of the fluxes from the five pp chain reactions. In what follows we refer to the neutrino 
fluxes by the corresponding source reaction, so, for instance, the neutrinos produced from 
8 B decay are called 8 B neutrinos. Most reactions produce a neutrino spectrum characteristic 
of j3 decay. For 8 B neutrinos the energy distribution presents deviations with respect to the 
maximum allowed energy because the final state, 8 Be, is a wide resonance. On the other 
hand, the 7 Be neutrinos are almost monochromatic, with an energy width of about 2 keV 
which is characteristic of the temperature in the core of the sun. 

As discussed in Sec. Ill, to describe the evolution of neutrinos in the solar matter, one 
needs to know other quantities that are predicted by the SSM, such as the density and 
composition of solar matter, which give the solar electron number density. As discussed 
in Sec. Ill and shown in Fig. 5, the solar matter density decreases monotonically and can 
be approximated an exponential profile. Furthermore, in order to precisely determine the 
evolution of the neutrino system one also needs to know the production point distribution 
for the different neutrino fluxes. In Fig. 10 we show the BP00 prediction for the distribution 
probability for the pp chain. 

The standard solar models have had notable successes. In particular, the comparison 
between the observed and the theoretically predicted sound speeds and pressures has tested 
the SSM and provided accurate information on the solar interior. The solar quantities that 
have been determined by helioseismology include the sound velocity and matter density 
as a function of solar radius, the depth of the convective zone, the interior rotation rate, 
and the surface helium abundance. The excellent agreement between the helioseismological 
observations and the solar model calculations has made a convincing case that the large 
discrepancies between solar neutrino measurements and solar model predictions cannot be 
due to errors in the solar models. 

A. Solar Neutrino Experiments 

1. Chlorine experiment: Homestake 

The first result on the detection of solar neutrinos was announced by Ray Davis Jr and 
his collaborators from Brookhaven in 1968 (Davis, Harmer and Hoffman, 1968). In the gold 
mine of Homestake in Lead, South Dakota, they installed a detector consisting of ~ 615 
Tons of C 2 C1 4 . Solar z/ e 's are captured via 

37 CI (u,e~) 37 Ar. 

The energy threshold for this reaction is 0.814 MeV, so the relevant fluxes are the 7 Be and 
8 B neutrinos. 7 Be neutrinos excite the Gamow- Teller (GT) transition to the ground state 
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with strength known from the lifetime of the electronic capture of 37 Ar. 8 B neutrinos can 
excite higher states of 37 Ar, including the dominant transition to the isobaric state with 
transition energy 4.99 MeV. The strength of the GT transitions to these excited states are 
determined from the isospin related (3 decay, 37 Ca(/3 + ) 37 K. The final result is that, for the 
SSM fluxes, 78% of the expected number of events are due to 8 B neutrinos while 13% arise 
from 7 Be neutrinos. The produced 37 Ar is extracted radiochemically every three months 
approximately and the number of 37 Ar decays (ti=34.8 days) is measured in a proportional 
counter. 

The average event rate measured during the more than 20 years of operation is 

.Rci = 2.56 ± 0.16 ± 0.16 SNU, 

(1 SNU = 10~ 36 captures/atom/sec) which corresponds to approximately one third of the 
SSM prediction (for the latest publications see Cleveland et al, 1998 and Lande et al, 1999). 

2. Gallium experiments: SAGE and GALLEX/GNO 

In January 1990 and May 1991, two new radiochemical experiments using a 71 Ga tar- 
get started taking data, SAGE and GALLEX. The SAGE detector is located in Baksan, 
Kaberdino-Balkaria, Russia, with 30 Tons (increased to 57 Tons from July 1991) of liquid 
metallic Ga. GALLEX is located in Gran Sasso, Italy, and consists of 30 Tons of GaCLj-HCl. 
In these experiments the solar neutrinos are captured via 

71 Ga(z/,e-) 71 Ge. 

The special properties of this target include a low threshold (0.233 MeV) and a strong 
transition to the ground level of 71 Ge, which gives a large cross section for the lower energy pp 
neutrinos. According to the SSM, approximately 54% of the events are due to pp neutrinos, 
while 26% and 11% arise from 7 Be and 8 B neutrinos, respectively. The extraction of 71 Ge 
takes places every 3-4 weeks and the number of 71 Ge decays (ti =11.4 days) is measured in 
a proportional counter. 

The event rates measured by SAGE (Abdurashitov et al, 1999; 2002) and GALLEX 
(Hampel et al, 1999) are 

#sage = 70.8T1 3 ljj;£SNU, 

#gallex = 77.5 ±6.2l£? SNU. 

while the prediction of the SSM is 130 SNU. 

The GALLEX program was completed in fall 1997 and its successor GNO started taking 
data in spring 1998. The latest combined GALLEX/GNO result is (Kirsten, 2002) 

-Rgallex+gno = 70.8 ± 5.9 SNU . 

Since the pp flux is directly constrained by the solar luminosity, in all stationary so- 
lar models there is a theoretical minimum of the expected number of events of 79 SNU. 
Furthermore, the largest uncertainties on the capture cross section for 7 Be neutrinos were 
considerably reduced after direct calibration using neutrinos from 51 Cr decay as a source. 
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3. Water Cerenkov experiments: Kamiokande and SuperKamiokande 

Kamiokande and its successor SuperKamiokande (SK) in Japan are water Cerenkov 
detectors that are able to detect in real time the electrons scattered from the water by 
elastic interaction of the solar neutrinos, 

v a + e~ -> v a + e" , (96) 

The scattered electrons produce Cerenkov light which is detected by photomultipliers. Notice 
that, while the detection process in radiochemical experiments is purely a charged-current 
(W-exchange) interaction, the detection process of Eq. (96) goes through both charged- 
and neutral-current (Z-exchange) interactions. Consequently, the detection process (96) is 
sensitive to all active neutrino flavors, although z/ e 's (which are the only ones to scatter via 
^-exchange) give a contribution that is about 6 times larger than that of v^s or z/ T 's. 

Kamiokande, with 2140 tons of water, started taking data in January 1987 and was 
terminated in February 1995. SK, with 45000 tons of water (of which 22500 are usable in 
solar neutrino measurements) started in May 1996 and it has analyzed so far the events 
corresponding to 1258 days. The detection threshold in Kamiokande was 7.5 MeV while SK 
started at a 6.5 MeV threshold and is currently running at 5 MeV. This means that these 
experiments are able to measure only the 8 B neutrinos (and the very small hep neutrino 
flux). Their results are presented in terms of measured 8 B flux. 

The final result of Kamiokande (Fukuda et al, 1996) and the latest result of SK (Smy 
et al, 2002) are 

$ K am = (2.80 ± 0.19 ± 0.33) x 10 6 cnrV 1 , 
$ SK = (2.35 ± 0.02 ± 0.08) x 10 6 cm^s" 1 , 

corresponding to about 40-50% of the SSM prediction. 

There are three features unique to the water Cerekov detectors. First, they are real time 
experiments. Each event is individually recorded. Second, for each event the scattered elec- 
tron keeps the neutrino direction within an angular interval which depends on the neutrino 
energy as \j2m e / E v . Thus, it is possible, for example, to correlate the neutrino detection 
with the position of the Sun. Third, the amount of Cerenkov light produced by the scattered 
electron allows a measurement of its energy. In summary, the experiment provides informa- 
tion on the time, direction and energy for each event. As we discuss below, signatures of 
neutrino oscillations might include distortion of the recoil electron energy spectrum, differ- 
ence between the night-time solar neutrino flux and the day-time flux, or a seasonal variation 
in the neutrino flux. Observation of these effects would be strong evidence in support of solar 
neutrino oscillations independent of absolute flux calculations. Conversely, non-observation 
of these effects can constrain oscillation solutions to the solar neutrino problem. 

Over the years the SK collaboration has presented the information on the energy and 
time dependence of their event rates in different forms. In Fig. 11 we show their spectrum, 
corresponding to 1258 days of data, relative to the predicted spectrum (Ortiz et al., 2000) 
normalized to BP00. As seen from the figure, no significant distortion is observed and the 
data is compatible with a horizontal straight line. 

The SK Collaboration has also measured the event rates as a function of the time of day 
or night or, in other words, of the position of the sun in the sky. The results are presented 
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as a nadir angle (the angle which the sun forms with the vertical) distribution of events. In 
Fig. 12 we show their nadir angle distribution, corresponding to 1258 days of data. SK has 
also presented their results on the day-night variation in the form of a day-night asymmetry, 

N — D 

An ' D " 2 D + N = °'° 21 ± °- 020 ( stat -) ± 0.013(syst.) (97) 

where D (N) is the event rate during the day (night) period. The SK results show a small 
excess of events during the night but only at the 0.8 a level. 

In order to simultaneously account for the energy and day-night variation, SK present 
the observed energy spectrum during the day and during the night separately as well as 
in several night bins. This is the most convenient observable for including both time and 
energy dependence simultaneously in the analysis. The reason for that is that the separate 
distributions - energy spectrum averaged in time and time dependence averaged in energy- 
do not correspond to independent data samples and therefore have non-trivial correlations. 

Finally, SK has also measured the seasonal dependence of the solar neutrino flux. Their 
results after 1258 days of data are shown in Fig. 13. The points represent the measured 
flux, and the curve shows the expected variation due to the orbital eccentricity of the Earth 
(assuming no new physics, and normalized to the measured total flux). As seen from the 
figure, the data are consistent with the expected annual variation. 



4. SNO 

The Sudbury Neutrino Observatory (SNO) was first proposed in 1987 and it started tak- 
ing data in November 1999 (Mcdonald et ai, 2000). The detector, a great sphere surrounded 
by photomultipliers, contains approximately 1000 Tons of heavy water, D2O, and is located 
at the Creighton mine, near Sudbury in Canada. It is immersed in ultra pure light water 
to provide shielding from radioactivity. SNO was designed to give a model independent 
test of the possible explanations of the observed deficit in the solar neutrino flux by having 
sensitivity to all flavors of active neutrinos and not just to v e . 

This sensitivity is achieved because energetic neutrinos can interact in the D2O of SNO 
via three different reactions. Electron neutrinos may interact via the Charged Current (CC) 
reaction 

u e + d^p + p + e- , (98) 

and can be detected above an energy threshold of a few MeV (presently T e > 5 MeV) . All 
active neutrinos [y a = z/ e , z/^, u T ) interact via the Neutral Current (NC) reaction 

u a + d^n + p + u a , (99) 

with an energy threshold of 2.225 MeV. The non-sterile neutrinos can also interact via Elastic 
Scattering (ES), v a + — > v a + , but with smaller cross section. 

The experimental plan of SNO consists of three phases of approximately one year dura- 
tion each. In its first year of operation, SNO has concentrated on the measurement of the 
CC reaction rate while in a following phase, after the addition of MgCl 2 salt to enhance the 
NC signal, it will also perform a precise measurement of the NC rate. In the third phase, 
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the salt will be eliminated and a network of proportional counters filled with 3 He will be 
added with the purpose of directly measuring the neutral current rate 3 He(n,p) 3 H. 

SNO can also perform measurements of the energy spectrum and time variation of the 
event rates. But the uniqueness of SNO lies in its ability to directly test if the deficit of 
solar u e is due to changes in the flavor composition of the solar neutrino beam, since the 
ratio CC/NC compares the number of v e interactions with those from all active flavors. This 
comparison is independent of the overall flux normalization. 

In June 2001, SNO published their first results (Ahmad, 2001) on the CC measurement 
and in April 2002 they published their results from the first phase of the experiment which 
include the day-night spectrum data in the full energy range above a threshold T e > 5 MeV. 
Their spectrum is the result of the combination of the three possible signals. Assuming an 
undistorted energy spectrum they extract the individual rates: 

(1.761^5 ± 0-09) x 10 6 cni'V 1 , 
(2.39 ± 0.24 ± 0.12) x 10 6 cm" V 1 , 
(5.091°;^) xlO 6 cm-V\ 

and a day-night asymmetry A N _ D = 0.07 ± 0.049(stat.) ± 0.013(syst.). 

5. Future: Borexino and Low Energy experiments 

The Borexino experiment (Oberauer, 1999) is designed to detect low-energy solar neu- 
trinos in real-time through the observation of the ES process v a + e~ — » v a + e~ . The energy 
threshold for the recoil electrons is 250 keV. It will use 300 tons of liquid scintillator in an 
unsegmented detector with 2000 photomultiplier tubes. The event rate predicted by the 
SSM for a fiducial volume of about 100 tons is about 50 events per day, mostly generated 
by the 0.86 MeV monoenergetic line of 7 Be solar neutrinos. Since this line gives a charac- 
teristic spectral signature in the ES process, the flux on the earth of 7 Be solar neutrinos 
will be determined and it will be possible to check if it is suppressed with respect to the 
one predicted by the SSM as suggested by the results of current experiments. The Borexino 
experiment is under construction in the Laboratori Nazionali del Gran Sasso in Italy and is 
scheduled to start data taking in the near future. 

A new generation of experiments aiming at a high precision real time measurement of 
the low energy solar neutrino spectrum is now under study (de Bellfon, 1999; Lanou et 
al, 1993; Arpesella, Broggini, and Cattadori, 1996; Raghavan, 1997). Some of them, such 
as HELLAZ, HERON and SUPER-MuNu, intend to detect the elastic scattering of the 
electron neutrinos with the electrons of a gas and measure the recoil electron energy and its 
direction of emission. The proposed experiment LENS plans to detect the electron neutrino 
via its absorption in a heavy nuclear target with the subsequent emission of an electron 
and a delayed gamma emission. The expected rates at these experiments for the proposed 
detector sizes are of the order of ~ 1 — 10 pp neutrinos a day. Consequently, with a running 
time of two years, they can reach a sensitivity of a few percent in the total neutrino rate at 
low energy, provided that they can achieve sufficient background rejection. 
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B. The Solar Neutrino Problem 



Table III summarizes our present knowledge of the solar neutrino fluxes, their contribu- 
tion to the expected rates and the data from measurements of solar neutrino experiments. 
The predicted event rates are linear functions of the seven important neutrino fluxes: p-p, 
pep, hep, 7 Be, 8 B, 13 N, and 15 0. We can make the following statements (see the last row of 
the table): 

• Before the NC measurement at SNO all experiments observed a flux that was smaller 
than the SSM predictions, $ obs /$ssM _ Q 3 _ Q Q 

• The deficit is not the same for the various experiments, which may indicate that the 
effect is energy dependent. 

These two statements constitute the solar neutrino problem (Bahcall, Bahcall and Shaviv, 
1968; Bahcall and Davis, 1976). More quantitatively, a fit to the data using BPOO shows 
a disagreement of 6.7 a, which means that the probability of explaining this results as a 
consequence of a statistical fluctuation of standard particle physics is P = 3 x 10~ n (we 
follow here the latest discussion by Bahcall, 2002a). 

One may wonder about the solar model dependence of the solar neutrino problem. To 
answer this question one can allow all the solar neutrino fluxes, with undistorted energy 
spectra, to be free parameters in fitting the measured solar neutrino event rates, subject 
only to the condition that the total observed luminosity of the sun is produced by nuclear 
fusion. This luminosity constraint can be written as a convenient linear equation in the 
neutrino fluxes: 



where fa are the individual neutrino fluxes {i = pp,pep, hep, 7 Be, 8 B, 13 N, and 15 0) and a, is 
the luminous energy released in the corresponding reaction (the most recent determination of 
these parameters can be found in Bahcall, 2002b). Since there are at present six experiments, 
one cannot use all seven of the neutrino fluxes as free parameters and some additional 
constraints must be used. For instance, the ratio of the pep to pp fluxes is usually taken to 
be the same as in the SSM. Also the CNO nuclear reactions are assumed to be in equilibrium 
and the 13 N and 15 neutrino fluxes are taken to be exactly equal (or, in some analyses, 
both are taken to be zero). Finally the contribution from hep neutrinos is usually neglected. 

A fit of this type (Bahcall, 2002b) shows that, previous to the SNO measurement, the 
best fit corresponded always to an unphysical, negative 7 Be flux. Forcing all fluxes to be 
positive, the hypothesis of no-new-physics was rejected at the effective 2.5a level (99% C.L.). 
Furthermore, it required that 7 Be neutrinos be entirely missing, a result which many authors 
have argued is not physically or astrophysically reasonable (for instance, Bahcall and Bethe, 
1990; Hata, Bludman, and Langacker, 1994; Parke, 1995; Heeger and Robertson, 1996; 
Bahcall, Krastev and Smirnov, 1998). 

The results of SNO have provided further model independent evidence of the problem. 
Both SNO and SK are sensitive mainly to the 8 B flux. Without NP, the measured fluxes in 
any reaction at these two experiments should be equal. Villante, Fiorentini and Lisi (1999) 
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showed how one can choose the energy thresholds for the SK and SNO experiments in such 
a way that the response functions for the two experiments are made approximately equal 
(see also Fogli et al, 2001). The advantage of this method is that some of the systematic 
errors are reduced, but there is a slight loss of statistical power (Bahcall, 2002a, Berezinski, 
2001). The first reported SNO CC result compared with the ES rate from SK showed that 
the hypothesis of no flavour conversion was excluded at ~ 3c. 
Finally, with the NC measurement at SNO one finds that, 

= nc,sno _ cc,sno = ^ ± q.45^) x 10 6 mTV 1 . (101) 

This result provides evidence for neutrino flavor transition (from v e to v^t) at the level of 
5.3o". This evidence is independent of the solar model. 

C. Solar Neutrino Oscillation Probabilities 

The most generic and popular explanation of the solar neutrino anomaly is given by the 
introduction of neutrino masses and mixing, leading to oscillations of v e into an active [y^ 
and/or v T ) or a sterile (u s ) neutrino. We now discuss some issues that have been raised in the 
literature concerning the computation of the corresponding neutrino survival probabilities 
in the full range of mass and mixing relevant to the solar neutrino problem. 

1. Quasivacuum oscillations and the Dark Side 

The presence of neutrino mass and mixing implies the possibility of neutrino oscillations. 
Solar electron neutrinos can undergo oscillations either in vacuum or via the matter-enhanced 
MSW mechanism, depending on the actual values of mass-squared differences and mixing 
angles. However, this simplified picture of solar neutrino oscillations contains a set of ap- 
proximations which are not always valid in the context of solar neutrinos. In order to clarify 
this issue let us first review the calculation of the solar neutrino survival probability in the 
two-neutrino case. 

The survival amplitude for a solar v e of energy E at a terrestrial detector can be written 
as follows: 

A ee = j2A s ei Af e exp[- im ^L-r)/(2E)] . (102) 

Here is the amplitude of the transition v e — > z/« from the production point to the surface 
of the Sun, Af e is the amplitude of the transition Vi — > v e from the surface of the Earth to the 
detector. The propagation amplitude in vacuum from the surface of the Sun to the surface 
of the Earth is given by the exponential term, where L is the distance between the center 
of the Sun and the surface of the Earth, and r is the radius of the Sun. The corresponding 
survival probability P ee is then given by: 

P ee = P 1 P le + P 2 P 2e + 2 y/P 1 P 2 P u P 2e cos^. (103) 
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Here Pj = |^4fj| 2 is the probability that the solar neutrinos reach the surface of the Sun as 
\ui), while Pi e = |A^,| 2 is the probability of z/j arriving at the surface of the Earth to be 
detected as a u e . Unitarity implies Pi + P 2 = 1 and P le + P 2e = 1. The phase £ is given by 

t = ^|z0 + ,, (104) 

where <5 contains the phases due to propagation in the Sun and in the Earth and can be safely 
neglected. In the evaluation of both Pj and Pj e the effect of coherent forward interaction 
with the matter of the Sun and the Earth must be taken into account. 

From Eq. (103) one can recover more familiar expressions for P ee that hold in certain 
limits: 

(1) For Am 2 /E < 5 x 10~ 17 eV, the matter effect suppresses flavor transitions in both the 
Sun and the Earth. Consequently, the probabilities Pi and P 2e are simply the projections 
of the v e state onto the mass eigenstates: Pi = cos 2 9, P 2e = sin 2 9. In this case we are left 
with the standard vacuum oscillation formula, 

P™ c = 1 - sin 2 29 sin 2 [Am 2 (L - r)/4P], (105) 

which describes the oscillations on the way from the surface of the Sun to the surface of the 
Earth. The probability is symmetric under the change of octant, 9 <-> f — 9, and change 
of sign Am 2 <-> —Am 2 (see Sec. III). This symmetry implies that each point in the (Am 2 , 
sin 2 29) parameter space corresponds to two physically independent solutions, one in each 
octant. 

Averaging Eq. (105) over the Earth Orbit, L(t) = L [l — ecos 2 ^], one gets : 



(PIT) =l~\ sin 2 29 



Am 2 L \ T f eAm 2 L 



1 — cos — — — — J, 



2E r \ 2E 



(106) 



where e = 0.0167 is the orbit eccentricity and J is the Bessel function. In Fig. 14. a we 
display the value of (P™°) as a function of AE/ Am 2 . As seen in the figure, for large values 
of Am 2 the probability averages out to a constant value 1 — \ sin 2 {29). 

(2) For Am 2 /E > 10~ 14 eV, the last term in Eq. (103) vanishes and we recover the 
incoherent MSW survival probability. In this case Pi and P 2e must be obtained by solving 
the evolution equation of the neutrino states in the Sun and the Earth respectively. As 
discussed in Sec. III.D, the approximate solution for the evolution in the Sun takes the 
well-known form 

Pi = g + (I ~ p lz) cos 29 mfi . (107) 

Here Plz denotes the standard Landau-Zener probability of Eq. (82) which, for the exponen- 
tial profile, takes the form shown in Eq. (89), and 9 m $ is the mixing angle in matter at the 
neutrino production point given in Eq. (61). For the approximation of exponential density 
profile, 7 = 7r-^-r , which is independent of the point in the Sun where the resonance takes 
place. Improvement over this constant-slope exponential density profile approximation can 
be obtained by numerically deriving the exact N e (r) profile at the resonant point. In this 
case 7 = 7r^f-r (r res ). 
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During the day there is no Earth matter effect. The survival probabilities P{ e are obtained 
by simple projection of the v e state onto the mass eigenstates: P2 e ,D — 1 — Pie,D — sin 2 #. 
One obtains the expression in Eq. (83) 

P^d =l+(\- P lz) cos 26 mfi cos 26. (108) 

In Fig. 14. b (14. c) we plot the this survival probability as a function of AE/Am 2 for a large 
(small) mixing angle. 

Let us make some remarks concerning Eq. (108): 

(i) In both the limits of large and very small E/Am 2 , P^ sw -> \-\ sin 2 26 (see Fig. 14.b) 
which is the same expression as for averaged vacuum oscillations. This result, however, 
comes from very different reasons in the two regimes. For large E/Am 2 , Plz = cos 2 6 and 
cos2# m) o = — 1. For small E/Am 2 , Plz = and cos2# mj o = cos 2$. 

(ii) Due to matter effects, -P e ^ sw is only symmetric under the simultaneous transforma- 
tion (Am 2 , 6) — > (—Am 2 , | — 9). For Am 2 > 0, the resonance is only possible for 6 < j and 
MSW solutions are usually plotted in the (Am 2 , sin 2 26) plane assuming that now each point 
on this plane represents only one physical solution with 6 in the first octant. But in principle 
non- resonant solutions are also possible for 6 > |, the so called dark side (Fogli, Lisi and 
Montanino, 1994 and 1996; de Gouvea, Friedland and Murayama, 2000; Gonzalez-Garcia 
and Pena-Garay, 2000; Fogli, Lisi, Montanino and Palazzo, 2000a). 

(3) The analysis of the survival probability of solar neutrinos would be simplified if there 
were a region where E/Am 2 is small enough that vacuum oscillations are averaged out and 
large enough to be in the extremely non-adiabatic MSW regime. This, however, is not the 
case, as can be seen form comparison of panels (a) and (b) of Fig. 14. In the intermediate 
range, 2 x 10 14 < AE/Am 2 < 10 16 eV _1 , adiabaticity is violated and the cos£ coherent term 
should be taken into account. The result is similar to vacuum oscillations but with small 
matter corrections. We define this case as quasi-vacuum oscillations (QVO) (Petcov, 1988b; 
Petcov and Rich, 1989; Pantaleone, 1990; Pakvasa and Pantaleone, 1990; Friedland 2000; 
Fogli, Lisi, Montanino, and Palazzo 2000b; Lisi, Marrone, Montanino, Palazzo and Petcov, 
2001). The range of E/Am 2 for the QVO regime depends on the value of E/Am 2 for which 
the MSW probability in Eq. (108) acquires the asymptotic value 1 — \ sin 2 26: the smaller 
E/Am 2 the more separated the MSW and vacuum regimes are, and the narrower the QVO 
region is. 

It is clear from these considerations that in order to compute the survival probability for 
solar neutrinos that would be valid for any experiment and any value of the neutrino mass 
and mixing, the full expression (103) has to be evaluated. The results that we show here 
were obtained using the general expression for the survival probability in Eq. (103) with Pi 
and P2e found by numerically solving the evolution equation in the Sun and the Earth. In 
the Sun the evolution equation has to be solved from the neutrino production point to the 
edge of the Sun and averaged over the production point distribution shown in Fig. 10. 

2. Evolution in the Earth 

In this section we discuss matter effects in the Earth (for some of the original works see 
J. Bouchez et. al. 1986, Mikheyev and Smirnov, 1987, A. J. Baltz and J. Weneser 1987 and 
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1988). Unlike the Sun, the Earth's matter density is not strongly varying. It consists mainly 
of two layers, the lower- density mantle and the higher-density core, each with approximately 
constant density. The core radius is L = 2896 km. 

When a neutrino crosses only the mantle, the Earth matter effects are well approximated 
by evolution with a constant potential. In this case, the resulting probability simplifies to: 

P 2e = sin 2 9 + ^4 sin 2 28 E sin 2 — , (109) 
Am 2 L rn 

where $e is the mixing angle in the Earth [Eq. (61)], L is the distance traveled by the 
neutrino within the Earth and L m is the oscillation length in matter [Eq (64)]. 

The right hand side of Eq. (109) consists of two terms. The first term gives a simple 
projection of the mass to flavor state in vacuum (which corresponds to the probability during 
the day time). The second is the regeneration term and will be denoted by / reg (for details 
see Gonzalez-Garcia, Pena-Garay and Smirnov, 2001). It contains the Earth matter effects 
and it is always positive. Averaging out the distance dependent term, we can now rewrite 
Eq. (103): 

Pee = Pee ~ (1 ~ ^Lz) COS 26 mfi f reg , (110) 

where (P^) is the v e survival probability during the day (night). Since / rcg > and, for 
the interesting cases of MSW transitions in the Sun, cos 2# mi o < 0, the relative magnitude 
of P® e and P^ e depends on the sign of (1 — IPlz)- For Plz < 0.5 (and, in particular, for 
adiabatic MSW transitions, P LZ = 0) the survival probability during the night is larger than 
during the day. The opposite holds for large adiabaticity violations, Plz > 0.5. 

When neutrinos cross the core, adiabaticity may be violated in their evolution. The 
abrupt density change between the mantle and the core may induce a new form of reso- 
nance, different from the MSW resonance (Petcov, 1998; Akhmedov, 1999) This parametric 
resonance (Ermilova, Tsarev and Chechin, 1986; Akhmedov, 1987) is relevant mainly for 
small mixing angles. 



D. Two-neutrino Oscillation Analysis 

1. Predictions 

The expected event rate in the presence of oscillations in the experiment i, Rf 1 , can be 
written as follows: 

Rf 1 = E & jdE v \ k {E v ) x [a e ,{E v ){P ee {E v ,t)) + a x ,(E u )(l - (P ee (E v ,t)))\, (111) 

fc=l,8 

where E v is the neutrino energy, 0^ is the total neutrino flux and is the neutrino energy 
spectrum (normalized to 1) from the solar nuclear reaction k. o~ e ^ (o~ Xt i) is the v e (u x , 
x = fi, r) interaction cross section in the SM with the target corresponding to experiment 
i, and (P ee (E v ,t)) is the time-averaged u e survival probability. The expected signal in the 
absence of oscillations, Rf P0 °, can be obtained from Eq. (Ill) by substituting P ee = 1. In 
Table III we give the expected rates at the different experiments which we obtain using the 
fluxes of Bahcall, Pinsonneault, and Basu (2001). 
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For the Chlorine, Gallium and SNO(CC) measurements, only the electron neutrino con- 
tributes and the cr^-term in Eq. (Ill) vanishes. For ES at SK or SNO there is a possible 
contribution from the NC interaction of the other active neutrino flavors present in the beam. 
For the NC rate at SNO all active flavours contribute equally. For absorption in chlorine 
and gallium, the cross sections a eji (E) can be found in Bahcall (1997). The cross sections 
for SNO and SK(ES) can be obtained from the differential cross sections of (Nakamura et 
al., 2001) and (Bahcall, Kamionkowsky, and Sirlin, 1995), respectively, by integrating with 
the corresponding detector resolutions and for the given detection thresholds. In particular, 
for the SK and SNO energy spectrum data, to obtain R th in a given energy bin one has 
to integrate within the corresponding electron recoil energy bin and to take into account 
that the finite energy resolution implies that the measured kinetic energy T of the scattered 
electron is distributed around the true kinetic energy T 1 according to a resolution function 
Res(T, T") of the form (from Bahcall, Krastev and Lisi, 1997): 



Res(T, T') = - 7 ^exp 



(T — T 
2^ 



l\2 



;ii2) 



where s = s y/T'/M.eV, and s = 0.47 MeV for SK (from Faid et al, 1997). On the other 
hand, the distribution of the true kinetic energy T' for an interacting neutrino of energy E v 
is dictated by the differential cross section da a (E u , T')/dT' and the kinematic limits are: 

0<T'<T'(E l/ )= f u . (113) 
1 + m e /2E u 

For assigned values of sq, T min , and T max , the corrected cross section a a (E) (a = e, x) is 
given as 



/■■j max (T\e v ) da JR. T') 
a a (E u ) = dT dT' Res(T, T') ] . (114) 

Jr min Jo dl 

For data taken within a given zenith angle period i, the expected number of events in 
the presence of oscillations is 

1 /"r(cos<I> maXi i) r 

Rf = — drJ2 ^ dE u \ k {E u ) x [a e ,(E u )(P ee (E u ,r)} (115) 

ZATj JT(cos$> m[ ^i) k=18 

+<T X7i (E„)(l-(P ee (E v ,T)))], 

where r measures the annual averaged length of the period % normalized to 1: Arj = 
r(cos $ max ,i) _ r(cos $ m i n ,i) =• For instance, if the period i corresponds to the entire day 
(i = D) or night (i = N) then At D ( N ) = 0.5. From these predictions one can easily build 
the corresponding expected day-night asymmetry. In Fig. 15 we show the isocontours of 
expected A^-d at SK in the (Am 2 , tan 2 6) plane for active neutrino oscillations (the results 
are very similar for sterile neutrinos). As discussed in Sec. IV. C, in most of the parameter 
space, the effect of traveling across the Earth is the regeneration of the v e component, 
resulting in a positive day-night asymmetry with the exception of the region where the non- 
adiabaticity of the oscillations in the Sun is important (tan 2 9 ~ 10~ 3 , Am 2 ~3x 10~ 6 eV 2 ). 
As discussed above, SK has observed a very small day-night asymmetry. This observation 
implies that some of the oscillation parameters space can be excluded as we discuss below. 
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In the same fashion, integrating Eq. (115) for the different electron recoil energy bins 
and for the day and night periods separately one can obtain the predictions for SK and 
SNO day-night spectrum. 

2. Analysis of total event rates: allowed masses and mixing 

The goal of the analysis of the solar neutrino data in terms of neutrino oscillation is to 
determine which range of mass-squared difference and mixing angle can be responsible for the 
observed deficit (see for instance Hata and Langacker, 1997; Fogli, Lisi and Montanino, 1998; 
Bahcall, Krastev and Smirnov, 1998; Gonzalez-Garcia et al, 2000). In order to answer this 
question in a statistically meaningful way one must compare the predictions in the different 
oscillation regimes with the observations, including all the sources of uncertainties and their 
correlations, by building, for instance, the x 2 function 

xl = E(*f - R?*W(RT - R e D- ( 116 ) 

i 

Here Rf 1 is the theoretical expectation (111), Rf p is the observed number of events in the 
experiment i, and is the error matrix which contains both the theoretical uncertainties 
and the experimental statistical and systematic errors. The main sources of uncertainty are 
the theoretical errors in the prediction of the solar neutrino fluxes for the different reactions. 
These errors are due to uncertainties in the twelve basic ingredients of the solar model, which 
include the nuclear reaction rates (parametrized in terms of the astrophysical factors Sn, 
S33, S34, 5*1,14 and 6*17), the solar luminosity, the metalicity Z/X, the sun age, the opacity, 
the diffusion, and the electronic capture of 7 Be, Cse- These errors are strongly correlated 
as the same astrophysical factor can affect the different neutrino production rates. Another 
source of theoretical error arises from the uncertainties in the neutrino interaction cross 
section for the different detection processes. For a detailed description of the way to include 
all these uncertainties and correlations in the construction of we refer to the work of 
Fogli and Lisi (1995a). Updated uncertainties can be found in Fogli et al. (2000a) and 
Gonzalez-Garcia and Pena-Garay (2001). 

The results of the analysis of the total event rates are shown in Fig. 16 where we plot the 
allowed regions which correspond to 90%, 95%, 99% and 99.73% (3a) CL for v e oscillations 
into active neutrinos. 

As seen in the figure, for oscillations into active neutrinos there are several oscillation 
regimes which are compatible within errors with the experimental data. These allowed 
parameter regions are denoted as MSW small mixing angle (SMA), MSW large mixing 
angle (LMA), MSW low mass (LOW) and vacuum oscillations (VAC). Before including 
the SNO(CC) data, the best fit corresponded to the SMA solution, but after SNO the 
best fit corresponds to the LMA solution. For the LMA solution, oscillations for the 8 B 
neutrinos occur in the adiabatic regime and the survival probability is higher for lower 
energy neutrinos. This situation fits well the higher rate observed at gallium experiments. 
For the LOW solution, the situation is opposite but matter effects in the Earth for pp and 7 Be 
neutrinos enhance the average annual survival probability for these lower energy neutrinos. 
The combination of these effects still allows a reasonable description of the Gallium rate. 
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Oscillations into pure sterile neutrinos are strongly disfavour by the SNO data since they 
would imply very similar neutrino fluxes for the NC, CC and ES rates in SNO as well as the 
ES rate at SK. Schematically, in presence of oscillations, 

$ CC = $ e , 

$ ES = $ e + r^, (117) 
$nc = $e + ^ 

where r = o^j o e ~ 0.15 is the ratio of the the v e — e and z/ M — e elastic scattering cross-sections. 
The flux $ MT of active no-electron neutrinos in the beam would vanish if v e oscillate into a 
purely sterile state. Thus, if the beam comprises of only z/ e 's and z/ s 's, the three observed 
rates should be equal (up to effects due to spectral distortions) , an hypothesis which is now 
ruled out at more than 5-sigma by the SNO data. Oscillations into an admixture of active 
and sterile states are still allowed provided that the 8 B neutrino flux is allowed to be larger 
than the SSM expectation (Barger, Marfatia, and Whisnant, 2002; Bahcall, Gonzalez-Garcia 
and Peha-Garay, 2002b). 



3. Day-Night spectra: excluded masses and mixing 

Further information on the different oscillation regimes can be obtained from the analysis 
of the energy and time dependence data from SK (for an early suggestion of the use of the 
day-night spectra see Maris and Petcov, 1997), which is currently presented in the form of 
the observed day-night or zenith spectrum (Fukuda et al, 2001: Smy et al, 2002). The 
way to statistically treat these data has been discussed for instance in Fukuda et al. (2001); 
Bahcall, Krastev, and Smirnov (2000); Gonzalez-Garcia et al. (2000); Gonzalez-Garcia and 
Peha-Garay (2001). 

The observed day-night spectrum in SK is essentially undistorted in comparison to the 
SSM expectation and shows no significant differences between the day and the night periods. 
Consequently, a large region of the oscillation parameter space where these variations are 
expected to be large can be excluded. In Fig. 11 we show the SK spectrum corresponding to 
1258 days of data relative to the Ortiz et al. (2000) spectrum normalized to BP00, together 
with the expectations from the best fit points for the LMA, SMA, LOW and VAC solutions. 

The various solutions give different predictions for the day-night spectrum For LMA 
and LOW, the expected spectrum is very little distorted. For SMA, a positive slope is 
expected, with larger slope for larger mixing angle within SMA. For VAC, large distortions 
are expected. The details are dependent on the precise values of the oscillation parameters. 

In Fig. 17 we show the excluded regions at the 99% CL from the analysis of SK day-night 
spectrum data, together with the contours corresponding to the 95% and 99.73% (3a) CL. 
In particular, the central region (2 x 10~ 5 < Am 2 < 3 x 10~ 7 , tan 2 9 > 3 x 10~ 3 ) is excluded 
due to the small observed day-night variation (compare to Fig. 15). The rest of the excluded 
region is due to the absence of any observed distortion of the energy spectrum. 

Superimposing Fig. 17 and Fig. 16 we can deduce the main consequences of adding the 
day-night spectrum information to the analysis of the total event rates: 

• SMA: within this region, the part with larger mixing angle fails to comply with the 
observed energy spectrum, while the part with smaller mixing angles gives a bad fit 
to the total rates. 
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• VAC: the observed flat spectrum cannot be accommodated. 

• LMA and LOW: the small Am 2 part of LMA and the large Am 2 part of LOW are 
reduced because they predict a day-night variation that is larger than observed. Both 
active LMA and active LOW solutions predict a flat spectrum in agreement with the 
observation. 

4. Global analysis 

In order to quantitatively discuss the combined analysis of the full bulk of solar neutrino 
data, one must define a global statistical function. How to implement such a program is a 
question that is being discussed intensively in the literature (Fogli, Lisi, Marrone, Montanino 
and Palazzo 2002; Bahcall, Gonzalez-Garcia, and Pena-Garay, 2001a, 2002a, 2002c; Garzelli 
and Giunti, 2001; de Holanda and Smirnov, 2002; Creminelli, Segnorelli and Strumia, 2001; 
Barger, Marfatia, Whisnant and Wood, 2002; Maltoni, Tortola, Schwetz, Valle 2002). The 
situation at present is that, although one finds slightly different globally allowed regions, 
depending on the particular prescription used in the combination, the main conclusions are 
independent of the details of the analysis. 

We show in Fig. 18 and Table IV the results from the global analysis of Bahcall, Gonzalez- 
Garcia and Pena-Garay, (2002c). In Fig. 18 we show the allowed regions which correspond 
to 90%, 95%, 99% and 99.73% (3a) CL for v e oscillations into active neutrinos. In the table 
we list the local minima of the allowed regions, the value of Xmin lYy eac h local minimum. 

The results show that at present the most favoured solution is the LMA oscillation 
while the LOW solution provides a worse fit. There are some small allowed islands for 
VAC oscillations. The active SMA solution does not appear at 3a as a consequence of the 
incompatibility between the observed small CC rate at SNO, which would favour larger 
mixing, and the flat spectrum which prefers smaller mixing. Oscillations of solar neutrinos 
into pure sterile state are also disfavoured at the 5.4<r level. As explained below, in Sec. VII, 
this difference between active and sterile oscillations has important implications on 4-z/ 
mixing schemes. 

V. ATMOSPHERIC NEUTRINOS 

Atmospheric neutrinos are produced in cascades initiated by collisions of cosmic rays 
with the Earth's atmosphere. Some of the mesons produced in these cascades, mostly 
pions and some kaons, decay into electron and muon neutrinos and anti-neutrinos. The 
predicted absolute fluxes of neutrinos produced by cosmic-ray interactions in the atmosphere 
are uncertain at the 20% level. The ratios of neutrinos of different flavor are, however, 
expected to be accurate to better than 5%. Since v e is produced mainly from the decay 
chain n — > fiu^ followed by fj, — > eu^u ei one naively expects a 2 : 1 ratio of to u e . (For 
higher energy events the expected ratio is smaller because some of the muons arrive to 
Earth before they had time to decay.) In practice, however, the theoretical calculation of 
the ratio of muon-like interactions to electron-like interactions in each experiment is more 
complicated. 
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Atmospheric neutrinos are observed in underground experiments using different tech- 
niques and leading to different type of events, which we briefly summarize here. They can 
be detected by the direct observation of their charged current interaction inside the detec- 
tor. These are the contained events. Contained events can be further classified into fully 
contained events, when the charged lepton (either electron or muon) that is produced in 
the neutrino interaction does not escape the detector, and partially contained muons, when 
the produced muon exits the detector. For fully contained events the flavor, kinetic energy 
and direction of the charged lepton can be best determined. As discussed later, some ex- 
periments further divide the contained data sample into sub-GeV and multi-GeV events, 
according to whether the visible energy is below or above 1.2 GeV. On average, sub-GeV 
events arise from neutrinos of several hundreds of MeV while multi-GeV events are origi- 
nated by neutrinos with energies of the order of several GeV. Higher energy muon neutrinos 
and antineutrinos can also be detected indirectly by observing the muons produced in their 
charged current interactions in the vicinity of the detector. These are the so called upgoing 
muons. Should the muon stop inside the detector, it is classified as a stopping muon (which 
arises from neutrinos E u ~ 10 GeV), while if the muon track crosses the full detector the 
event is classified as a through-going muon (which is originated by neutrinos with energies 
of the order of hundreds of GeV). Downgoing muons from interactions above the detec- 
tor cannot be distinguished from the background of cosmic ray muons. Higher energy z/ e 's 
cannot be detected this way as the produced e showers immediately in the rock. In Fig. 19 
we display the characteristic neutrino energy distribution for these different type of events 
(taken from Engel, Gaisser and Stanev, 2000). 

Atmospheric neutrinos were first detected in the 1960's by the underground experiments 
in South Africa (Reines et al, 1965) and the Kolar Gold Field experiment in India (Achar 
et al, 1965). These experiments measured the flux of horizontal muons (they could not 
discriminate between downgoing and upgoing directions) and although the observed total 
rate was not in full agreement with theoretical predictions (Zatsepin and Kuzmin, 1962; 
Cowsik et al, 1965; Osborne, Said and Wolfendale, 1965), the effect was not statistically 
significant. 

A set of modern experiments were proposed and built in the 1970's and 1980's. The 
original purpose was to search for nucleon decay, for which atmospheric neutrinos constitute 
background, although the possibility of using them to search for oscillation was also known 
(Ayres et al, 1984). Two different detection techniques were employed. In water Cerenkov 
detectors the target is a large volume of water surrounded by photomultipliers which detect 
the Cerenkov-ring produced by the charged leptons. The event is classified as an electron-like 
(muon-like) event if the ring is diffuse (sharp). In iron calorimeters, the detector is composed 
of a set of alternating layers of iron which act as a target and some tracking element (such as 
plastic drift tubes) which allows the reconstruction of the shower produced by the electrons 
or the tracks produced by muons. Both types of detectors allow for flavor classification of 
the events. 

The two oldest iron calorimeter experiments, Frejus (Daum et al, 1995) and NUSEX 
(Aglietta et al, 1989), found atmospheric neutrino fluxes in agreement with the theoretical 
predictions. On the other hand, two water Cerenkov detectors, 1MB (Becker-Szendy et al, 
1992) and Kamiokande, detected a ratio of z/^-induced events to z/ e -induced events smaller 
than the expected one by a factor of about 0.6. Kamiokande performed separate analyses 
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for both sub-GeV neutrinos (Hirata et al., 1992) and multi-GeV neutrinos (Fukuda et al, 
1994), which showed the same deficit. This was the original formulation of the atmospheric 
neutrino anomaly. In Fig. 20 we show the values of the ratio R^/ e / Rff/e 5 which denotes the 
double ratio of experimental-to-expected ratio of muon-like to electron-like events. Whether 
R^/e/R^e is small because there is z/ M disappearance or v e appearance or a combination 
of both could not be determined. Furthermore, the fact that the anomaly appeared only 
in water Cerenkov and not in iron calorimeters left the window open for the suspicion of a 
possible systematic problem as the origin of the effect. 

Kamiokande also presented the zenith angular dependence of the deficit for the multi- 
GeV neutrinos (Fukuda et al, 1994). The zenith angle, parametrized in terms of cos#, 
measures the direction of the reconstructed charged lepton with respect to the vertical 
of the detector. Vertically downgoing (upgoing) particles correspond to cos 6 = +1(— 1). 
Horizontally arriving particles come at cos 9 = 0. Kamiokande results seemed to indicate 
that the deficit was mainly due to the neutrinos coming from below the horizon. Atmospheric 
neutrinos are produced isotropically at a distance of about 15 km above the surface of the 
Earth. Therefore neutrinos coming from the top of the detector have traveled approximately 
those 15 kilometers before interacting while those coming from the bottom of the detector 
have traversed the full diameter of the Earth, ~ 10 4 Km before reaching the detector. The 
Kamiokande distribution suggested that the deficit increases with the distance between the 
neutrino production and interaction points. 

In the last five years, the case for the atmospheric neutrino anomaly has become much 
stronger with the high precision and large statistics data from SK and it has received im- 
portant confirmation from the iron calorimeter detectors Soudan2 and MACRO. 

SK is a 50 kiloton water Cerenkov detector constructed under Mt. Ikenoyama located 
in the central part of Japan, giving it a rock over-burden of 2,700 m water-equivalent. The 
fiducial mass of the detector for atmospheric neutrino analysis is 22.5 kilotons. In June 
1998, in the Neutrino98 conference, SK presented evidence of z/ M oscillations (Fukuda et al, 
1998) based on the angular distribution for their contained event data sample. Since then 
SK accumulated more statistics and has also studied the angular dependence of the upgoing 
muon sample (Fukuda et al, 1999a and 1999b). In Fig. 21 (from Toshito et al, 2001) we 
show their data at summer 2001, corresponding to 79 kiloton year (1289 days) exposure 
for their contained sub-GeV (2864 1-ring e-like events and 2788 1-ring /i-like events) and 
multi-GeV (626 1-ring e-like events, 558 1-ring /i-like events and 754 PC events), as well as 
the stopping-muons and through-going muon samples (1416 events). 

In the figure we show the angular zenith distribution of the different samples. Comparing 
the observed and the expected (MC) distributions, we can make the following statements: 

(i) v e distributions are well described by the MC while z/ M presents a deficit. Thus the 
atmospheric neutrino deficit is mainly due to disappearance of and not the appearance 
of v e . 

(ii) The suppression of contained /i-like events is stronger for larger cos 6, which implies 
that the deficit grows with the distance traveled by the neutrino from its production point 
to the detector. This effect is more obvious for multi-GeV events because at higher energy 
the direction of the charged lepton is more aligned with the direction of the neutrino. It can 
also be described in terms of an up-down asymmetry: 
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U-D 
U + D 



0.316 ± 0.042(stat.) ± 0.005(syst.) 



(118) 



where U (D) are the contained /i-like events with zenith angle in the range — 1 < cos 9 < —0.2 
(0.2 < cos# < 1). It deviates from the SM value, = 0, by 7.5 standard deviations. 

(iii) The overall suppression of the flux of stopping-muons, $st, is by a factor of about 
0.6, similar to contained events. However, for the flux of through-going muons, $th, the 
suppression is weaker, which implies that the effect is smaller at larger neutrino energy. This 
effect is also parametrized in terms of the double flux ratio: 



which deviates from the SM value of 1 by about 3 standard deviations. 

These effects have been confirmed by the results of the iron calorimeters Soudan2 and 
MACRO which removed the suspicion that the atmospheric neutrino anomaly is simply a 
systematic effect in the water detectors. In particular Soudan2, which has mainly analyzed 
contained events (Allison et al, 1999), have measured a ratio R^/ e /R^ = 0.68 ±0.11 ±0.06, 
in good agreement with the results from the water Cerenkov experiments. The main results 
from MACRO concern the angular distribution for through-going muons (see Ambrosio et 
al. (2001) for their latest data) and shows good agreement with the results from SK. The 
Baksan experiment (Boliev et al, 1999) has also reported results on the angular distribution 
of through-going muons but their data is less precise. 

To analyze the atmospheric neutrino data in terms of oscillations one needs to have a 
good understanding of the different elements entering into the theoretical predictions of the 
event rates: the atmospheric neutrino fluxes and their interaction cross section, which we 
describe next. 

A. Atmospheric Neutrino Fluxes 

Modern calculations of atmospheric neutrino fluxes consist of a Monte Carlo procedure 
that combines the measured energy spectra and chemical composition of the cosmic ray flux 
at the top of the atmosphere with the properties of their hadronic interaction with the light 
atmospheric nuclei followed by the neutrino production from secondary n, K and \i decay. 

Present experiments use the neutrino flux calculations from mainly Honda et al. 
(1990,1995,1996), the Bartol group (Gaisser, Stanev and Barr, 1988; Barr, Gaisser and 
Stanev , 1989; Agraval et al., 1996; Lipari, Gaisser, and Stanev, 1998) and Fiorentini, Nau- 
mov and Villante (2001). These calculations have in common a one-dimensional picture, in 
which all secondary particles in the showers, neutrinos included, are considered collinear with 
the primary cosmic rays (see also Bugaev and Naumov, 1989). The three-dimensional pic- 
ture was first considered by Lee and Koh (1990) and most recently by Batistoni et al. (2000) 
and Tserknyav et al. (1999). [For detailed comparisons between the various simulations, see 
Gaisser et al. (1996), Lipari (2000) and Battistoni (2001)]. 

The flux of neutrinos of flavor % coming from the direction Q can be schematically written 

as 



$ ST /$ TH | 

obs 

$ gT /$ TH | MC 



0.635 ± 0.049(stat.) ± 0.035(syst.) ± 0.084(theo.) 



(119) 



(120) 



A 



41 



Here §a is the primary cosmic-ray spectrum, Ra is the geomagnetic cutoff for the protons 
or light nuclei incident on the atmosphere from the direction Q, and Y PjTl ^ Vi is the yield 
per nucleon of v { . The separation to different nuclear species of the cosmic ray spectrum 
is necessary because the neutrino yield depends on the energy-per-nucleon of the incident 
cosmic rays but the geomagnetic cutoff depends on the magnetic rigidity (1Z = pc/Ze) which, 
at the same energy per nucleon, depends on A/Z. 

We now briefly summarize the uncertainties in each of these three ingredients of the 
calculation, namely the primary cosmic ray flux, the geomagnetic effects and the hadronic 
interactions on light nuclei. 

1. Cosmic ray spectrum 

The cosmic radiation incident at the top of the atmosphere includes all stable charged 
particles and nuclei. Apart from particles associated with solar flares, the cosmic rays are 
assumed to come from outside the solar system. In Fig. 22 (from Simpson, 1983) we show 
a compilation of the major components of the cosmic rays as a function of the energy per 
nucleon. The following features of the cosmic ray spectrum are relevant to the atmospheric 
neutrino fluxes: 

(i) Composition: Most of the primaries are protons although the chemical composition 
varies with energy at low energies. Above a few GeV of energy-per-nucleon, about 79% of 
the primaries are protons while helium, the next most abundant component, is about 3%. 
This fraction remains constant till beyond 100 TeV. 

(ii) Energy Dependence: For energy above a few GeV, the cosmic ray flux is a a steeply 
falling function of the energy: dip/dE oc E~ 3J . 

(iii) Solar modulation: At energies below ~ 10 GeV, the incoming charged particles are 
modulated by the solar wind which prevents the lower energy cosmic rays from the inner 
solar system to reach the Earth. Because of this effect there is an anticorrelation between 
the eleven-year-cycle solar activity and the intensity of cosmic rays below 10 GeV, which 
makes the cosmic ray fluxes time dependent. The flux difference at solar maximum and 
solar minimum is more than a factor of two for 1 GeV cosmic rays, and it decreases to ~ 
10 % for 10 GeV cosmic rays. The fluxes shown in Fig. 22 correspond to a particular epoch 
of the solar cycle. 

2. Geomagnetic effects 

Lower energy cosmic rays are affected by the geomagnetic field which they must penetrate 
to reach the top of the atmosphere. This effect depends on the point of the Earth where the 
cosmic rays arrive (being stronger near the geomagnetic equator), and on their direction. 
Near the geomagnetic poles almost all primary particles can reach the atmosphere moving 
along field lines. In contrast, close to the geomagnetic equator the field restricts the flux 
at the top of the atmosphere to particles with energy larger than a few GeV, the exact 
value depending on the direction of the particle trajectory. The relevant quantity is the 
magnetic rigidity, 1Z = pc/Ze, which parameterizes the characteristic radius of curvature of 
the particle trajectory in the presence of a magnetic field: particles with TZ below a local 
cutoff are bent away by the Earth magnetic field and do not reach the atmosphere. 
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Given a map of the magnetic field around the Earth, the value of the rigidity cutoff 
can be obtained from a computer simulation of cosmic ray trajectories using a backtracking 
technique in which one determines if a given primary particle three-momentum and position 
correspond to an allowed or a forbidden trajectory. To do so, one integrates back the equation 
of motion of the cosmic ray particle (actually one integrates the anti-particle equation) in 
the geomagnetic field, and finds if the past trajectory remains confined to a finite distance 
from the Earth, in which case it is a forbidden trajectory, or it originates from large enough 
(> 10i? ffi ) distances. The rigidity cutoff is calculated as the minimum 1Z for which the 
trajectory is allowed or, in other words, the backtracked trajectory escapes from the magnetic 
field of the Earth. 

The following features of the geomagnetic effects are relevant to the atmospheric neutrino 
fluxes: 

(i) For a fixed direction, the cutoff rigidity grows monotonically from zero at the magnetic 
pole to a maximum value at the magnetic equator. 

(ii) The cutoff rigidities for particles traveling toward the magnetic west are higher than 
those for particles traveling toward the magnetic east. 

(iii) Thus, the highest cutoff corresponds to westward going, horizontal particles reaching 
the surface of the Earth at the magnetic equator and it is approximately 60 GV. 

3. The neutrino yield 

As cosmic rays propagate in the atmosphere and interact with air nuclei, they create n- 
and if-mesons, which in turn decay and create atmospheric z/'s: 

A cr + A a { r — > 7r^~ , K^, K°, • • • ; 

n^K^fi* + upfr) (121) 
fJ^ -> + v e {v e ) + v^v^). 

We have only displayed the dominant channels. Since the decay distributions of the sec- 
ondary mesons and muons are extremely well known, the largest source of differences between 
the various calculations for the neutrino yield is the use of different models for the hadronic 
interactions (Gaisser et al, 1996). Since neutrinos produced by the decays of pions, kaons 
and muons have different energy spectra, the main features of the interaction model which 
affect the v^j v e composition and the energy and angular dependence of the neutrino fluxes 
are the K/n ratio and their momentum distribution. 

4. The neutrino fluxes 

In Figs. 23 and 24 we show the atmospheric neutrino fluxes from the calculation of the 
Bartol group (Gaisser, Stanev and Barr, 1988; Barr, Gaisser, and Stanev, 1989; Agraval et 
al., 1996; Lipari, Gaisser, and Stanev, 1998) for the location of SK and for maximum solar 
activity. In the first figure we show the flux as a function of the neutrino energy averaged 
over the arrival direction. In the second figure we show the fluxes as a function of the zenith 
angle for various neutrino energies. We emphasize the following points: 
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(i) Energy dependence: For E v ^ 1 GeV, the fluxes obey an approximate power law, 
dQ/dE oc E^ 1 where 7 ~ 3(3.5) for muon (electron) neutrinos and antineutrinos. For 
E v < 1 GeV, the dependence on energy is weaker as a consequence of the bending of the 
cosmic ray spectrum by geomagnetic effects and solar modulation. 

(ii) v/v ratio: For E < 1 GeV, when all pions and subsequent muons decay before 
reaching the Earth, the ratio v e jv e = tt~/tt + < 1. At these energies, the cosmic rays are 
mainly protons. The positively charged protons produce, on average, more 7r + 's than 7r~'s 
in their interactions. On the other hand, since n + produces a (z/ M ) in its decay, we 
expect v^jvp = 1. 

As energy increases, the secondary /^'s do not have time to decay before reaching the 
surface of the Earth and consequently v^j decreases. 

(iii) (z/ M + Vfj) j (v e + v e ) ratio: At E £ 1 GeV, the ratio is very close to 2 as expected from 
the chain decays in Eq. (121). At higher energies, the ratio decreases because, as mentioned 
above, the /z's do not have time to decay before reaching the Earth and less u' e s are produced. 

(iv) Up-down asymmetry: At E < a few GeV, the fluxes are up-down asymmetric due to 
geomagnetic effects. Geomagnetic effects are very small already at E = 2 GeV (see Fig. 24) 
and the flux becomes up-down symmetric above that energy. 

(v) Horizontal-vertical ratio: For E ^ a few GeV, the fluxes are maximal for neutrinos 
arriving horizontally and minimal for neutrinos arriving vertically (see Fig. 24). This is due 
to the difference in the atmosphere density which determines whether the pions decay before 
reinteracting with the air (thereby producing neutrinos). Pions arriving horizontally travel 
longer time in less dense atmosphere and are more likely to decay before reinteracting. 



B. Interaction Cross Sections 



In order to determine the expected event rates at the experiment, one needs to know 
the neutrino-nucleon interaction cross sections in the detector. The standard approach is 
to consider separately the contributions to this cross section from the exclusive channels 
of lower multiplicity: quasi-elastic scattering and single pion production, and include all 
additional channels as part of the deep inelastic (DIS) cross section: 



^CC* = &QE + 0"1tt + GDIS ■ 

In Fig. 25 we plot the cross sections for these processes. 

The cross section for the quasi-elastic interaction is given by (Smith, 1972): 
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where s — u = 4ME U + q 2 
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M is the proton mass, rri£ is the charged lepton mass, and 



q 2 is the momentum transfer. For up — > £ + n, a similar formula applies with the only change 
A 2 — > — A 2 . The functions Ai, A 2 , and A 3 can be written in terms of axial and vector form 
factors: 
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where we neglected second class currents and assumed CVC. With this assumption all form 
factors are real and can be written as follows: 
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Here /x p and //„ are the proton and neutron anomalous magnetic moments and M v = 
0.71 GeV 2 is the vector mass which is measured with high precision in electron scattering 
experiments. The largest uncertainties in this calculation are associated with the axial 
form factor. For example, the axial mass used by various collaborations varies in the range 
M\ = 0.71 - 1.06 GeV 2 . 

So far we have neglected nuclear effects. The most important of these effects is related 
to the Pauli principle and can be included by using a simple Fermi gas model (Smith, 1972). 
In this approximation, the cross section of a bound nucleon is equal to the cross section of 
a free nucleon multiplied by a factor of (1 — D/N). For neutrons, 
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(126) 



with z = ^{q 2 + m 2 ) 2 /(4M 2 ) - q 2 /{2k f ), w = {2N/A) 1 / 3 , and v = {2Z/A) 1 / 3 . Here A, Z, N 
are, respectively, the nucleon, proton and neutron numbers and kf is the Fermi momentum, 
kf = 0.225(0.26) GeV for oxygen (iron). For protons, the same formula applies with the 
exchange N <-> Z. The effect of this factor is to decrease the cross section. The decrease is 
larger for smaller neutrino energy. For energies above 1 GeV the nuclear effects lead to an 
8% decrease on the quasi-elastic cross section. As seen in Fig. 25, quasi-elastic interactions 
dominate for E v < 1 GeV and induce most of the observed contained events. 

An important role in the interpretation of the zenith angle distribution of the contained 
events is played by the relative angle between the direction of the incoming neutrino (carrying 
the information on the neutrino path length) and the direction of the produced charged 
lepton (which is measured). In Fig. 26 we plot this angle as a function of the measured 
charged lepton energy for quasi-elastic interactions. For energies below 1 GeV, the opening 
angle is rather large: for sub-GeV events the correlation between the measured i direction 
and the distance traveled by the neutrinos is weak. For E •> a few GeV, the two directions 
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are almost aligned (for example, it is in average 1° for upgoing muons) and the t direction 
gives a very good measurement of the neutrino path length. 

Single pion production, uN — > £~N'ir, is dominated by the A(1232) resonance (Fogli and 
Nardulli, 1979; Nakahata et. al, 1986). It is most relevant at E u ~ 1 GeV (see Fig. 25). 

Deep inelastic processes, vN — > £~X where X represents any hadronic system, dominate 
atmospheric neutrino interactions for E u > a few GeV. The parton model cross section is 
given by 



do-pis{v) _ G 2 F sx 
dxdy Arc 



F 1 -F 3 + (F 1 + F 3 )(l-y) 2 , (127) 



where y — 1 — Ei/E v and x = —q 2 /(2ME u y). For 9, a similar formula applies, with 
F 3 — > — F 3 . F 1 and F 3 are given in terms of the parton distributions. For isoscalar targets 
F\ = 2^2(qi + q~i) and F 3 = ^2(q~i — <ft). In order to avoid double counting of the single 

i i 

pion production, the deep inelastic contribution must be integrated in the region of hadronic 
masses W x > W c (W c = 1.4 GeV), which implies 2ME 1/ y{\ - x) > W 2 - M 2 . 



C. Two-Neutrino Oscillation Analysis 

The simplest and most direct interpretation of the atmospheric neutrino anomaly is that 
of muon neutrino oscillations (Learned, Pakvasa, and Weiler, 1988; Barger and Whisnant, 
1988; Hidaka, Honda and Midorikawa, 1988). The estimated value of the oscillation param- 
eters can be easily derived in the following way: 

• The angular distribution of contained events shows that, for E ~ 1 GeV, the deficit 
comes mainly from L ~ 10 2 — 10 4 km. The corresponding oscillation phase must be 
maximal, Am ^^j km) ~ 1, which requires Am 2 ~ 10" 4 - 10" 2 eV 2 . 

• Assuming that all upgoing u^s which would lead to multi-GeV events oscillate into a 
different flavor while none of the downgoing ones do, the up-down asymmetry is given 
by \Afj,\ = sin 2 20/(4 — sin 2 29). The present one sigma bound reads \A^\ > 0.27 [see 
Eq. (118)], which requires that the mixing angle is close to maximal, sin 2 26 l > 0.85. 

In order to go beyond these rough estimates, one must compare in a statistically meaningful 
way the experimental data with the detailed theoretical expectations. We now describe how 
to obtain the allowed region. We consider two neutrino cases, where oscillates into either 
v e or v T or a sterile neutrino v s . 



1. Predicted number of events 



For a given neutrino oscillation channel, the expected number of /i-like and e-like con- 
tained events, N a (a = //, e), can be computed as: 

= iV + N e» , N e = N ee + iV Me , (128) 

where 
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/<i 2 $ d<7 
— — j- ^ K a (h, cos 9 v ,E u )P a p— e(Ep)dE u dEpd(cos9 u )dh . (129) 

Here P a p is the conversion probability of v a — > vp for given values of E v , cos9 u and h, i.e., 
P a p = P(v a — > up; E u , cos 9 u ,h). In the SM, the only non-zero elements are the diagonal 
ones, i.e. P aa = 1 for all a. In Eq. (129), n t denotes the number of target particles, T is the 
experiment running time, E v is the neutrino energy, <& a is the flux of atmospheric v^s, Ep 
is the final charged lepton energy, s(Ep) is the detection efficiency for such a charged lepton, 
a is the neutrino-nucleon interaction cross section, and 9 U is the angle between the vertical 
direction and the incoming neutrinos (cos 0^=1 corresponds to the down-coming neutrinos). 
In Eq. (129), h is the slant distance from the production point to the sea level for a-type 
neutrinos with energy E v and zenith angle 9 V , and n a is the slant distance distribution, 
normalized to one. (We use in our calculations n a from Gaisser and Stanev, 1998). 

To obtain the expectation for the angular distribution of contained events one must 
integrate the corresponding bins for cos9p, where 9p is the angle of the detected lepton, 
taking into account the opening angle between the neutrino and the charged lepton directions 
as determined by the kinematics of the neutrino interaction. 

As discussed above, the neutrino fluxes, in particular in the sub-GeV range, depend on 
the solar activity. In order to take this fact into account, one uses in Eq. (129) a linear 
combination of atmospheric neutrino fluxes $™ ax and $™ in which correspond to the most 
(solar maximum) and least (solar minimum) active phases of the Sun, respectively, with 
different weights which depend on the specific running period. 

Experimental results on upgoing muons are presented in the form of measured muon 
fluxes. To obtain the effective muon fluxes for both stopping and through-going muons, 
one must convolute the survival probabilities for v^s with the corresponding muon fluxes 
produced by the neutrino interactions with the Earth. One must further take into account 
the muon energy loss during propagation both in the rock and in the detector, and also the 
effective detector area for both types of events, stopping and through-going. Schematically, 

§ ' w -i(zb)Cw v,£ '^' • (130) 



where 



POO POO POO POO 

N A dE^o dE v \ dX dhK„ (h,cos9,E„) 
Je u Je uQ Jo Jo 



dE^d cos 9 Je^ Je^ 

d^(E v ,9) dcrjE^E^) 

dKd cos 9 P ^ dE, F ^ E ^ M ' ( 131 ) 

Here Na is the Avogadro number, E^ is the energy of the muon produced in the neutrino 
interaction, E^ is the muon energy when entering the detector after traveling a distance 
X in the rock, and cos# labels both the neutrino and the muon directions which, at the 
relevant energies, are collinear to a very good approximation. We denote by .Frock (£^0) E^, X) 
the function which characterizes the energy spectrum of the muons arriving the detector. 
The Standard practice is to use an analytical approximation obtained by neglecting the 
fluctuations during muon propagation in the Earth. In this case the three quantities E^ , 
E^ and X are not independent: 
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f°° 1 

Jo F ^(E, ,E„X)dX= {dg ^ )/dxy (132) 

where (dS^E^/dX) is the average muon energy loss due to ionization, bremsstrahlung, 
e + e~ pair production and nuclear interactions in the rock (Lohmann, Kopp and Voss, 1985). 

For SK, the pathlength traveled by the muon inside the detector is given by the muon 
range function in water: 

L ^=C wmm iK - (133) 

In Eq. (130), A(L min , 9) = A S (E^, 9) + A T (E^, 9) is the projected detector area for internal 
path- lengths longer than L min (— 7 m in SK). Here A s and A T are the effective areas for 
stopping and through-going muon trajectories. These effective areas can be computed using 
a simple geometrical picture given, for instance, in Lipari and Lusignoli (1998). For a given 
angle, the threshold energy cut for SK muons is obtained by equating Eq. (133) to L min , 
i.e., L{Ef) = L min . 

In contrast to SK, MACRO present their results as muon fluxes for > 1 GeV, after 
correcting for detector acceptances. Therefore in this case we compute the expected fluxes 
as in Eqs. (130) and (131) but without the inclusion of the effective areas. 



2. Conversion probabilities 

We consider a two-flavor scenario, — > u x (X = e, r or s). The oscillation probabilities 
are obtained by solving the evolution equation of the v^—vx system in the matter background 
of the Earth (see section III.C): 

d / i/ M \ _ / H^x \(vp 



where 



H» = V»-—cos29, 
Am 2 

H x =Vx + -—cos 20, 
H fMX =-r E ^sin29. (135) 



The various neutrino potentials in matter are given by 

M v 2 ' 



" T M 2 nh 

V s = 0. (136) 
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Here p is the matter density in the Earth (Dziewonski and Anderson, 1981), M is the 
nucleon mass, and Y e (Y n ) is the electron (neutron) fraction. For anti- neutrinos, the signs of 
the potentials are reversed. 

For X = t, we have = V T and consequently these potentials can be removed from the 
evolution equation. The solution of Eq. (134) is then straightforward and the probability 
takes the well-known vacuum form [Eq. (35)], which is equal for neutrinos and anti-neutrinos. 
For X = e or s, the effect of the matter potentials requires a numerical solution of the 
evolution equations in order to obtain P a p which, furthermore, is different for neutrinos and 
anti-neutrinos. As first approximation, one can use a constant Earth matter density. Then 
(see section III) the solutions take the same form as the vacuum probability [Eq. (35)], but 
with the mixing angle and the oscillation length replaced by their effective values in matter 
[Eqs. (61) and (64) with A = 2E(V fl - V x )\. In Fig. 27 (from Lipari and Lusignoli, 1989) we 
show the survival probability of for the different oscillation channels for Am 2 = 5 x 10 -3 
eV 2 , sin 2 29 = 1 and various values of E v . As seen in the figure, matter effects damp 
the oscillation amplitude. For the chosen mass difference, they are important for neutrino 
energies of few 10 GeV and therefore are relevant mainly for upgoing muons. We return to 
this point below when describing the results of the analysis. 

3. Statistical analysis 

In order to define in a statistically meaningful way the regions of neutrino flavor param- 
eters that are allowed by a given set of atmospheric neutrino observables, one can construct, 
for example, a x 2 function: 

X 2 ee J2(N^ - N* ) • (er 2 ata + al)j} ■ (N$^ - N?), (137) 
i, J 

where / and J stand for any combination of experimental data sets and event-types con- 
sidered, i.e, I = (A, a) and J = (B,(3). The latin indices A, B stand for the different 
experiments or different data samples in a given experiment. The greek indices denote 
electron-type or muon-type events, i.e, a, j3 — e, p. In Eq. (137), Nf 1 stands for the pre- 
dicted number of events (or for the predicted value of the flux, in the case of upgoing muons) 
calculated as discussed above, whereas jVf ata is the corresponding experimental measure- 
ment. In Eq. (137), cr 2 ata an d ofh are ^ ne error matrices containing the experimental and 
theoretical errors, respectively. They can be written as 

a 2 u ee a a (A) Pa p(A, B) a p (B), (138) 

where p a p(A,B) is the correlation matrix containing all the correlations between the a- 
like events in the A-tjpe experiment and /3-like events in 5-type experiment, whereas 
<y a {A) [erg (5)] are the errors for the number of a-like events in the A [/3-like events in 
the B] experiment. The dimension of the error matrix varies depending on the combination 
of experiments included in the analysis. Details on the statistical analysis of the atmo- 
spheric neutrino data and results from the analyses with various data samples can be found 
in Fogli, Lisi and Montanino (1994,1995); Fogli and Lisi (1995b); Gonzalez- Garcia et al. 
(1998,1999); Gonzalez- Garcia, Fornengo and Valle (2000); Foot, Volkas and Yasuda (1998); 
Yasuda (1998); Akhmedov et al. (1999). 
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Using the definitions given above, one calculates \ 2 °f Eq. (137) as a function of the 
neutrino parameters. By minimizing x 2 with respect to sin 2 29 and Am 2 , one determines 
the best fit results, while the allowed regions are determined by the conditions: \ 2 = 
XL + 4.61 (6.1) [9.21] for a confidence level (CL) of 90 (95) [99] %. 

4. i>n -> i/ e 

At present z/^ — > z/ e is excluded with high CL as the explanation to the atmospheric 
neutrino anomaly for two different reasons: 

(i) SK high precision data show that the v e contained events are very well described by 
the SM prediction both in normalization and in their zenith angular dependence [see Fig. 21]. 
The distribution, however, shows an angle-dependent deficit, — > v e oscillations can 
explain the angular dependence of the flux only at the price of introducing angular 
dependence of the v e flux, in contrast to the data. Furthermore, even the best fit point for 
v n — > v e oscillations does not generate the observed up-down asymmetry in the Multi-GeV 
muon sample. This is illustrated in Fig. 28, where we show the predicted angular distribution 
of contained events at SK for the best fit points of the different oscillation channels. For 

— > v e oscillations, the asymmetry in the Multi-GeV muon distribution is much smaller 
than in the — > v T or — > v s channels. 

(ii) Explaining the atmospheric data with — > v e transition has direct implications for 
the v e — > transition. In particular, there should be a v e deficit in the CHOOZ reactor 
experiment. Thus the neutrino parameters not only give a poor fit to the atmospheric data 
but are actually excluded by the negative results from the CHOOZ reactor experiment (see 
Fig. 28). 

5. Ufj, -»■ v T and -> v s 

In Fig. 29 we show the values of the oscillation parameters (Am 2 , sin 2 9) which describe 
various sets of data for these two oscillation channels. The upper panels (labeled as FINKS) 
and the central panels refer to contained events. The upper panels take into account only 
the total rates, and the central ones only the angular distribution (from both Kamiokande 
and SK). The lower panels correspond to the angular distribution for upgoing muons from 
SK. The three shaded regions are allowed at the 90, 95 and 99% CL. In Fig. 30 we plot 
the allowed regions from the global analysis, including all the atmospheric neutrino data. 
Also shown are the expected sensitivity from disappearance at the long baseline (LBL) 
experiments K2K and MINOS discussed in VI. D. (These figures are an update of the results 
presented in Gonzalez- Garcia et al, 2000.) We emphasize the following points: 

(i) The allowed regions from the various data samples overlap. The oscillation hypothesis 
can then consistently explain the atmospheric neutrino data. 

(ii) The information from the total event rates alone is consistent with arbitrarily high 
Am 2 values. The reason is that no information on the minimum oscillation length can be 
inferred from the data [see section III. A]. 

(iii) The allowed regions for — * v T transition are symmetric with respect to maximal 
mixing. This must be the case because the corresponding probabilities take the vacuum 
expression and therefore depend on sin 2 29. 
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(iv) The allowed regions for — * v s transition are asymmetric due to Earth matter 
effects. These effects are more pronounced when the condition for maximal matter effect, 
Am 2 cos 29 ~ 2EV S , is fulfilled. Since in the chosen convention the potential difference 
A = 2E(V /X — V s ) < 0, the matter effects enhance neutrino oscillations for cos 29 < 
(sin 2 9 > 0.5). The opposite situation holds for antineutrinos, but neutrino fluxes are larger 
and dominate in the resulting effect. As a result the allowed regions are wider in the 
sin 2 9 > 0.5 side of the plot. 

(v) The best fit to the full data for — > u T corresponds to Am 2 = 2.6 x 10~ 3 eV 2 and 
sin 2 29 = 0.97. The best fit for ^ -> v s lies at Am 2 = 3 x 10~ 3 eV 2 and sin 2 # = 0.61. 

In order to discriminate between the — > v T and — > v s options, one can use the 
difference in the survival probabilities due to the presence of matter effects for oscillations 
into sterile neutrinos (Lipari and Lusignoli, 1998). As discussed above, the effect is important 
mainly for the higher energy neutrinos which lead to through-going muons events. In Fig. 31 
we plot the expected distributions for the best fit points for the two channels. As seen in the 
figure, the distribution is steeper for v T while for — > v a a flattening is observed for 

neutrinos coming close to the vertical due to the damping of the oscillation amplitude [see 
Fig. 27]. The data favours the steeper distributions and this translates into a better global fit 
for oscillations into v T . For the global analysis, an update of the results of Gonzalez-Garcia, 
et al. (2000) shows that for v T oscillations Xmin = 56/63 d.o.f while for v s oscillations 
Xmin = 72/63 d.o.f. As we will see in Sec. VII, this difference has an important role in A-u 
mixing schemes. 

SK has also used other methods to distinguish between the u T and v s hypotheses for ex- 
plaining atmospheric disappearance. One is to examine events likely to have been caused 
by NC interactions. While z/ r 's readily undergo such interactions, z/ s 's do not, resulting in 
a relative suppression of the NC signal (Vissani and Smirnov, 1998; Hall and Murayama, 
1998). Another method attempts to observe appearance of the newly created u T , even if only 
on a statistical basis, by selecting enriched samples. All methods strongly favour <-> v T 
oscillations over <-> v s (Fukuda et al, 2000; Toshito et al, 2001). 

VI. LABORATORY EXPERIMENTS 

Laboratory experiments to search for neutrino oscillations are performed with neutrino 
beams produced at either accelerators or nuclear reactors. In disappearance experiments, 
one looks for the attenuation of a neutrino beam primarily composed of a single flavor due 
to the mixing with other flavors. In appearance experiments, one searches for interactions 
by neutrinos of a flavor not present in the original neutrino beam. 

Most of the past and present laboratory experiments did not have an oscillation signal. 
In such discussed in Sec. III. A, the experiment sets a limit on the corresponding 

oscillation probability. Appearance experiments set limits {Pap) < Pl for given flavors 
a ^ j3. Disappearance experiments set limits (P aa ) > 1 — Pl for a given flavor a which, in 
the two neutrino case, can be translated into (Pap) < Pl for (3 ^ a. The results are usually 
interpreted in a two neutrino framework as exclusion regions in the (Am 2 , sin 2 29) plane. 
One can take the upper bound on the mixing angle in the asymptotic large Am 2 range and 
translate it back into the value of Pl. sin 2 29\ hn = 2P^ [see discussion in Sec. III. A]. The 
probability Pl is the relevant quantity when interpreting the results in the more-than-two 



51 



neutrino framework. 



A. Short Baseline Experiments at Accelerators 

Conventional neutrino beams from accelerators are mostly produced by n decays, with 
the pions produced by the scattering of the accelerated protons on a fixed target: 

p + target -> n ± +X (139) 

Thus the beam can contain both \x- and e-neutrinos and antineutrinos. The final composition 
and energy spectrum of the neutrino beam is determined by selecting the sign of the decaying 
7r and by stopping the produced \i in the beam line. 

Most oscillation experiments performed so far with neutrino beams from accelerators 
have characteristic distances of the order of hundreds of meters. We call them short baseline 
(SBL) experiments. With the exception of the LSND experiment, which we discuss below, 
all searches have been negative. In table V we show the limits on the various transition 
probabilities from the negative results of the most restricting SBL experiments. In Fig. 32 
(from Astier et al, 2001) we show the excluded regions corresponding to (the absence of) 
i/p — > v T and v e — > v T oscillations. Due to the short path length, these experiments are not 
sensitive to the low values of Am 2 invoked to explain either the solar or the atmospheric 
neutrino data but they are relevant for 4-z/ mixing schemes as we will see in Sec. VII. B. 



B. LSND and KARMEN 

The only positive signature of oscillations at a laboratory experiment comes from the 
Liquid Scintillator Neutrino Detector (LSND) running at Los Alamos Meson Physics Facility 
(Athanassopoulos et al, 1995, 1996, 1998). The primary neutrino flux comes from 7r + 's 
produced in a 30-cm-long water target when hit by protons from the LAMPF linac with 
800 MeV kinetic energy. The detector is a tank filled with 167 metric tons of dilute liquid 
scintillator, located about 30 m from the neutrino source. 

Most of the produced 7r + 's come to rest and decay through the sequence 7r + — > 
followed by ji + — > e + u e u^. The P M 's so produced have a maximum energy of 52.8 MeV. This 
is called the decay at rest (DAR) flux and is used to study — > V e oscillations. The energy 
dependence of the flux from decay at rest is very well known, and the absolute value is 
known to 7%. The open space around the target is short compared to the pion decay length. 
Thus only 3% of the 7r + 's decay in flight (DIF). The DIF flux is used to study — > u e 
oscillations. 

In 1995, the LSND experiment published data showing candidate events that are con- 
sistent with 9^ — > V e oscillations. Further supporting evidence was provided by the signal 
in the — > v e channel. We summarize here their main results in both the DAR and DIF 
channels. 
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For DAR related measurements, z/ e 's are detected in the quasi elastic process v e p — > e + n, 
in correlation with a monochromatic photon of 2.2 MeV arising from the neutron capture 
reaction np — > c/7. The main background is due to the v e component in the beam that is 
produced in the decay chain starting with 7r~'s. This background is suppressed by three 
factors. First, the n + production rate is about eight times the n~ production rate in the 
beam stop. Second, 95% of the 7r _ 's come to rest and are absorbed before decay in the beam 
stop. Third, 88% of the /x _ 's from 7r~'s DIF are captured from atomic orbit, a process which 
does not give a v e . Thus, the relative yield, compared to the positive channel, is estimated 
to be ~ (1/8) x 0.05 x 0.12 = 7.5 x 10~ 4 . 

In Athanassopoulos et al. (1995) LSND report a total of 22 events with e + energy 
between 36 and 60 MeV when 4.6 ±0.6 background events are expected. They fit the full e + 
event sample in the energy range 20 < E e < 60 MeV and the result yields 64.3^1^7 beam- 
related events. Subtracting the estimated neutrino background with a correlated gamma of 
12.5 ± 2.9 events results in an excess of 51.2tjJ;g ± 8.0 events. The interpretation of this 
anomaly in terms of — > u e oscillations requires P eM = (3.1 ± 1.2 ± 0.5) x 10 -3 . 

For DIF related measurements, the z/ e 's are observed via the detection of electrons pro- 
duced in the process v e C — > e~X with energy 60 < E e < 200 MeV. Using two independent 
analyses, 27.7 ±6.4 events are observed. The neutrino induced backgrounds are dominated 
by fi + — > e + u^u e and tt + — > e + u e decays in flight in the beam stop area and are estimated to 
be 9.6 ± 1.9 events. The excess above the expected background from conventional processes 
is then 18.1 ± 6.6 ± 4.0 events. The excess events are consistent with — > v e oscillations 
with P Me = (2.6 ± 1.0 ± 0.5) x lO- 3 . 

The LSND results have been recently updated to include the runs till 1998 (Aguilar 
et al, 2001) and the total fitted excess is of 87.9 ± 22.4 ± 6 events, corresponding to an 
oscillation probability of (2.64 ± 0.67 ± 0.45) x 10~ 3 . In the two-family formalism these 
results lead to the oscillation parameters shown in Fig. 33. The shaded regions are the 90 % 
and 99 % likelihood regions from LSND. The best fit point corresponds to Am 2 = 1.2 eV 2 
and sin 2 26 = 0.003. 

The region of parameter space which is favoured by the LSND observations has been 
partly tested by other experiments like the old BNL E776 experiment (Borodovsky et al, 
1992) and more recently by the KARMEN experiment (Gemmeke et al, 1990). The KAR- 
MEN experiment is performed at the neutron spallation facility ISIS of the Rutherford 
Appleton Laboratory. Neutrinos are produced by stopping the 800 MeV protons in a mas- 
sive beam stop target, thereby producing pions. The 7r~'s are absorbed by the target nuclei 
whereas 7r + 's decay at rest producing muon neutrinos via tc + — > The low momen- 

tum yU + 's are also stopped within the massive target and decay at rest, /i + — > e + v e p il . The 
tt + — /i + decay chain at rest gives a neutrino source with identical intensities for , v e 
and emitted isotropically. There is a minor fraction of ir~ decaying in flight (DIF) with 
the following fi~ DAR in the target station which leads to a very small contamination of 
v e jv e < 6.2 ■ 10~ 4 . The energy spectra of the z/'s are well defined due to the DAR of both 
the 7T + and 

The neutrinos are detected in a rectangular tank filled with 56 1 of a liquid scintillator. 
The signature for the detection of v e is a spatially correlated delayed coincidence of positrons 
from p{p e , e + )n with energies up to E e + = E Pe — Q = 51.0 MeV and 7 emission of either of 
the two neutron capture processes: p(n, 7)0? with one 7 of £'(7) = 2.2 MeV or Gd(n, 7) with 
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37-quanta on average and J2E(l) — 8MeV. 

The raw data presented in the EPS HEP 2001 conference (Wolf et al, 2001) correspond 
to ~9400C protons on target. Analyzing the data results in 11 sequential events which 
satisfy all cuts. This number is in good agreement with the total background expectation of 
12.3 ± 0.6. Applying a Bayesian renormalization procedure, an upper limit of iV(osc) < 6.3 
at 90%CL can be extracted. However using the spectral information and a maximum like- 
lihood analysis, KARMEN find a best fit value iV(osc) = within the physically allowed 
range of parameters, which can be translated into an upper limit of 3.8 and 3.1 oscillation 
events for Am 2 < 1 eV 2 and Am 2 > 20 eV 2 , respectively. These numbers are based on 
a complete frequentist approach as suggested by G. Feldman and R. Cousins. The corre- 
sponding exclusion curve in the two-neutrino parameter space is given in Fig. 33 together 
with the favoured region for the LSND experiment (from Wolf et al, 2001). At high Am 2 , 
KARMEN results exclude the region favored by LSND. At low Am 2 , KARMEN leaves 
some allowed space, but the reactor experiments at Bugey and CHOOZ add stringent limits 
for the larger mixing angles. This figure represents the final status of the LSND oscillation 
signal. 

The MiniBooNE experiment (Bazarko et al., 2000) searches for — > v e oscillations and 
is specially designed to make a conclusive statement about the LSND's neutrino oscillation 
evidence. They use a beam of energy 0.5 - 1.0 GeV initiated by a primary beam of 8 
GeV protons from the Fermilab Booster, which contains only a small intrinsic v e component 
(less than 0.3%). They search for an excess of electron neutrino events in a detector located 
approximately 500 m from the neutrino source. The MiniBooNE neutrino detector consists 
of 800 tons of pure mineral oil contained in a 12.2 m diameter spherical tank. A structure in 
the tank supports phototubes, which detect neutrino interactions in the oil by the Cerenkov 
and scintillation light that they produce. 

The L/E ratio is similar to that of LSND, giving MiniBooNE sensitivity to the same 
mode of oscillations. However, neutrino energies are more than an order of magnitude higher 
than at LSND, so that the search at MiniBooNE employs different experimental techniques. 
In Fig. 33 we show the 90% CL limits that MiniBooNE can achieve. Should a signal be 
found then the next step would be the BooNe experiment. 

C. Disappearance Experiments at Reactors 

Neutrino oscillations are also searched for using neutrino beams from nuclear reactors. 
Nuclear reactors produce V e beams with E v ~ MeV. Due to the low energy, e's are the 
only charged leptons which can be produced in the neutrino CC interaction. If the v e 
oscillated to another flavor, its CC interaction could not be observed. Therefore oscillation 
experiments performed at reactors are disappearance experiments. They have the advantage 
that smaller values of Am 2 can be accessed due to the lower neutrino beam energy. In table 
VI we show the limits on P ee from the negative results of the reactor experiments Gosgen 
(Zacek, G., 1986), Krasnoyarsk (Vidyakin et al, 1994), Bugey (Achkar et al, 1995), and 
CHOOZ (Apollonio et al, 1999). Gosgen, Krasnoyarks and Bugey have relatively short 
baselines while CHOOZ is the first long baseline (LBL) reactor experiment. 

In Fig. 34 we show the corresponding excluded regions in the parameter space for two 
neutrino oscillations. Bugey sets the strongest constraint on the allowed mixing in the Am 2 
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range that is interesting for the LSND signal. Due to its longer baseline, CHOOZ is sensitive 
to the lowest values of Am 2 , low enough to be in the range of interest for atmospheric 
neutrinos oscillations and the upper sector of the LMA solution for solar neutrinos. 

The CHOOZ experiment searches for disappearance of P e 's produced in a power station 
with two pressurized-water nuclear reactors with a total thermal power of 8.5 GW. At the 
detector, located at L ~ 1 km from the reactors, the v e reaction signature is the delayed 
coincidence between the prompt e + signal and the signal due to the neutron capture in the 
Gd-loaded scintillator. The ratio between the measured and expected fluxes averaged over 
the neutrino energy spectrum is given by 

R = 1.01 ± 2.8 %(stat) ±2.7 %(syst). (140) 

Thus no evidence was found for a deficit in the flux. The negative result is translated into 
exclusion regions in the (Am 2 , sin 2 26) plane shown in Fig. 34. The resulting 90% CL limits 
include Am 2 < 7 x 10~ 4 eV 2 for maximal mixing, and sin 2 29 < 0.10 for large Am 2 . The 
CHOOZ results are significant in excluding part of the region that corresponds to the LMA 
solution of the solar neutrino problem (see Sec. IV). Furthermore, the CHOOZ bound rules 
out with high significance the possibility that — > u e oscillations explain the atmospheric 
neutrino deficit. The constraint on the mixing angle is also relevant to the interpretation of 
the atmospheric neutrino anomaly in the framework of three-neutrino mixing. We return to 
this issue in Sec. V. 

Smaller values of Am 2 can be accessed at future reactor experiments using longer base- 
line. Pursuing this idea, the KamLAND experiment (Piepke et al, 2001), a 1000 ton liquid 
scintillation detector, is currently in operation in the Kamioka mine in Japan. This under- 
ground site is conveniently located at a distance of 150-210 km from several Japanese nuclear 
power stations. The measurement of the flux and energy spectrum of the z/ e 's emitted by 
these reactors will provide a test to the LMA solution of the solar neutrino anomaly. In 
Fig. 34 we plot the expected 90% sensitivity for the KamLAND experiment after 3 years 
of data taking (from Piepke et al., 2001). The experiment will, for the first time, provide 
a completely solar model independent test of this particle physics solution of the solar neu- 
trino problem. After a few years of data taking, it should be capable of either excluding the 
entire LMA region or, not only establishing u e <-> u ot \ lcv oscillations, but also measuring the 
oscillation parameters with unprecedented precision Data taking is expected to commence 
in 2002. 

D. Long Baseline Experiments at Accelerators 

Smaller values of Am 2 can also be accessed using accelerator beams at long baseline 
(LBL) experiments. In these experiments the intense neutrino beam from an accelerator is 
aimed at a detector located underground at a distance of several hundred kilometers. The 
main goal of these experiments is to test the presently allowed solution for the atmospheric 
neutrino problem by searching for either disappearance or v T appearance. At present there 
are three such projects approved: K2K (Nishikawa et al., 1997; Ahn et al, 2001; Nishikawa, 
et al, 2001), which runs with a baseline of about 235 km from KEK to SK, MINOS (Abies 
et al., 1995; Wojcicki, 2001) under construction with a baseline of 730 km from Fermilab 
to the Soudan mine where the detector will be placed, and OPERA (Shibuya et al, 1997; 
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Cocco et ai, 2000), under construction with a baseline of 730 km from CERN to Gran Sasso. 
With their expected sensitivities, these experiments can cover either some fraction or all of 
the parameter region suggested by the atmospheric neutrino anomaly discussed in Sec. V. 

In the K2K experiment, a wide-band, almost pure beam from 7r + decays is generated 
in the KEK 12-GeV/c Proton Synchrotron and a neutrino beam-line. The detector is in 
SK at a distance of 250 km. Various beam monitors along the beam line and two different 
types of front detectors (FDs) have also been constructed at the KEK site. The FDs are 
a lkt water Cerenkov detector, which is a miniature of the SK detector, and a so-called 
fine-grained detector which is composed of a scintillating fiber tracker trigger counters, lead 
glass counters and a muon range detector. The characteristics of the neutrino beam in the 
KEK site - direction, intensity, stability, energy spectrum and v e — composition - are 
examined using FDs and beam monitors. They are then extrapolated to the SK site and 
used to obtain the expected number of events and the energy spectrum. 

The K2K experiment had a successful start in early 1999, and data were recorded during 
several periods in 2000/1. The accumulated beam intensity during the 2000 runs was 22.9 x 
10 18 protons on target (p.o.t.) (Ahn et ai, 2001), increased to ~ 38 x 10 18 p.o.t. with the 
2001 run till the summer (Nishikawa, et ai, 2001) which was about 40% of the goal of the 
experiment, 10 20 p.o.t.. The no-oscillation prediction for this sample, based on the data 
from the FDs, is 63.9l|g events while a total of 44 events have been observed in 22.5kt 
of the fiducial volume. The statistical probability that the observation would be equal to 
or smaller than 44 is ~ 8%. In the presence of oscillations the expected number of events 
would be 41.5 for Am 2 = 3 x 10~ 3 eV 2 and maximal mixing. Although the central value 
of the number of observed events is consistent with the data from atmospheric neutrinos, 
the discrepancy with the no-oscillation prediction is still within the statistical error. K2K 
has also studied the energy distribution of these events compared to the expectation based 
on the pion monitor data and a Monte-Carlo simulation and they find that the data are 
consistent within statistics with an oscillation signal. 

MINOS is designed to detect neutrinos delivered by the Main Injector accelerator at 
Fermilab (NuMI) with average energies of ~ 5-15 GeV depending on the beam configuration. 
Two detectors, functionally identical, will be placed in the NuMI neutrino beam: one at 
Fermilab and the second one in Soudan iron mine, 732 km away. The MINOS detectors are 
iron/scintillator sampling calorimeters with a toroidal magnetic field in the iron. Observed 
interactions of z/ M can be divided into two classes: CC-like, with an identified fi track, and 
NC-like, muonless. The ratio of the observed numbers of the CC and NC-like events in the 
two detectors provides a sensitive test for oscillations. With its expected sensitivity MINOS 
will be able to precisely measure (roughly at the level of 10%) the oscillation parameters in 
the z/ M — > v T channel. The primary measurement for this is the comparison of the rate and 
spectrum of the CC events in the Far Detector with those in the Near Detector. Comparing 
the NC/CC ratios in the two detectors, the experiment can also be sensitive to the presence 
of z/ M — > v s oscillations. MINOS is scheduled to start data taking at the end of 2004. 

Both K2K and MINOS have also the capability for detecting the appearance of z/ e 's due 
to — > v e oscillations. This signal however suffers from large backgrounds due to both the 
z/ e 's in the beam and the NC events with a topology similar to the v e interaction. These 
backgrounds can be partially suppressed using the information from the near detectors. In 
particular, MINOS may be able to improve the CHOOZ bounds. 
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OPERA is designed to search for — > v T oscillations in the Gran Sasso Laboratory. It 
will study the interaction of 20 GeV neutrinos produced at CERN. The goal is to observe 
the appearance of z/ r 's in a pure beam. The detector is based on a massive lead/nuclear 
emulsion target. Nuclear emulsions are exploited for the direct observation of the decay of 
the r in a very low background environment. 

E. Direct Determination of Neutrino Masses 

Oscillation experiments have provided us with important information on the differences 
between the neutrino masses-squared, Am?-, and on the leptonic mixing angles, Uij. But 
they are insensitive to the absolute mass scale for the neutrinos, m,. 

Of course, the results of an oscillation experiment do provide a lower bound on the 
heavier mass in Am^, \rrii\ > ^ Am? • for Am?- > 0. But there is no upper bound on this 
mass. In particular, the corresponding neutrinos could be approximately degenerate at a 
mass scale that is much higher than Amfj . Moreover, there is neither upper nor lower 
bound on the lighter mass rrij. 

Information on the neutrino masses, rather than mass differences, can be extracted from 
kinematic studies of reactions in which a neutrino or an anti-neutrino is involved. In the 
absence of mixing the present limits are (Groom et al, 2000) 

m„ T < 18.2 MeV (95% CL) from r~ -> nn + v T (141) 
m„ M < 190 keV (90% CL) from tt - -> yT (142) 
m Ue < 2.2 eV (95% CL) from 3 H -> 3 He + e" + V e , (143) 

where for the bound on m Ue we take the latest limit from the Mainz experiment (Bonn et 
al, 2001). A similar bound is obtained by Troitsk experiment (Lobashev et al, 2001). A 
new experimental project, KATRIN, is under consideration with an estimated sensitivity 
limit: m Ve ~ 0.3 eV. 

In the presence of mixing these limits have to be modified and in general they involve 
more than a single flavor parameter. The limit that is most relevant to our purposes is the 
most sensitive one from tritium beta decay. In presence of mixing, the electron neutrino 
is a combination of mass eigenstates and the tritium beta decay spectrum is modified as 
(Shrock, 1980): 

^ = R{E) J2 \U et \ 2 [(Eo - Ef - m t 2 ]^Q(E -E- m,), (144) 

i 

where E is the energy of electron, E is the total decay energy and R(E) is m^-independent. 
The step function, 0(-Eo — E — m^), reflects the fact that a given neutrino can only be 
produced if the available energy is larger than its mass. According to Eq. (144), there 
are two important effects, sensitive to the neutrino masses and mixings, on the electron 
energy spectrum: (i) Kinks at the electron energies E® = E ~ Eo-rrii with sizes that are 
determined by |£4«| 2 ; ( u ) A shift of the end point to E ep = Eq — mi, where mi is the lightest 
neutrino mass. The situation simplifies considerably if we are interested in constraining the 
possibility of quasi-degenerate neutrinos with mass ~ m v . In this case the distortion of 
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the spectrum can be described by a single parameter (Vissani, 2001a; Farzan, Peres and 

m-\U I 2 

Smirnov, 2001), ma = ^ ' % ~ m v . So the limit in Eq. (143) applies to the unique 
neutrino mass scale. 

Direct information on neutrino masses can also be obtained from neutrinoless double 
beta decay (2(30v) searches: 

(A,Z) -> (A,Z + 2) + e- +e~. (145) 
The rate of this process is proportional to the effective Majorana mass of u e , 



(146) 



which, in addition to five parameters that affect the tritium beta decay spectrum, depends 
also on the three leptonic CP violating phases. Notice that in order to induce the 2(30u 
decay, v e must be a Majorana particle. 

The present strongest bound from 2/30z/-decay is obtained by Heidelberg-Moscow group 
(Klapdor-Kleingrothaus et al, 2001): 

m ee < 0.34 (0.26) eV, 90 % (68%) C.L.. (147) 

Taking into account systematic errors related to nuclear matrix elements, the bound may 
be weaker by a factor of about 3. A sensitivity of m ee ~ 0.1 eV is expected to be reached 
by the currently running NEM03 experiment (Marquet et al, 2000), while a series of new 
experiments (CUORE, EXO, GENIUS) is planned with sensitivity of up to m ee ~ 0.01 eV. 

The knowledge of m ee will provide information on the mass and mixing parameters that 
is independent of the Am?'s. However, to infer the values of neutrino masses, additional 
assumptions are required. In particular, the mixing elements are complex and may lead 
to strong cancellation, m ee <C m\. Yet, the combination of results from 2/50z/ decays and 
Tritium beta decay can test and, in some cases, determine the mass parameters of given 
schemes of neutrino masses (Vissani, 1999; Farzan, Peres and Smirnov, 2001; Bilenky, Pascoli 
and Petcov, 2001a, 2001b; Klapdor-Kleingrothaus, Pas and Smirnov, 2001; Pascoli and 
Petcov 2002). 



VII. THREE- AND FOUR-NEUTRINO MIXING 

In the previous sections we discussed the three pieces of evidence for neutrino masses and 
mixing (solar neutrinos, atmospheric neutrinos and the LSND results) as usually formulated 
in the framework of two-neutrino oscillations. The results are summarized in Fig. 35 where 
we show the ranges of masses and mixing implied by these signals at 90 and 99% CL for 2 
d.o.f., as well as relevant constraints from negative searches in laboratory experiments. 

The three pieces of evidence correspond to three values of mass-squared differences of 
different orders of magnitude. Consequently, there is no consistent explanation to all three 
signals based on oscillations among the three known neutrinos. The argument for this 
statement is very simple. With three neutrinos, there are only two independent mass- 
squared differences, since the following relation must hold: 
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Am 2 21 + Am 2 32 + Am\ 3 = 0. (148) 

This relation cannot be satisfied by three Am 2 - that are of different orders of magnitude. 
One may wonder if this naive extrapolation from the two-neutrino oscillation picture holds 
once the full mixing structure of the three-neutrino oscillations in taken into account or, 
on the contrary, some special configuration of the three-neutrino parameters could fit the 
three pieces of evidence. The combined fit to the data performed by Fogli, Lisi, Montanino 
and Scioscia (1999a) shows that this is not the case (see also Gonzalez-Garcia and Maltoni, 
2002). 

Whereas in the case of the solar and atmospheric neutrino indications, several experi- 
ments agree on the existence of the effect, the third indication is presently found only by the 
LSND experiment. Therefore, in many studies the LSND result is left out and the analysis 
of the solar and atmospheric data is performed in the framework of mixing between the 
three known neutrinos. In Sec. VILA we discuss the derived masses and mixing in these 
scenarios. On the other hand, if all three indications in favour of neutrino oscillations are 
confirmed, one needs a minimum of four neutrinos. The present phenomenological status 
of the possibility of mixing between four neutrinos is discussed in Sec. VII. B. Alternatively 
a possible explanation of the four pieces of evidence with only three neutrinos has been 
proposed assuming that CPT is violated in the neutrino sector (see for instance, Murayama 
and Yanagida, 2001; Barenboim, Borissov, Lykken and Smirnov, 2001). 



A. Three-Neutrino Mixing 

The combined description of both solar and atmospheric anomalies requires that all three 
known neutrinos take part in the oscillations. The mixing parameters are encoded in the 
3x3 lepton mixing matrix (Maki, Nakagawa and Sakata, 1962; Kobayashi and Maskawa, 
1973). The two Majorana phases do not affect neutrino oscillations (Bilenky, Hosek and 
Petcov, 1980; Langacker et ai, 1987) . The Dirac phase (that is the analog of the KM phase 
of the quark sector) does affect neutrino oscillations in general, but for the purposes of this 
section we can set it to zero. 

In this case the mixing matrix can be conveniently parametrized in the standard form 
(Groom et al, 2000): 



U = 



C13C12 S12C13 S13 



-S12C23 - S23S13C12 C23C12 - S23S13S12 S23C13 I , (149) 

\ S23S12 ~~ s 13 c 23 c 12 — S23 c 12 — Sl3 s 12C23 C23 c 13 



where c^- = cos % and Sij = sin % . 

As we have seen in the previous sections, in most of the parameter space of solutions for 
solar and atmospheric oscillations, the required mass differences satisfy 

Ami « A ™atm- (150) 

In this approximation the angles #y can be taken without loss of generality to lie in the first 
quadrant, G [0,7r/2]. There are two possible mass orderings which we chose as 
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Am 2 21 = Ami < A ™32 - ^ = Am* tm > 0; (151) 
A< = Am| < -Am 2 31 ~ -Am^ = |Am" tm | > 0. (152) 

We refer to the first option, Eq. (151), as the direct scheme, and to the second one, Eq. (152), 
as the inverted scheme. The direct scheme is naturally related to hierarchical masses, mi <C 
rri2 <C W3, for which m 2 — \J Am 21 and m 3 ~ \j Am^, or to quasi-degenerate masses, 
m! ~ m 2 — m 3 ^> Am^Am^. On the other hand, the inverted scheme implies that 
m 3 < mi ~ m 2 . 

One may wonder how good an approximation it is to set the CP violating phases to 
zero. It turns out to be an excellent approximation for the analysis of solar, atmospheric 
and laboratory data if Eq. (150) holds. In this discussed below, no simultaneous 

effect of the two mass differences is observable in any //-appearance transition. 



1. Probabilities 

The determination of the oscillation probabilities for both solar and atmospheric neu- 
trinos requires that one solves the evolution equation of the neutrino system in the matter 
background of the Sun or the Earth. In a three-flavor framework, this equation reads: 

i^ = Hu, H = U -H^-U^ + V , (153) 
at 

where U is the lepton mixing matrix, z7= (z/ e , v T ) , H{j is the vacuum hamiltonian, 

H$ = ^diag (-Am 2 21 , 0, Am 2 32 ) , (154) 

and V is the effective potential that describes charged-current forward interactions in matter: 

V = diag (±V2G F N e , 0, o) = diag (V e , 0, 0) . (155) 

In Eq. (155), the sign + (— ) refers to neutrinos (antineutrinos), and N e is electron number 
density in the Sun or the Earth. 

In what follows we focus on the direct scheme of Eq. (151), for which the five relevant 
parameters are related to experiments in the following way: 

Ami = Am 2i, Am^n = Am\ 2) (156) 

= #12, #atm = #23) ^reactor = #13 • (157) 

For transitions in vacuum, the results apply also to the inverted scheme of Eq. (152). In 
the presence of matter effects, the direct and inverted schemes are no longer equivalent, 
although the difference is hardly recognizable in the current solar and atmospheric neu- 
trino phenomenology as long as Eq. (150) holds (Fogli, Lisi, Montanino and Scioscia, 1997; 
Gonzalez-Garcia and Maltoni, 2002) . Under this approximation, the results obtained for 
the direct scheme can be applied to the inverted scheme by replacing Am\ 2 — > —Am 32 . 

In general the transition probabilities present an oscillatory behaviour with two oscil- 
lation lengths. However, the hierarchy in the splittings, Eq. (150), leads to important 
simplifications. 
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Let us first consider the analysis of solar neutrinos. A first simplification occurs because 
Lg| c = 47iE/Am 32 is much shorter than the distance between the Sun and the Earth. 
Consequently, the oscillations related to Lg| c are averaged in the evolution from the Sun to 
the Earth. A second simplification occurs because, for the evolution in both the Sun and the 
Earth, Am^ 3> 2y/2GpN e E l/ sin 2 29± 3 . Consequently, matter effects on the evolution of v 3 
can be neglected. The result of these two approximations is that the three-flavor evolution 
equations decouple into an effective two-flavor problem for the basis (Kuo and Pantaleone, 
1986; Shi and Schramm, 1992) 

V e , = Cx2 V\ + Sl2^2 , iy = ~S\2V\ + c 12 v 2 , (158) 

with the substitution of N e by the effective density 

N e =>■ N e cos 2 #13 . (159) 

Thus the survival probability takes the following form: 

P£ = sin 4 9 13 + cos 4 13 P e 2 £, , (160) 

where P 2 >e',o is the two-flavor survival probability in the (Am 2 ^,^) parameter space but 
with the modified matter density of Eq. (159). We conclude that the analysis of the solar 
data constrains three of the five independent oscillation parameters: Am^, 9\ 2 and #13. 

Eq. (160) reveals what is the dominant effect of a non- vanishing #13 in the solar neutrino 
survival probability: the energy dependent part of the probability, P 2u e i Q , gets damped by 
the factor cos 4 6*13, while an energy independent term, sin 4 #13, is added. Thus increasing #13 
makes the solar neutrino survival probability more and more energy independent. 

Let us now consider the analysis of atmospheric neutrinos. Here Lfi = AnE/ Am^ is 
much larger than the relevant distance scales. Consequently, the corresponding oscillating 
phase is negligible. In this approximation one can rotate away the corresponding angle 9\ 2 . 
Thus the resulting survival probabilities do not depend on Am 2 x and 9 12 . For instance for 
constant Earth matter density, the various P a p can be written as follows: 

-fee — 1 — ^ s 13,m C 13,m ^31; 

E> _ 1 A 2 2 4 o /I 2 2 2 o /I 2 2 2 o 

r m — 1 - 4Sl3,m C 13,m S 23 ^31 - 4s 13,m S 23 C 23 D 21 - ^C\ 3 , m S 23 C 23 b 32 , 

P TT = 1- 4s 2 3 m c?3 im c4 - 4:S 2 13tm s 2 23 c 2 23 S 21 - ^c\ 3 m s\ 3 (? 23 S 32 , (161) 

P - ZLe 2 r 2 e 2 <7 
r e/i — ^*13,m L 13,m*23 °31) 

Per = 4Si 3 ,m C 13,m C 23 ^31, 

P MT = -4s 2 3 im c 2 3 im s^c^3 531 + 4:s 2 13 m s 2 23 c 2 23 S 21 + 4:C 2 13tm s 2 23 c 2 23 S 32 . 
Here #i 3 , m is the effective mixing angle in matter: 

• oa sin 2#i3 

sin2fli 3im = (162) 

v /(cos2#i3 - 2E u V e /Am 2 32 y + (sin2^ 3 ) 2 

and Sij are the oscillating factors in matter: 
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S i3 = sin 2 (j^l) . (163) 



V 



In Eq. (163), A/i^- are the effective mass-squared differences in matter: 

A/4, - =^f- I —— -1) I I . 



l 32 



l 3l 



Am 2 2 


/ sin2# 13 


2 


\ v sin26'i3 im 


_ Am 2 2 


/ sin 2^3 


2 


\sin26>i 3im 


= Am 2 2 


sin 26*i 3 


sin26»i 3>m ' 



+ l J + ^K, 

(164) 



and L is the pathlength of the neutrino within the Earth, which depends on its direction. We 
conclude that the analysis of the atmospheric data constrains three of the five independent 
oscillation parameters: Aro^, 6> 23 and 9 13 . 

So we find that in the approximation of Eq. (150) the mixing angle 6> 13 is the only 
parameter common to both solar and atmospheric neutrino oscillations and which may 
potentially allow for some mutual influence. The main effect of the three-neutrino mixing 
is that now atmospheric neutrinos can oscillate simultaneously in both the — > v T and 

— > v e (and, similarly, v e — > v T and u e — > v^) channels. The oscillation amplitudes for 
channels involving v e are controlled by the size of sin 2 9 13 = \U e3 \ 2 . We learn that in the 
approximation of Eq. (150) solar and atmospheric neutrino oscillations decouple in the limit 
$13 = 0. This angle is constrained by the CHOOZ reactor experiment. To analyze the 
CHOOZ constraints we need to evaluate the survival probability for v e of average energy 
E ~ few MeV at a distance of L ~ 1 Km. For these values of energy and distance, one can 
safely neglect Earth matter effects. The survival probability takes the analytical form: 

pCHOOZ = x _ cog 4 sin 2 20i2 gin 2 (165) 



- sin 2 26 



13 



cos 2 9 12 sin 2 f^Usin 2 W^ Am - r 



AE V \ 4£„ 



~l-sin 2 2^3sin 2 ^ Am ^ L> 



V 4^ 



where the second equality holds under the approximation Am^ E u /L which can only be 
safely made for Am 2 x < 3 x 10 _4 eV 2 . Thus in general the analysis of the CHOOZ reactor date 
involves four oscillation parameters: Am 2 2 , 6> 13 , Am^, and 9 12 (Gonzalez- Garcia, Maltoni, 
2002; Bilenky, Nicolo and Petcov, 2001) 



2. Allowed masses and mixing 

There are several three-neutrino oscillation analyses in the literature which include either 
solar (Fogli, Lisi and Montanino, 1996; Fogli, Lisi, Montanino and Palazzo, 2000a, 2000b; 
Gago, Nunokawa and Zukanovich, 2001) or atmospheric (Fogli, Lisi, Montanino and Scioscia, 
1997; Fogli, Lisi, Marrone and Scioscia, 1999b; Fogli, Lisi, Marrone, 2001b; de Rujula, Gavela 
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and Hernandez, 2001; Teshima and Sakai 1999) neutrino data. Combined studies have also 
been performed (Barger, Whistnant and Phillips, 1980; Fogli, Lisi and Montanino, 1994, 
1995; Barbieri et ai, 1998; Barger and Whisnant, 1999). We follow and update here the 
results from the analysis of Gonzalez- Garcia, Maltoni, Pena-Garay and Valle (2001). 

As discussed above, in the approximation of Eq. (150), solar and atmospheric neutrino 
oscillations decouple in the limit |£/e 3 | = sin #13 = 0. In this limit the allowed values of the 
parameters can be obtained directly from the results of the analyses in terms of two-neutrino 
oscillations presented in Sec. IV and Sec. V. Deviations from the two-neutrino scenario are 
then determined by the size of 9 13 . Thus the first question to answer is how the presence of 
this additional angle affects the analysis of the solar and atmospheric neutrino data. 

In Fig. 36 we show the allowed regions for the oscillation parameters Am^ and tan 2 812 
from the global analysis of the solar neutrino data in the framework of three-neutrino oscil- 
lations for different values of sin 2 6*13 (updated from Gonzalez-Garcia, Maltoni, Peha-Garay 
and Valle, 2001). The allowed regions for a given CL are defined as the set of points satisfy- 
ing the condition x 2 (Am 2 2 , tan 2 9 12 , tan 2 13 ) — Xmin ^ Ax 2 (CL, 3 d.o.f) where, for instance, 
Ax 2 (CL, 3 d.o.f)=6.25, 7.83, and 11.36 for CL=90, 95, and 99%, respectively. The global 
minimum used in the construction of the regions lies in the LMA region and corresponds to 
tan 2 9 13 = 0, that is, to the decoupled scenario. 

As seen in the figure, the modifications to the decoupled case are significant only if 9 13 is 
large. As sin 2 6*13 increases, all the allowed regions disappear, leading to an upper bound on 
sin 2 #i3 that is independent of the values taken by the other parameters in the three-neutrino 
mixing matrix. For instance, no region of parameter space is allowed (at 99% C.L. for 3 
d.o.f) for sin 2 ^13 = |t/ e3 | 2 > 0.80. This fact is also illustrated in Fig. 37 where we plot the 
shift in x 2 as a function of sin 2 # 13 when the mass and mixing parameters Amf 2 and tan 2 12 
are chosen to minimize \ 2 ■ 

In Fig. 38 we show the allowed regions for the oscillation parameters Am 32 and tan 2 #23 
from the analysis of the atmospheric neutrino data in the framework of three-neutrino oscil- 
lations for different values of sin 2 #13 (updated from Gonzalez-Garcia, Maltoni, Pena-Garay 
and Valle, 2001). In the upper-left panel tan 2 13 = which corresponds to pure — > v T 
oscillations. Note the exact symmetry of the contour regions under the transformation 
$23 7r /4 — #23- This symmetry follows from the fact that in the pure — > v T channel 
matter effects cancel out and the oscillation probability depends on #23 only through the 
double- valued function sin 2 26*23 (see Sec. V.C). For 9\% 7^ this symmetry breaks due to 
the three-neutrino mixing structure even if matter effects are neglected. The analysis of 
the full atmospheric neutrino data in the framework of three-neutrino oscillations clearly 
favours the — > v T oscillation hypothesis. As a matter of fact, the best fit corresponds 
to a very small value, #i 3 = 6°, but it is statistically indistinguishable from the decoupled 
scenario, #i 3 = 0°. No region of parameter space is allowed (at 99% C.L. for 3 d.o.f) for 
sin 9i 3 = \U e3 \ 2 > 0.40. The physics reason for this limit is clear from the discussion of the 
v n v e oscillation channel in Sec. V.C: large values of 13 imply a too large contribution 
of the — > v e channel and would spoil the otherwise successful description of the angular 
distribution of contained events. This situation is illustrated also in Fig. 37 where we plot 
the shift of x 2 f° r the analysis of atmospheric data in the framework of oscillations between 
three neutrinos as a function of sin 2 #13 when the mass and mixing parameters Am 32 and 
tan 2 6*25 are chosen to minimize x 2 . 
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For any value of the mixing parameters, the mass-squared difference relevant for the 
atmospheric analysis is restricted to lie in the interval 1.25 x 1CT 3 < Am^/eV 2 < 8 x 1CT 3 
at 99 % CL. Thus it is within the range of sensitivity of the CHOOZ experiment. Indeed, 
as illustrated in Fig. 37, the limit on sin 2 9 13 becomes stronger when the CHOOZ data are 
combined with the atmospheric and solar neutrino results. 

One can finally perform a global analysis in the five dimensional parameter space com- 
bining the full set of solar, atmospheric and reactor data. Such analysis leads to the following 
allowed 3a ranges for individual parameters (that is, when the other four parameters have 
been chosen to minimize the global x 2 ): 



2.4 x 1(T 5 < Am 2 21 /eV 2 < 2.4 x 10~ 4 LMA 
0.27 < tan 2 9 12 < 0.79 LMA 
1.4 x 10~ 3 < Am 2 2 /eV 2 < 6.0 x 10~ 3 
0.4 < tan 2 #23 < 3.0, 

sin 2 < 0.06 (166) 

These results can be translated into our present knowledge of the moduli of the mixing 
matrix U : 



\U\ 



( 0.73 - 0.89 0.45 - 0.66 < 0.24 

0.23- 0.66 0.24-0.75 0.52 - 0.87 | . (167) 
V 0.06 - 0.57 0.40 -0.82 0.48 -0.85, 



In conclusion, we learn that at present the large mixing-type solutions provide the best fit. 
As concerns \U e3 \, both solar and atmospheric data favour small values and this trend is 
strengthened by the reactor data. 



B. Four-Neutrino Mixing 

If all three indications of neutrino oscillations are confirmed, then a minimum of four 
neutrinos will be required to accommodate the data. As discussed in Sec. II, the measure- 
ment of the decay width of the Z° boson into neutrinos makes the existence of three, and 
only three, light (that is, m v <m z /2) active neutrinos an experimental fact (Eq. 6). There- 
fore, a fourth neutrino must not couple to the standard electroweak current, that is, it must 
be sterile. 

One of the most important issues in the context of four-neutrino scenarios is the four- 
neutrino mass spectrum. There are six possible four-neutrino schemes, shown in Fig. 39, 
that can accommodate the results from solar and atmospheric neutrino experiments as well 
as the LSND result. They can be divided in two classes: (3+1) and (2+2) (Barger et al, 
2001). In the (3+1) schemes, there is a group of three close-by neutrino masses that is 
separated from the fourth one by a gap of the order of 1 eV 2 , which is responsible for the 
SBL oscillations observed in the LSND experiment. In (2+2) schemes, there are two pairs 
of close masses separated by the LSND gap. The main difference between these two classes 
is the following: if a (2+2)-spectrum is realized in nature, the transition into the sterile 
neutrino is a solution of either the solar or the atmospheric neutrino problem, or the sterile 
neutrino takes part in both, whereas with a (3+l)-spectrum the sterile neutrino could be 
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only slightly mixed with the active ones and mainly provide a description of the LSND 
result. 

As concerns the mixing parameters, we emphasize that the mixing matrix that describes 
CC interactions in these schemes is a 3 x 4 matrix. The reason is that there are three 
charged lepton mass eigenstates (e,/z,r) and four neutrino mass eigenstates (^1,^2,^3,^4)- 
As discussed in Sec. II. D, if we choose an interaction basis where the charged leptons are the 
mass eigenstates, then the CC mixing matrix U is a sub-matrix of the 4x4 unitary matrix 
V v that rotates the neutrinos from the interaction basis to the mass basis, where the line 
corresponding to v s is removed. 

1. Status of (3+1) schemes 

It has been argued in the literature that the (3+l)-spectra are strongly disfavoured by 
the data from SBL laboratory experiments (Okada and Yasuda, 1997; Bilenky, Giunti, and 
Grimus, 1998; Barger, Pakvasa, Weiler and Whisnant, 1998a; Bilenky, Giunti, Grimus and 
Schwetz, 1999). At what level they can actually be ruled out has been a matter of debate and 
the answer has also changed with the changes in the analysis of the LSND and KARMEN 
collaborations (Barger et al, 2000; Giunti and Laveder, 2001; Peres and Smirnov, 2001). 

The arguments which disfavour the (3+l)-mass spectra are based on exclusion curves 
from SBL experiments, and on some simplified treatment of the solar and atmospheric 
neutrino data. We summarize here these arguments and the present phenomenological 
status as given in Grimus and Schwetz (2001), Maltoni, Schwetz, and Valle (2001a, 2001b) 
and Maltoni, Tortola, Schwetz, and Valle (2002b). 

The probability that is relevant for LSND (and also for KARMEN and NOMAD) is 
given by 

A ™2 T 

P,e = Pne = MU e4 U^\ 2 sin 2 , (168) 

where L is the distance between source and detector. In this expression solar and atmospheric 
splittings have been neglected as they are too small to give any observable effect at the L/E 
that is relevant for LSND. In this approximation, for schemes (3+l)A, Am 2 = Am| x = 
Amf 2 = Am 2 3 . 

The LSND experiment gives then an allowed region in the Am 2 — \U e JJ^\ 2 plane which 
can be directly obtained from the two-neutrino oscillation region shown in Fig. 33 with the 
identification sin 2 29 — > 4\U e 4U^\ 2 . In the same way, the KARMEN experiment gives an 
excluded region in the same plane which can be directly obtained from Fig. 33. 

Further constraints on \U e4i U^\ 2 can be obtained by combining the bounds on |?7 e4 | and 
I £7^4 1 from reactor and accelerator experiments in combination with the information from 
solar and atmospheric neutrinos. The strongest constraints in the relevant Am 2 region are 
given by the Bugey (Achkar et al, 1995) and CDHS (Dydak et al, 1984) experiments. Using 
their limits on the survival probabilities one finds 

4|t/ e4 | 2 (l - \U e4 \ 2 ) < D^fAm 2 ), 

D e Bugcy < 0.001 - 0.1 (for 0.1 < Am 2 /eV 2 < 10), (169) 
4|^ 4 | 2 (l-|^ 4 | 2 )<^ CDHS (Am 2 ), 

^cdhs < g 05 _ 0.1 (for Am 2 > 0.5 eV 2 ), (170) 
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while for lower Am 2 the CDHS bound weakens considerably. 

In principle, the inequalities in Eq. (169) and Eq. (170) can be satisfied with either small 
mixing parameters \U a ^\ 2 < D a /A, or close to maximal mixing, \U ai \ 2 ^ 1 — D a /4. This 
is the generalization to these schemes of the symmetry 9 <-> | — 9 of two-neutrino vacuum 
oscillations (see Sec. III. A). Solar and atmospheric data are invoked in order to resolve this 
ambiguity. Since in any of the allowed regions for solar oscillations P® ^ 0.5, only the small 
values of \U e 4\ are possible. For atmospheric neutrinos one can use the fact that oscillations 
with the large Am 2 would wash-out the up-down asymmetry and this wash-out grows with 
the projection of the over states separated by the large Am 2 , which is controlled by the 
mixing parameter |L^ 4 |. Consequently only the small values of |£/^ 4 | are allowed. Thus 
naively one obtains the bound 

4|C/ e4 CVr < 0.25P>^ ugey (Am 2 ) J Dj DHS (Am 2 ), (171) 

which further strengthens the exclusion curve corresponding to the KARMEN bound. A 
detailed and statistically meaningful evaluation of the final combined limit, including also 
the results from other experiments like NOMAD and CHOOZ, leads to the results presented 
in Fig. 40 (Grimus and Schwetz, 2001; Maltoni, Schwetz and Valle, 2001a). The figure 
shows that there is no overlap of the region allowed by the combined bound at 95% CL 
with the region allowed by LSND at 99% CL. For the combined bound at 99% CL there are 
marginal overlaps with the 99% CL LSND allowed region at Am 2 ~ 0.9 and 2 eV 2 , and a 
very marginal overlap region still exists around 6 eV 2 . 

2. (2+2) schemes: active-sterile admixtures 

The main feature of (2+2)-spectra is that either solar or atmospheric oscillations must 
involve the sterile neutrino. Such oscillations are, however, disfavored for both the solar (see 
Sec. IV) and atmospheric (see Sec. V) neutrinos. One expects then that the (2+2) schemes 
are disfavoured. However, as first discussed by Dooling, Giunti, Kang and Kim (2000), 
within (2+2) schemes, oscillations into pure active or pure sterile states are only limiting 
cases of the most general possibility of oscillations into an admixture of active and sterile 
neutrinos. One can wonder then whether some admixture of active-sterile oscillations gives 
an acceptable description of both solar (Giunti, Gonzalez-Garcia and Pena-Garay, 2000) 
and atmospheric (Yasuda, 2000; Fogli, Lisi and Marrone, 2001a, 2001b) data. Combined 
analysis have been performed by Gonzalez-Garcia, Maltoni and Pena-Garay (2001); Maltoni, 
Schwetz and Valle, 2001b; Maltoni, Tortola, Schwetz and Valle, 2002b; 

For the phenomenology of neutrino oscillations, the (2+2)A and (2+2)B schemes are 
equivalent up to the relabeling of the mass eigenstates (or, equivalently, of the mixing 
angles). Thus in what follows we consider the scheme B, where the mass spectrum presents 
the following hierarchy: 

Am 2 G = Am^ < Am 2 tm = Am 43 < Am 2 SND = Am 2 x ~ Am 42 ~ Am 2 31 ~ Am 2 2 . (172) 

Neglecting possible CP phases, and choosing a convention that is convenient for the 
study of solar and atmospheric neutrinos, the matrix (a = e, s, /i, r) can be written as 
follows: 



66 



v u = v 24 v 23 V M V l3 V u V 12 , 



(173) 



where Vij represents a rotation of angle % in the ij plane. Since the parametrizations of 
the leptonic mixing matrix U and of the corresponding lines in V v are the same, and in 
particular involve six mixing angles, we will concentrate below on the allowed values of the 
4x4 neutrino mixing matrix V u . 

This general form can be further simplified by taking into account the negative results 
from the reactor experiments (in particular the Bugey experiment) which in the range of 
Amy relevant to the LSND experiment imply that 

ra 2 +ra 2 = 44+4;$i(r 2 . (174) 

For our purposes, the two angles #13 and #14 can then be safely neglected and the U matrix 
takes the effective form: 



V v 



( Ci 2 Si 2 \ 

— Sl2C23 c 24 C12C23C24 S23C24C34 — S24S34 S23C24S34 + S24C34 

S12S23 -C12S23 C23C34 C23S34 

V S12C23S24 -C12C23S24 -S23S24C34 - C24S34 -S23S24S34 + C24C34 J 



(175) 



The full parameter space relevant to solar and atmospheric neutrino oscillation can be cov- 
ered with three of the angles 6ij in the first quadrant, 9^ G [0, 7r/2] while one (which we 
choose to be 9 U ) is allow to vary in the range — f < #34 < f . 

In this scheme solar neutrino oscillations are generated by the mass-squared difference 
between v 2 and v\ while atmospheric neutrino oscillations are generated by the mass-squared 
difference between u 3 and 1/4. It is clear from Eq. (175) that the survival of solar z/ e 's depends 
mainly on 9 12 while atmospheric z/ e 's are not affected by the four-neutrino oscillations in the 
approximation 6> 13 = 9 U = and neglecting the effect of Am^ in the range of atmospheric 
neutrino energies. The survival probability of atmospheric z/^'s depends mainly on the #34 . 

Thus solar neutrino oscillations occur with a mixing angle 9 12 between the states 

v e -> v a with v a = C23C24 v s + C2 3 ci 4 v a , (176) 
where v a is a linear combination of and v T , 

(«23^ + C 23 S24^t)- (177) 



\/l — C 23 C 24 

We remind the reader that and v T cannot be distinguished in solar neutrino experiments, 
because their matter potential and their interaction in the detectors are equal, due to only 
NC weak interactions. Thus solar neutrino oscillations cannot depend on the mixing angle 
#34 and depend on # 23 and #24 through the combination C23C24. 

Atmospheric neutrino oscillations, i.e. oscillations with the mass difference Am\ A and 
mixing angle #34, occur between the states 

z/3 — > v 1 with z/3 = s 2 3C24^ s + c 2 3^V ~~ s 23S24^r and v 1 = s 24 u s + C 2A V T . (178) 

We learn that the mixing angles #23 and #24 determine two projections. First, the pro- 
jection of the sterile neutrino onto the states in which the solar u e oscillates is given by 
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44 = i - ici 2 - IK2I 2 = IK1I 2 + ra 2 . (179) 

Second, the projection of the z/^ over the solar neutrino oscillating states is given by 

4 = \v^\ 2 + |v; 2 | 2 = 1 - |v^| 2 - 1\/; 4 | 2 (iso) 

One expects s 23 to be small in order to explain the atmospheric neutrino deficit. We will see 
that this is indeed the case. Furthermore, the negative results from the CDHS and CCFR 
searches for //^-disappearance also constrain such a projection to be smaller than 0.2 at 90% 
CL for Am£ SND > 0.4 eV 2 . 

We distinguish the following limiting cases: 

(i) If c 23 = 1 then = V" 2 = 0. The atmospheric = v$ state oscillates into a state 
^7 = C24^ r + S24^s- We will denote this case as restricted (ATM R ). In particular: 

• If c 23 = c 24 = 1, Va\ = Va2 = (Ks = = 0) and we have the limit of pure 
two-generation solar u e — > v s transitions and atmospheric v T transitions. 

• If c 2 4 = then = V" 2 = and V" 3 = V" 4 = 0, corresponding to the limit of pure 
two-generation solar v e — > v T transitions and atmospheric — > v s transitions. 

(ii) If c 2 3 = then = V" 2 = corresponding to the limit of pure two-generation solar 
- * v a with a = fi and there are no atmospheric neutrino oscillations as the projection of 

Vfj, over the relevant states cancels out (V^ 3 = V^ 4 = 0). 

Notice that in the restricted case 623 = 0, there is an additional symmetry in the relevant 
probabilities so that the full parameter space can be spanned by < 34 < |. The reason 
for this is that we now have effectively two-neutrino oscillations for both the solar {y e — > v a ) 
and atmospheric {y^ — > z/ 7 ) cases. 

To summarize, solar neutrino oscillations depend on the new mixing angles only through 
the product c 23 c 2 4 and therefore the analysis of the solar neutrino data in four-neutrino 
mixing schemes is equivalent to the two-neutrino analysis but taking into account that the 
parameter space is now three-dimensional (Am 21 , tan 2 #i 2 , c 2 3 c 24 ). Atmospheric neutrino 
oscillations are affected independently by the angles 823 and # 2 4, and the analysis of the 
atmospheric neutrino data in the four-neutrino mixing schemes is equivalent to the two- 
neutrino analysis, but taking into account that the parameter space is now four- dimensional 
(Am 2 3 , 634, c 23 , c 24 ). Furthermore the allowed ranges of active-sterile oscillations depends on 
the assumed 8 B neutrino flux (Barger, Marfatia, and Whisnant, 2002; Bahcall, Gonzalez- 
Garcia and Pena-Garay, 2002b). Allowing for 8 B neutrino flux larger than the SSM expec- 
tation result into a less stringent limit on the active-sterile admixture. 

As an illustration we show in Fig. 41 the shift in x 2 a s a function of the active-sterile 
admixture \V S \\ 2 + \V S U 2 \ 2 = 

c 23 c 24- From the figure we conclude that the solar neutrino data 
favour pure v e — > u a oscillations but sizable active-sterile admixtures are still allowed. In 
this curve the 8 B flux allowed to take larger values than in the SSM which, as discussed 
above so the active-sterile bound is as model independent as possible. 

Similar analysis can be performed for the atmospheric neutrino data (Gonzalez- Garcia, 
Maltoni and Pena-Garay, 2001) to obtain the allowed regions for the oscillation parameters 
Am 43 and tan 2 # 3 4 from the global analysis for different values of 823 and # 2 4 (or, equiva- 
lently, of the projections |VTi| 2 + \ V" 2 \ 2 and |V^| 2 + |K, 2 | 2 ). The global minimum corresponds 
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to almost pure atmospheric — v T oscillations and the allowed regions become consider- 
ably smaller for increasing values of the mixing angle #23, which determines the size of the 
projection of over the neutrino states oscillating with Amf 3 , and for increasing values of 
the mixing angle #24, which determines the active-sterile admixture in which the almost-u^ 
oscillates. Therefore the atmospheric neutrino data give an upper bound on both mixings 
which further implies a lower bound on the combination C2 3 c 2 4 = |Kil 2 + IKSI 2 - The same 
combination is limited from above by the solar neutrino data. 

Fig. 41 shows the shift in x 2 f° r the analysis of the atmospheric data as a function of 
the active-sterile admixture |Kil 2 + IK2I 2 = c 23 c 24 m the general case (in which the analysis 
is optimized with respect to the parameter sl 3 = \V" X \ 2 + \V^ 2 \ 2 ) as taken from the latest 
article by Maltoni, Tortola, Schwetz, and Valle (2002b). 

In summary, the analysis of the solar data favours the scenario in which the solar oscilla- 
tions in the 1 — 2 plane are u e — u a oscillations, and gives an upper bound on the projection 
of v s on this plane. On the other hand, the analysis of the atmospheric data prefers the 
oscillations of the 3 — 4 states to occur between a close-to-pure and an active (u T ) neutrino 
and gives an upper bound on the projection of the v s over the 3 — 4 states or, equivalently, a 
lower bound on its projection over the 1 — 2 states. From Fig. 41 we see that the exclusion 
curves from the solar and atmospheric analyses are only marginally compatible at 3a CL. 

Thus the conclusion concerning the four-neutrino scenarios is that both 3+1 and 2+2 
schemes seem to be unsatisfactory as explanations of the three experimental pieces of evi- 
dence (see also Maltoni, Tortola, Schwetz and Valle, 2002b). 

VIII. IMPLICATIONS OF THE NEUTRINO MASS SCALE AND FLAVOR STRUCTURE 
A. New Physics 

The simplest and most straightforward lesson of the evidence for neutrino masses is also 
the most striking one: there is NP beyond the SM. This is the first experimental result that 
is inconsistent with the SM. 

Most likely, the NP is related to the existence of singlet fermions at some high energy 
scale that induce, at low energies, the effective terms (13) through the see-saw mechanism. 
The existence of heavy singlet fermions is predicted by many extensions of the SM, most 
noticeably GUTs (beyond SU(5)) and left-right symmetric theories. 

There are of course other possibilities. One could induce neutrino masses without intro- 
ducing any new fermions beyond those of the SM. This requires the existence of a scalar 
A L (l,3) + i, that is a triplet of SU(2) L . The smallness of neutrino masses is then related to 
the smallness of the VEV (A° L ) (required also by the success of the p = 1 relation) and does 
not have a generic natural explanation. 

In left-right symmetric models, however, where the breaking of SU(2)r x U(1)b~l —> 
U(l)y is induced by the VEV of an SU(2) R triplet, A R , there must exist also an SU(2) L 
triplet scalar. Furthermore, the Higgs potential leads to an order of magnitude relation 
between the various VEVs, (A° L ){A R ) ~ v 2 (where v is the electroweak breaking scale), and 
the smallness of (A°) is related to the high scale of SU(2)r breaking. This situation can be 
thought of as a see-saw of VEVs. In this model there are, however, also singlet fermions. 
The light neutrino masses arise from both the see-saw mechanism and the triplet VEV. 
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Neutrino masses could also be of the Dirac type. Here, again, singlet fermions have to be 
introduced, but lepton number conservation needs to be imposed by hand. This possibility is 
disfavored by theorists since it is likely that global symmetries are violated by gravitational 
effects. Furthermore, the lightness of neutrinos (compared to charged fermions) is again 
unexplained. 

Another possibility is that neutrino masses are generated by mixing with singlet fermions 
but the mass scale of these fermions is not high. Here again the lightness of neutrino 
masses remains a puzzle. The best known example of such a scenario is the framework of 
supersymmetry without R-parity. 

Let us emphasize that the see-saw mechanism or, more generally, the extension of the 
SM with non-renormalizable terms, is the simplest explanation of neutrino masses. Models 
where neutrino masses are produced by NP at low energy imply a much more dramatic 
modification of the SM. Furthermore, the existence of see-saw masses is an unavoidable 
prediction of various extensions of the SM. In contrast, many (but not all) of the low energy 
mechanisms are introduced for the specific purpose of generating neutrino masses. 

In this and in the next section, where we discuss the implications of the experimental 
data for theories beyond the SM, we choose to focus on models that explain the atmospheric 
and solar neutrino data through mixing among three active neutrinos. In other words, we 
assume three-neutrino mixing with the oscillation parameters derived in Sec. VII. We do not 
review models that try to incorporate the LSND data by adding light sterile neutrinos and 
we only comment on this possibility in the context of theories with new extra dimensions. 

B. The Scale of New Physics 

Given the relation (14), m v ~ i; 2 /Anp, it is straightforward to use measured neutrino 
masses to estimate the scale of NP that is relevant to their generation. In particular, if there 
is no quasi- degeneracy in the neutrino masses, the heaviest of the active neutrino masses 
can be estimated, 

m h = m 3 ~ V Am atm ~ 0.05 eV. (181) 

(In the case of inverted hierarchy the implied scale is m h = m 2 ~ ^J\Aml tm \ ps 0.05 eV). It 
follows that the scale in the non-renormalizable term (13) is given by 

A N p ~ v 2 /m h « 10 15 GeV. (182) 

We should clarify two points regarding Eq. (182): 

1 . There could be some level of degeneracy between the neutrino masses that are relevant 
to the atmospheric neutrino oscillations. In such a case (see Sec. VI. E), Eq. (181) is modified 
into a lower bound and, consequently, Eq. (182) becomes an upper bound on the scale of 
NP. 

2. It could be that the Zij couplings of Eq. (13) are much smaller than one. In such a 
case, again, Eq. (182) becomes an upper bound on the scale of NP. On the other hand, in 
models of approximate flavor symmetries, there are relations between the structures of the 
charged lepton and neutrino mass matrices that give quite generically Z 33 > rn^jv 2 ~ 10" 4 . 
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We conclude that the likely range of Anp that is implied by the atmospheric neutrino results 
is given by 

10 11 GeV < Anp < 10 15 GeV. (183) 

The estimates (182) and (183) are very exciting. First, the upper bound on the scale of 
NP is well below the Planck scale. This means that there is a new scale in Nature which 
is intermediate between the two known scales, the Planck scale m P \ ~ 10 19 GeV and the 
electroweak breaking scale, v ~ 10 2 GeV. 

It is amusing to note in this regard that the solar neutrino problem does not necessarily 
imply such a new scale. If its solution is related to vacuum oscillations with Am^ ~ 10~ 10 
eV 2 , it can be explained by A N p ~ mp\. However (see Sec. IV.D.4), the favoured explanation 
for the solar neutrino deficit is the LMA solution which again points towards NP scale in 
the range of Eq. (183). 

Second, the scale Anp ~ 10 15 GeV is intriguingly close to the scale of gauge coupling 
unification. We will say more about this fact when we discuss GUTs. 



C. Implications for Flavor Physics 

1. The flavor parameters 

Flavor physics refers to interactions that are not universal in generation space. In the 
Standard Model interaction basis, this is the physics of the Yukawa couplings. In the mass 
basis, this term refers to fermion masses and mixing parameters. There are two main reasons 
for the interest in flavor physics: 

1. The SM has thirteen flavor parameters. These are the nine charged fermion masses 
or, equivalently, Yukawa couplings, Y; L = -^j^'- 

-5 

1 

-4 

i 

" 6 , (184) 

and the four CKM parameters: 

|Vy~0.2, |V/ c6 |~0.04, iKfcl- 0.004, sin<$KM~l- (185) 

One can easily see that the flavor parameters are hierarchical (that is, they have very different 
magnitudes from each other) and all but two (the top- Yukawa and the CP violating phase) 
are small. The unexplained smallness and hierarchy pose the flavor puzzle of the SM. Its 
solution may direct us to physics beyond the SM. 

2. The smallness of flavor changing neutral current (FCNC) processes, such as Amx 
and n — > cy, is explained within the SM by the GIM mechanism. In many extensions of the 
SM there is, generically, no such mechanism. In some cases, experimental bounds on flavor 
changing processes are violated. The solution of such flavor problems leads to refinement of 
the models. In other cases, the model predictions are within present bounds but well above 
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the SM range. Then we can hope to probe this NP through future measurements of flavor 
changing processes. 

Many mechanisms have been proposed in response to either or both of these flavor 
aspects. For example, approximate horizontal symmetries, broken by a small parameter, 
can lead to selection rules that both explain the hierarchy of the Yukawa couplings and 
suppress flavor changing couplings in sectors of new physics. 

In the extension of the SM with three active neutrinos that have Majorana-type masses, 
there are nine new flavor parameters in addition to those of the SM with massless neutrinos. 
These are three neutrino masses, three lepton mixing angles and three phases in the mixing 
matrix. The counting of parameters is explained in Sec. II. D. Of these, four are determined 
from existing measurements of solar and atmospheric neutrino fluxes: two mass-squared 
differences and two mixing angles. This adds significantly to the input data on flavor physics 
and provides an opportunity to test and refine flavor models. 

If neutrino masses arise from effective terms of the form (13), then the overall scale of 
neutrino masses is related to the scale A NP and, in most cases, does not teach us about flavor 
physics. The more significant information for flavor models can then be written in terms 
of three dimensionless parameters whose values can be read from the results of the global 
analysis in Eq. (166); First, the mixing angles that are relevant to atmospheric neutrinos: 

|EV/; 3 |~0.4-0.5; (186) 
Second, the mixing angles that are relevant to solar neutrinos: 

\U el U* 2 \ ~ 0.35 -0.5; (187) 
Third, the ratio between the respective mass-squared differences: 

Am| r lQ- 1 - 10- 3 LMA, 



Am§ 2 | 1 10~ 4 - 10 5 LOW. 



(188) 



In addition, the upper bound on the third mixing angle from the CHOOZ experiments often 
plays a significant role in flavor model building: 

\U e3 \<0.24. (189) 



2. Special features of the neutrino flavor parameters 

There are several features in the numerical values of these parameters that have drawn 
much attention and drove numerous investigations: 

(i) Large mixing and strong hierarchy: The mixing angle that is relevant to the at- 
mospheric neutrino problem, U^, is large, of order one. On the other hand, the ratio of 
mass-squared differences Am^/IAmSy is small. If there is no degeneracy in the neutrino 
sector then the small ratio of mass-squared differences implies a small ratio between the 
masses themselves, 7712/^3 *C 1. (In the case of inverted hierarchy, the implied hierarchy 
is m 3 /m 2 < 1.) It is difficult to explain in a natural way a situation where, in the 2 — 3 
generation sector, there is large mixing but the corresponding masses are hierarchical. Below 
we discuss in detail the difficulties and various possible solutions. 
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(ii) If the LMA solution to the solar neutrino problem holds, then the data can be 
interpreted in a very different way. In this case, the two measured mixing angles are of order 
one. Moreover, Am^/I Am§ 2 | ~ 1CT 2 means that the ratio between the masses themselves 
(which, for fermions, are after all the fundamental parameters) is not very small, m 2 /m 3 ~ 
1CT 1 . Such a value could easily be an accidentally small number, without any parametric 
suppression. If this is the correct way of reading the data, the measured neutrino parameters 
may actually reflect the absence of any hierarchical structure in the neutrino mass matrices 
(Hall, Murayama and Weiner (2000), Haba and Murayama (2001), Berger and Siyeon (2001), 
Hirsch and King (2001)). Obviously, this interpretation is plausible only if the LMA solution 
to the solar neutrino problem holds, and will be excluded if the SMA (small mixing angle) 
or LOW (small mass ratio) solutions hold. Another important test of this idea will be 
provided by a measurement of \U e3 \. If indeed the entries in M v have random values of the 
same order, all three mixing angles are expected to be of order one. If experiments measure 
|£7 e3 | ~ 10 -1 , that is close to the present CHOOZ bound, it can be again argued that its 
smallness is accidental. The stronger the upper bound on this angle will become, the more 
difficult it will be to maintain this view. 

(iii) A special case of large mixing is that of maximal mixing. In a two generation case, 
with a single mixing angle, maximal mixing is defined as sin 2 29 = 1. In the three generation 
case, what we mean by maximal mixing is that the disappearance probability is equivalent 
to that for maximal two neutrino mixing at the relevant mass scale (Barger et al, 1998b). 
Maximal atmospheric neutrino mixing corresponds to 4|£/ m3 | 2 (1 — \U^\ 2 ) = 1, which leads 
to 

\U^\ 2 = 1/2. (190) 

Maximal solar neutrino mixing corresponds to A\U e iU e 2\ 2 = 1, which leads to 

\U el \ 2 = \U e2 \ 2 = l/2, \U e3 \ 2 = 0. (191) 

As seen in Sec VII, present data are consistent with atmospheric neutrino mixing being 
near-maximal and solar neutrino mixing being large but not maximal. The possibility that 
both mixings are near-maximal is, however, not entirely excluded. This scenario, where 
both atmospheric and solar neutrino mixing are near-maximal, means that the structure of 
the leptonic mixing matrix is given to a good approximation by 
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(192) 



The case of equation (192) is commonly called bi-maximal mixing. We would like to make 
the following comments, regarding maximal or bi-maximal mixing: 

1. Theoretically, it is not difficult to construct models that explain near- maximal solar 
neutrino mixing in a natural way. We will encounter some examples below. Experimentally, 
it may be difficult to make a convincing case for near-maximal (rather than just order one) 
solar neutrino mixing (Gonzalez- Garcia, Pena-Garay, Nir and Smirnov, 2001). 

2. It is highly non-trivial to construct models that explain near-maximal atmospheric 
neutrino mixing in a natural way. The reason is that near-maximal mixing is often related 
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to quasi- degeneracy. So when solar neutrino data are also taken into account, approxi- 
mate degeneracy among all three neutrinos is required and models of three quasi-degenerate 
neutrinos are, in general, not easy to construct in a natural way. 

3. The case of bi-maximal mixing is also very challenging for theory. Many of the 
attempts in the literature involve fine-tuning. (Alternatively, the term bi-maximal mixing is 
sometimes used to denote the case of two large, rather than maximal, mixing angles.) 

3. Large mixing and strong hierarchy 

In this subsection, we focus our attention on the large mixing — strong hierarchy problem. 
We explain the challenge to theory from this feature and present various solutions that have 
been proposed to achieve this challenge. 

A large mixing angle by itself should not be a surprise. After all, the naive guess would be 
that all dimensionless couplings [such as the Zij couplings of Eq. (13)] are naturally of order 
one and consequently all mixing angles are of order one. However, the quark mixing angles 
are small, and this situation has led to a prejudice that so would the lepton mixing angles 
be. Second, and more important, flavor models have a built-in mechanism to naturally 
induce smallness and hierarchy in the quark parameters (and perhaps also in charged lepton 
masses). For example, a mechanism that has been intensively studied in the literature is 
that of selection rules due to an approximate symmetry. Within such frameworks, numbers 
of order one are as difficult (or as easy) to account for as small numbers: one can assign 
charges in such a way that small flavor parameters correspond to terms in the Lagrangian 
that carry charge different from zero, while the order one parameters correspond to terms 
that carry no charge. 

The combination of a large mixing angle and strong hierarchy is however a true puzzle. 
To understand the difficulty in this combination, let us assume that both the mixing and 
the hierarchy can be understood from the relevant 2x2 block in the mass matrix for light 
(Majorana) neutrinos. The generic form of this block is 

This matrix would lead to a large mixing angle if 

\B\>\C-A\, (194) 

and to strong mass hierarchy if (for simplicity we assume real entries): 

\AC-B 2 \ <max(A 2 ,5 2 ,C 2 ). (195) 

If we examine the two neutrino sector alone, these conditions mean fine-tuning. Hence the 
challenge to find models where the hierarchy is naturally induced. 

To make the problem even sharper, let us explicitly discuss models with a horizontal 
U{1)h symmetry (Froggatt and Nielsen, 1979). A horizontal symmetry is one where different 
generations may carry different charges, in contrast to the SM gauge group. We assume that 
the symmetry is broken by a small parameter A. To be concrete, we take A = 0.2, close 
to the value of the Cabibbo angle. One can derive selection rules by attributing charge to 



74 



the breaking parameter. Our normalization is such that this charge is H(X) — — 1. While 
coefficients of order one are not determined in this framework, one can derive the parameteric 
suppression of all couplings and consequently have an order of magnitude estimate of the 
physical parameters. Take, for example, the case that the //-charges of the left-handed 
lepton doublets Li are positive. The parameteric suppression of an entry in the neutrino 
mass matrix is determined by the charge that it carries, (M„)y ~ \H(L i )+H(L j )+2H(4>) _ rj^g 
parametric suppression of the mixing angles and mass ratios can then be estimated (see e.g. 
Leurer, Nir and Seiberg, 1993) 

|t^-|~A H ^)- H ^), ^^AW)-^ (i<j). (196) 

If the various generations of left-handed fields (quarks or leptons) carry different if -charges 
then, in general, the (quark or lepton) mixing angles are suppressed. For example, sin6> c ~ A 
is naturally induced if the charges of the first two quark-doublet generations are chosen 
appropriately: H(Qi) — H(Q2) = 1. A mixing angle of order one can be naturally induced 
if the charges of the corresponding lepton-doublet fields are equal to each other. Equation 
(196) shows however that in this class of models, independent of the charge assignments (as 
long as they are all positive), we have (Grossman and Nir, 1995) 

m(vi)/m(uj) ~ \Uij\ 2 . (197) 

Hence, for order one mixing, there is no parametric suppression of the corresponding neutrino 
mass ratio, and no hierarchy induced. 

There is another possibility to induce large lepton mixing which is unique to the case of 
Majorana neutrinos. Here one assigns, for example, opposite charges, H{Lj) = —H(Lj), to 
the relevant lepton-doublet fields. The selection rules lead to the following structure of the 
mass matrix: 

— ( \ cxm ), A,B,C = 0(1). (198) 

(In a supersymmetric framework, the combination of supersymmetry and the horizontal 
symmetry (Leurer et al, 1993, 1994) leads to a situation where either A or C vanish.) 
This mass matrix has a pseudo-Dirac structure and it leads to near-maximal mixing and 
near-degeneracy of masses (see, for example, Nir (2000), Joshipura and Rindani (2000b)): 

sin 2 = 1/2-0 (A 2 ^)l) , ^ = O (A 2 I^)I) . (199) 

Of course, a pseudo-Dirac structure is inconsistent with mass hierarchy. 

Here are a few mechanisms that have been suggested in the literature to induce strong 
hierarchy simultaneously with large mixing angle: 

(i) Accidental hierarchy: the mass matrix (193) has A,B,C = 0(1) and (195) is ac- 
cidentally fulfilled. In the context of Abelian horizontal symmetries, this means that the 
mass ratio is numerically but not parametrically suppressed (Binetruy et al, 1996; Irges et 
al, 1998; Elwood et al, 1998; Vissani, 1998, 2001b; Ellis et al, 1999; Sato and Yanagida, 
2000). We would like to emphasize that this is not an unlikely resolution if the LMA solution 
of the solar neutrino problem is correct. While the oscillation experiments give the ratio 
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of masses-squared, the relevant quantity for theories is the ratio of masses. For the LMA 
solution, m 2 /m 3 ~ 0.1 which is not a particularly small number and could easily arise from 
a combination of order one terms. 

(ii) Several sources for neutrino masses: the leading contribution to the neutrino mass 
matrix, which is responsible to the order one mixing, is rank one. The lighter generation 
masses come from a different, sub-dominant source. This is, for example, the generic sit- 
uation in supersymmetric models without R-parity (see, for example, Banks et ai, 1995; 
Borzumati et al, 1996; Davidson and Losada, 2000). There tree level mixing with neutrali- 
nos can lead to large mixing but gives a mass to only one neutrino generation, while the 
lighter masses arise at the loop level. Another realization of this principle is the scenario of 
a see-saw mechanism with a single right-handed neutrino (Davidson and King, 1998; King, 
1998, 1999, 2000; de Gouvea and Valle, 2001). For another example, within Supersymmetric 
theories, see Borzumati and Nomura, 2001. 

(iii) Large mixing from the charged lepton sector: it is possible that the neutrino mass 
matrix is hierarchical and nearly diagonal in the 2 — 3 sector, and the large mixing is coming 
from the diagonalization of the charged lepton mass matrix. A variety of models have 
been constructed that give this structure of mass matrices, based on discrete horizontal 
symmetries (Grossman et al, 1998; Tanimoto, 1999b), holomorphic zeros (Grossman et ai, 
1998) and U(\)h with two breaking parameters (Nir and Shadmi, 1999). For example, 
Grossman et al. (1998) work in the supersymmetric framework with a horizontal U{1) x 
U(l) symmetry. By an appropriate choice of horizontal charges, they obtain the following 
structure of lepton mass matrices (arbitrary coefficients of order one are omitted in the 
various entries): 

A 8 A 6 A 4 ' 

M y ~ A I , M i± ~ (cj) d ) I A 7 A 5 A 3 I , (200) 




Anp . , ,, , / \ ^7 ^5 ^3 

where the zero entries are a consequence of holomorphy. These matrices lead to 

Ar74/|Am 2 32 |~A 3 , |C/ M3 |~1, (201) 

where the dominant contribution to the mixing angle comes from diagonalization of M(±. 
The solar neutrino parameters of this model correspond to the LMA solution. 

(iv) Large mixing from the see-saw mechanism: it is possible that sterile neutrinos play 
a significant role in the flavor parameters of the light, dominantly active neutrinos. In other 
words, an effective mass matrix of the form (193) with the condition (195) can be a result 
of selection rules applied in the full high energy theory to a mass matrix that includes both 
active and sterile neutrinos and integrating out the heavy, dominantly sterile ones (Barbieri 
et al, 1998; Altarelli and Feruglio, 1998; Eyal and Nir, 1999). It is interesting to note 
that in a large class of such models, the induced hierarchy is too strong for the LMA and 
SMA solutions, unless at least three sterile neutrinos play a role in determining the low 
energy parameters (Nir and Shadmi, 1999). An interesting example for a model of this type 
is presented by Altarelli and Feruglio (1998). They consider a horizontal U(l) symmetry 
broken by two parameters of opposite charges. By an appropriate choice of horizontal 
charges for the lepton fields, they obtain the following neutrino Dirac mass matrix and 
Majorana mass matrix for the sterile neutrinos M N : 
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M D ~ (0 M > A A' 1 , M N ~ A NP 1 A' 2 A' . (202) 
V A A' 1 / V A A ' 1 / 

Integrating out of the three heavy neutrinos, the following Majorana mass matrix for the 
light neutrinos is obtained: 

M v ~^- A 2 A B\, A,B,C = 0(1), (203) 

(204) 

If one were to apply the selection rules directly on M v , the generic structure of Eq. (203) 
would be reproduced, but not the relation between the order one coefficients in Eq. (204). 

(v) A three generation mechanism: approximate L e — — L T . One of the more in- 
teresting frameworks that produces all the observed features of neutrino flavor parameters 
is intrinsically a three generation framework. One applies an approximate L e — — L T 
symmetry to the mass matrices of light, active neutrinos and of charged leptons (Barbieri et 
al., 1998, 1999; Frampton and Glashow, 1999; Joshipura and Rindani, 2000a; Mohapatra et 
al., 2000a; Cheung and Kong, 2000; Shafi and Tavartkiladze; 2000d; Nir and Shadmi,2000; 
Nir, 2000). For the most general case, the symmetry is broken by small parameters, e + and 
e_, of charges +2 and —2, respectively. The lepton mass matrices have the following form: 

/ , v 2 / e_ 1 1 \ / A e A M e_ A T e_ \ 

M v ~ ^p- 1 e+ e+ , M t ~ (fa) A e e+ X, A T , (205) 
ANP V 1 e + e + j VV+ K K J 

where the Aj's allow for generic approximate symmetry that acts on the S , f/(2)-singlet 
charged leptons. The resulting neutrino masses are as follows: 

mi ;2 = m (1 ± C[max(e + , e_)]) , m 3 = mO(e + ). (206) 

Note that the quasi-degenerate pair of mass eigenstates involved in the solar neutrino 
anomaly is heavier than the third, separated state. For the mixing angles, one finds 

\U^\ = 0(1), U e3 = G[max(e + , e_)], \U e2 \ 2 = 1/2 - G[max(e + , e_)]. (207) 

The overall picture is that, somewhat surprisingly, the lepton parameters (186), (187) 
and (188) are not easy to account for with Abelian flavor symmetries. The simplest and 
most predictive models have difficulties in accommodating the large 2 — 3 mixing together 
with the strong 2 — 3 hierarchy. One can find more complicated models that naturally induce 
these parameters, but often at the cost of losing predictive power. In particular, it may be 
the case that, specifically for neutrinos, one cannot ignore the existence of heavy degrees of 
freedom (sterile neutrinos), well beyond the reach of direct experimental production, that 
affect the flavor parameters of the low energy effective theory. If true, this situation would 
mean that measuring the low-energy neutrino parameters cannot by itself make a convincing 
case for the idea of Abelian flavor symmetries. (An exception are models of approximate 
L e - L M - L T symmetry.) 
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Similar difficulties are encountered in the framework of non-Abelian symmetries. Again, 
the simplest models do not work and have to be extended in rather complicated ways (Bar- 
bieri et al, 1999a,b). In some cases, the non-Abelian symmetries can give testable (almost) 
exact relations between masses and mixing angles. For example, the model of Barbieri et al. 
(1999c) predicts sin#i 2 = \Jm e /m^ cos #23 and sin #13 = ^m e /m^ sin 6*23. If it turns out that 
all three light neutrinos are quasi-degenerate, non-Abelian symmetries will become an un- 
avoidable ingredient in the flavor model building, but the task of constructing realistic models 
will be very challenging (see, for example, Carone and Sher, 1998; Fukugita, Tanimoto and 
Yanagida, 1998, 1999; Fritzsch and Xing, 1996,2000; Barbieri et al, 1999d; Wetterich, 1999; 
Tanimoto, 1999a,2000; Tanimoto, Watari and Yanagida, 1999; Wu, 1999a,b,c; Perez, 2000; 
Ma and Rajasekaran, 2001). Radiative corrections are an important issue when examining 
the naturalness of various models that account for quasi-degeneracy among the neutrinos 
(see, for example, Casas et al, 1999, 2000). 

IX. IMPLICATIONS FOR MODELS OF NEW PHYSICS 

In this section we review the implications of neutrino physics for various extensions of the 
SM. We do not attempt to describe all relevant models of NP, but take two representative 
frameworks. We choose to focus on well motivated models that were not constructed for 
the special purpose of explaining the neutrino parameters. Thus, the neutrino parameters 
either test these frameworks or provide further guidance in distinguishing among different 
options for model building within each framework. The two frameworks are the following: 

(i) GUTs: here the overall scale of the neutrino masses and some features of the mixing 
provide interesting tests. 

(ii) Large extra dimensions: in the absence of a high energy scale, the see-saw mechanism 
is not operative in this framework. The lightness of the neutrinos is a challenge. 

A. Grand Unified Theories 

There are two significant facts about the gauge symmetries of the Standard Model and 
the structure of its fermionic representations that motivate the idea of supersymmetric grand 
unification. First, GUTs provide a unification of the SM multiplets (for a recent review, see 
Wilczek, 2001). Second, in the framework of the supersymmetric Standard Model, the three 
gauge couplings unify (Dimopoulos, Raby and Wilczek, 1981; Ibahez and Ross, 1981). The 
unification scale is given by (see, for example, Langacker and Polonsky, 1995) 

A GU ~ 3 x 10 16 GeV. (208) 

Further support for grand unification comes from the flavor sector: the masses of the bottom 
quark and the tau lepton are consistent with equal masses at A GU , as predicted by SU(5) 
and its extensions. 

The evidence for neutrino masses from atmospheric neutrinos, the implied scale of Am^ tm 
and the required mixing of order one should be considered as three further triumphs of the 
grand unification idea. 
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First, as mentioned above, SO(10) theories and their extensions require that there exist 
singlet fermions. Neutrino masses are then a prediction of these theories. 

Second, SO(10) theories naturally give the singlet fermions heavy masses at the SO(10) 
breaking scale and, furthermore, relate the neutrino Dirac mass matrix to the up quark 
mass matrix M u . Specifically, the naive SO(10) relation reads Mjf = M u . It follows that 
the mass of the heaviest among the three light neutrinos can be estimated: 



It only takes that the lightest of the three singlet fermion masses is a factor of fifty below 
the unification scale to induce neutrino masses at the scale that is appropriate for the 
atmospheric neutrino anomaly. In other words, the intriguing proximity of Anp of Eq. (182) 
to Agu of Eq. (208) finds a natural explanation in the framework of SO(10). In this sense 
the see-saw mechanism of Gell-Mann, Ramond and Slansky (1979), and of Yanagida (1979), 
predicted neutrino masses at a scale that is relevant to atmospheric neutrinos. 

The third triumph, \U^\ = 0(1), is more subtle but still quite impressive (Harvey, Reiss 
and Ramond, 1982; Babu and Barr, 1996; Sato and Yanagida, 1998). Unlike the previous 
points that were related to SO(10) GUT, here the consistency relates to SU(5) GUT. SU(5) 
theories relate the charged lepton mass matrix Mu to the down quark mass matrix M^. 
Specifically, the naive SU(5) relation reads M e = Mj. It follows that 



Given the experimental values \V c b\ ~ 0.04 and m s /m b ~ 0.03 we conclude that the naive 
SU(5) relations predict \U^ 3 \ ~ 1. Of course, \U^ 3 \ is also affected by the diagonalization 
of the light neutrino mass matrix, but there is no reason to assume that this contribution 
cancels against the charged lepton sector in such a way that {U^] <C 1. 

Many of the other naive SO(10) relations fail, as do many of the naive SU(5) re- 
lations. Specifically, SO(10) predicts vanishing CKM mixing and mass ratios such as 
m c /m t = m s /m b . SU(5) predicts m s = and = m e . It is possible however that 
these bad predictions are corrected by subleading effects while all the successful predictions 
(particularly m b = m T , m^/m r ~ m s /m b and (t^V^I ~ m s /m b ) are retained since they 
depend on the leading contributions. This is demonstrated in a number of specific GUT 
models (Albright, Babu and Barr, 1998; Albright and Barr, 1998, 1999a, 1999b, 2000a, 
2000b, 2001; Altarelli and Feruglio, 1998, 1999a, 1999b; Altarelli, Feruglio and Masina, 
2000a, 2000b; Berezhiani and Rossi, 1999, 2001; Hagiwara and Okamura, 1999; Babu, Pati 
and Wilczek, 2000; Shafi and Tavartkiladze, 1999, 2000a, 2000b, 2000c; Maekawa, 2001). 

The flavor structure of the first two neutrino generations depends on both the Majorana 
mass matrix for the singlet fermions and the subdominant effects that correct the flavor 
parameters of the first two generation quarks and charged leptons. This part of GUTs 
is therefore much more model dependent. Explicit GUT models were constructed that 
accommodate the various solutions of the solar neutrino problem. We will not describe 
them in any detail here. 

B. Extra Dimensions 

New ideas concerning the possibility of large additional dimensions and the world on a 
brane can lead to sources of neutrino masses that are very different from the ones that we 



m 3 ~ m^/A S o(io) ~ 10 3 eV. 



(209) 



U^Vcb\ ~ m s /m b . 



(210) 
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discussed so far. For example, the small mass could be related to the large volume factor 
of the extra dimensions which suppresses the coupling to bulk fermions (Arkani-Hamed et 
al, 2002; Dienes et al, 1999; Mohapatra et al, 1999; Mohapatra and Perez-Lorenzana, 
2000; Lukas et al, 2000, 2001; Dienes and Sarcevic, 2001; Mohapatra et al, 2000b), to the 
breaking of lepton number on a distant brane (Arkani-Hamed et al, 2002), or to the warp 
factor in the Randall- Sundrum framework (Grossman and Neubert, 2000). In this section, 
we briefly describe these three mechanisms. 

The existence of large extra dimensions does not only provide new ways of generating 
small neutrino masses but can also lead to interesting phenomenological consequences. In 
particular, the phenomenology of matter oscillations in the Sun can be affected (Dvali and 
Smirnov, 1999; Barbieri et al, 2000; Lukas et al, 2000; Ma et al, 2000; Caldwell et al, 2001a, 
2001b). Other phenomenological implications can be used to constrain the parameters of 
the extra dimensions (Faraggi and Pospelov, 1999; Ioannisian and Pilaftsis, 2000; Das and 
Kong, 1999). 

1. Coupling to bulk fermions 

It is possible that we live on a brane with three spatial dimensions that is embedded in a 
spacetime with n additional large spatial dimensions (Arkani-Hamed et al, 1998). This idea 
has the potential of providing an understanding of the hierarchy between the gravitational 
mass scale, M P1 , and the electroweak scale, m Z - The hope is to solve the hierarchy problem 
by avoiding a fundamental high energy (that is Planck or GUT) scale altogether. The 
observed Planck scale, Mp\ = {Gn)~ 1 ^ 2 ~ 10 19 GeV, is related to the fundamental Planck 
scale (most likely the string scale) of the 4 + n dimensional theory, M*, by Gauss law: 

Ml, = M: +2 V n , (211) 

where V n is the volume of the n-dimensional extra space. This picture has dramatic phe- 
nomenological consequences for particle physics and cosmology. Such a situation also poses 
an obvious problem for neutrino physics. If there is no scale of physics as high as (183), 
the see-saw mechanism for suppressing the light neutrino masses cannot be implemented. 
Conversely, if there are singlet neutrinos that are confined to the three dimensional brane 
where the active neutrinos live, one expects that their mass is at the string scale, M*, which 
in these models is much smaller than the four dimensional M P1 , perhaps as small as a few 
TeV, and the resulting light neutrino masses are well above the scales of atmospheric or 
solar neutrinos. 

The first implication within this framework of the evidence for neutrino masses is that 
there better be no singlet fermions confined to the brane. Alternatively, the model must 
include some special ingredients to avoid ordinary see-saw masses. 

On the other hand, it is typical in this framework that there are singlet fermions in 
the bulk. This would be the case, for example, with modulinos, the fermionic partners of 
the moduli fields that are generic to string theories. The crucial point is that the Yukawa 
interaction between bulk singlet fermions and brane active neutrinos is highly suppressed by 
a large volume factor of the n extra dimensions, \/y/Vn- This suppression factor reflects the 
small overlap between the wave functions of the sterile neutrino in the bulk and the active 
one on the brane. By construction, this factor provides a suppression of the neutrino Dirac 
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mass by the ratio v/Mp\. Consequently, large extra dimensions provide a natural source of 
very small Dirac mass terms for the active neutrinos. 

The consequences for neutrino masses and mixing depend on the details of the physics 
of the bulk neutrinos. The possible scenarios are clearly described by Lukas et al. (2000). 
There an explicit derivation of the effective four dimensional action from a five dimensional 
one is given. The final result for the Dirac mass is as follows: 

c = = 



Here Yj is the Yukawa coupling between the lepton doublet Li and a bulk singlet M* is the 
string scale and V n is the volume of the n-dimensional extra space. A usual Dirac mass is of 
order Yv, but in Eq. (212) we see explicitly the (M*/Mp{) suppression factor in the effective 
four-dimensional Yukawa couplings, leading to highly suppressed Dirac mass terms. 

Being a bulk state, the singlet fermion has a whole tower of Kaluza-Klein (KK) associated 
states which are all coupled to the left-handed brane neutrinos. In the simplest scenario and 
for one extra dimension of radius R, the masses of all KK states are determined by the scale 
l/R: 

< = ( 213 ) 

In more general scenarios there can be other bulk mass terms and the masses of the KK 
states receive additional contributions. Then the lightest KK mass can be taken as an 
independent parameter while the mass splitting between the states is still determined by 
the scale l/R. 

Let us denote by M min the lightest mass in the Kaluza-Klein spectrum and by Am the 
mass scale that is relevant to atmospheric or to solar neutrino oscillations. Then one can 
distinguish three cases: 

1. 1/R^$> Am and M min ^> Am. The Kaluza-Klein states play no direct role in these 
oscillations. Their main effect is to give see-saw masses to the active brane neutrinos. 

2. l/R ^> Am and M min < Am. The situation is equivalent to conventional models 
containing a small number of sterile states. In particular, for M m j n = it is a possible 
framework for light Dirac neutrinos. For small but non-vanishing M min it may lead to the 
four-neutrino mixing scenarios. However these scenarios are not particularly favoured. 

3. l/R < Am and M min < Am. A large number of bulk modes can take direct part in 
the oscillation phenomena. This situation modifies in a very interesting way the solution to 
the solar neutrino problem that involves MSW resonance conversion into a sterile neutrino. 
Now the v e can oscillate into a whole set of Kaluza-Klein states. However, as discussed in 
Sec. IV, solar oscillations into sterile states are now strongly disfavoured by the comparison 
of the CC and NC fluxes at SNO. 

In summary, large extra dimensions with sterile bulk fermions provide a natural expla- 
nation for light neutrino masses. In the simplest realization, in which there are no bulk 
mass terms and lepton number conservation is imposed in the bulk-brane couplings, light 
Dirac masses are generated. In the general case, the light neutrinos can be either Dirac or 
Majorana particles, depending on the details of the physics of the bulk neutrinos. 
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2. Lepton number breaking on a distant brane 

In its simplest realization, the mechanism described above provides a natural way of 
generating light Dirac masses for the neutrinos. In the framework of extra dimensions one 
could also generate small Majorana masses by an alternative mechanism in which lepton 
number is broken on a distant brane (Arkani-Hamed et ai, 2002). Imagine that lepton 
number is spontaneously broken at the scale M* on a different brane located a distance r 
from our brane. Further assume that the information about this breaking is communicated 
to our brane by a bulk field with a mass m. For the sake of concreteness we take the case 
that mr ^> 1 and there are two extra dimensions. The resulting Majorana mass for the 
active neutrinos can be estimated as follows: 

m. ~ £e~. (214) 

We learn that the naive see-saw mass, of order v 2 /M*, is suppressed by a small exponential 
factor. The consequences for the neutrino spectrum depend on various model dependent 
features: the number of large extra dimensions, the string scale, the distance between our 
brane and the brane where lepton number is broken, and the mass of the mediating bulk 
field. 

A variant of the above two mechanisms can be implemented if the SM fields are con- 
fined to a thick wall (Arkani-Hamed and Schmaltz, 2000). Dirac masses are suppressed if 
left-handed and right-handed neutrinos are located at different points in the wall and con- 
sequently there is an exponentially small overlap of their wave functions. Majorana masses 
are suppressed if lepton number is spontaneously broken by a VEV that is localized within 
the wall but at some distance from the lepton doublet fields. 



3. The warp factor in the Randall-Sundrum [RS] scenario 

A different set-up for the extra dimensions (leading to a different explanation of the 
v/Mpi hierarchy) was proposed by Randall and Sundrum (1999) [RS]. They considered one 
extra dimension parameterized by a coordinate y = r c 0, where —n < <fi < +n. r c is the 
radius of the compact dimension, and the points (x, <fi) and (x, — 0) are identified. A visible 
brane is located at = n and a hidden one at = 0. This set-up leads to the following 
non-factorizable metric: 

ds 2 = e~ 2kr ^r]^dx^dx u - r 2 d<f) 2 . (215) 

The parameter k is the bulk curvature. All dimensionful parameters in the effective theory 
on the visible brane are suppressed by the warp factor, e = e~ kTc7r . With kr c xs 12, that is 
e ~ 10~ 16 , this mechanism produces physical masses of order v from fundamental masses of 
order M P1 . 

The two mechanisms described above to generate small neutrino masses do not work 
in this scenario because the extra dimensions are small and there is no volume suppression 
factor available. Grossman and Neubert (1999) proposed a different mechanism to generate 
small masses in the RS framework. (For a different mechanism, see Huber and Shafi, 2001.) 
With appropriate choice of orbifold boundary conditions, it is possible to locate the zero 
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mode of a right-handed bulk neutrino on the hidden brane. If the fundamental mass scale 
m of the bulk fermions is larger than half the curvature k of the compact dimension, the 
wave function of the right-handed zero mode on the visible brane is power-suppressed in the 
ratio v/Mpi. Coupling the bulk fermions to the Higgs and lepton doublet fields yields 



Note that generically the relation between the neutrino mass and the weak scale is different 
from the see-saw mechanism (except for the special case m/k = 3/2). 

To summarize, in the presence of large extra dimensions, neutrino masses could be sup- 
pressed by the large volume factor if the left-handed neutrinos couple to the bulk fermions 
(being Dirac for massless bulk fermions), or by the distance to a brane where lepton number 
is broken (which generate Majorana masses). In the RS framework, the suppression can be 
induced by a power of the warp factor. The detailed structure of neutrino mass hierarchy 
and mixing is often related to the parameters that describe the extra dimensions. 

C. Supersymmetry without i?-parity 

Supersymmetry is a symmetry that relates fermions and bosons. It is useful in solv- 
ing the fine-tuning problem of the Standard Model (that is the stability of the small ratio 
between the electroweak breaking scale and the Planck scale against radiative corrections) 
and has many other virtues. Extending the Standard Model by imposing Supersymmetry 
requires that the particle content is also extended. In particular, one must add scalar su- 
perpartners to the three generations of quarks and leptons. (The supersymmetric multiplet 
that contains fermions and bosons with the same gauge quantum number is called a super- 
field.) These modifications make a significant difference as far as the accidental symmetries 
of the Standard Model, Eq. (12), are concerned. One has the option of restoring baryon 
number conservation and total lepton number conservation by imposing a discrete symme- 
try, R p = (_i ) 3 - B + L + 25 ', Models where this symmetry is not imposed, and the accidental 
global symmetries of the Standard Model are all violated (or, at least, either baryon number 
or lepton number is not conserved) are therefore called models of Supersymmetry without 
R-parity. 

Supersymmetry without -R-parity provides a very rich phenomenology. In general, this 
framework has baryon number violation, lepton number violation, and many new sources of 
flavor changing and CP violating interactions. Since total lepton number and lepton flavor 
symmetries are violated, it is not surprising that this framework is relevant to the question 
of neutrino masses and mixing. 

In some sense, the smallness of neutrino masses poses a problem for Supersymmetry 
without R p . If dimensionful R p violating (RPV) couplings were at the electroweak breaking 
scale, and if dimensionless R p violating couplings were of order one, then one neutrino mass 
would be at the electroweak breaking scale, and the other two suppressed by a loop factor 
and bottom and tau Yukawa couplings. Thus the three masses would be well above the 
experimental bounds. This is part of the motivation for considering supersymmetry with 
R p . (Proton decay and neutron-antineutron oscillations provide another source of motivation 
for R p .) In other words, even if R p is not an exact symmetry of Nature, it better be an 




(216) 
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approximate one. Various theories have been considered that lead to suppression of R p 
violating couplings. The question of neutrino masses is usually considered in the framework 
of such special theories. (In particular, one usually assumes that baryon number is a good 
symmetry of the theory. While this assumption has no direct consequences for our discussion 
below, it implies that there is no need to consider constraints from proton decay on the lepton 
number violating couplings.) 

One way to think about lepton number violation in the context of supersymmetry without 
R p is to notice that in such a framework, there is no quantum number to distinguish between 
the down-Higgs and the lepton-doublet superfields. Consequently, there are four flavors of 
the (1, 2)_!/ 2 representation. It is convenient then to denote the four doublets by L a , where 
a = 0, 1, 2, 3. There are bilinear and trilinear couplings that lead to lepton number violation. 
In the superpotential, we have the bilinear /x-terms, and the trilinear A and A' terms: 

Wfi p = H a L a (f) u + Xaf3kL a Lf3£k + \' ajk L a Qjdk. (217) 

Here <p u is the Y = +1/2 Higgs superfield, £ k {k = 1,2,3) are the SU(2) L singlet charged 
lepton superfields, and Q k (d k ) are the SU{2)l doublet (down singlet) quark superfields. 
Among the soft supersymmetry breaking terms in the scalar potential, we have the bilinear 
B terms and m 2 terms and trilinear A terms: 

V^f = B a L a <f) u + m 2 ap L a Ll + A a p k L a Lpl k + A' ajk L a Qjd k + h.c, (218) 

where we use the same notation for superfields in Eq. (217) and for their scalar components 
in Eq. (218). 

There are many ways to define what one means by the down Higgs field, even when 
R p is an approximate symmetry. These ways represent different choices of basis in the 
four-flavor space of L a . Let us denote the Higgs fields by L and the lepton fields by Li 
(i = 1,2,3). Then, for example, one can define L by choosing a basis where /ij = 0, or 
by choosing a basis where (Lj) = (vanishing sneutrino VEVs), or by choosing Li to have 
the three lightest charged mass eigenstates. Physical results are of course independent of 
the choice of basis, but there is much confusion in the literature resulting from the use of 
different conventions. For a clarification of this subject and a very useful basis-independent 
formulation, see Davidson and Losada, 2000. 

The most striking feature of neutrino masses in the context of RPV supersymmetry is 
that there are many different sources for them. Consequently, the hierarchy in neutrino 
masses can be the consequence of the interplay between these different sources. Relations 
between mixing and mass hierarchies that often apply when there is a single source of 
masses are easily modified in this framework. Below is a list of some of the contributions to 
neutrino masses that appear in a generic RPV supersymmetry. The estimates of the sizes of 
the various contributions is given in the convention that (Lj) = and with the simplifying 
assumption that there is a single scale m susy that characterizes all supersymmetry breaking 
and the /z term. 

1. Tree level contribution: 

(M„) y ~ (219) 

''^susy 

2. Sneutrino-neutralino loops: 



84 



3. Slepton-lepton loops: 



(M v )ij ~ -Z-^\ink^jkn — — — — ■ (221) 
o7T 77l SUS y 

4. Squark-quark loops: 

~ (222) 



susy 



There are many additional loop diagrams that involve various two insertions of RPV 
couplings. These could be two bilinear couplings, as in Eq. (220), two trilinear couplings, as 
in Eqs. (221) and (222), or one bilinear and one trilinear coupling. For a complete list, see 
Davidson and Losada, 2000. 

A few comments are in order: 

(i) The relative importance of the various contributions is model dependent. 

(ii) The tree contribution vanishes if the following two conditions are fulfilled (Hall and 
Suzuki, 1984; Banks et al, 1995): 

B a oc m 2 ap np = rh 2 /j, a . (223) 

(iii) Even in case that (223) holds at a high energy scale, misalignment between the B 
and \x terms will be induced by RGE, with the likely consequence that the tree contribution 
is dominant at low energy (Nardi, 1997). 

(iv) The tree contribution gives mass to a single neutrino state. Thus at least two 
neutrino masses are induced by loop corrections. 

(v) If there is a single small parameter that suppresses the leading RPV couplings, that 
is //3///0 ~ B 3 /B ~ \3jk/^ojk ~ ^3jfcMojfc (this is the typical case with an abelian flavor 
symmetry), then the sneutrino related loop of Eq. (220) dominates over the other loop 
contributions (Grossman and Haber, 1999). 

There are many ways in which the lepton parameters deduced from atmospheric and so- 
lar neutrino measurements can be accommodated in the framework of RPV supersymmetry. 
The implications of the atmospheric neutrino results for the tree contributions were discussed 
by Bednyakov, Faessler and Kovalenko, 1998; Chun and Lee, 1999; Diaz et al, 2000; Datta, 
Mukhopadhyaya and Vissani, 2000; Joshipura, Vaidya and Vempati, 2000; Suematsu, 2001. 
Implications, in addition, of solar neutrino results for the A and A' dependent loops were 
investigated by Chun et ai, 1999; Joshipura and Vempati, 1999a, 1999b; Mukhopadhyaya, 
Roy and Vissani, 1998; Choi et ai, 1999; Bhattacharyya, Klapdor-Kleingrothaus and Pas, 
1999; Kong, 1999; Datta, Mukhopadhyaya and Roy, 2000; Takayama and Yamaguchi, 2000; 
Mira et al, 2000. Contributions from misalignment of bilinear terms at tree and loop levels 
was analyzed by Grossman and Haber, 1999, 2001; Hambye, Ma and Sarkar, 2000; Romao 
et al., 2000; Hirsch et ai, 2000. Contributions from trilinear terms only were considered in 
Rakshit, Bhattacharyya and Raychaudhuri, 1999; Drees et ai, 1998. Combinations of bilin- 
ear and trilinear couplings in the loop contributions were pointed out by Choi, Hwang and 
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Chun, 1999; Kaplan and Nelson, 2000. Model independent descriptions of all contributions 
were given by Davidson and Losada, 2000; Abada and Losada, 2000a, 2000b; Choi, Chun 
and Hwang, 2000; Kong, 2000. 

Since there are many independent RPV couplings, it is easy to make an ansatz that 
would be consistent with the data. A more interesting question is whether one can start 
from some theoretical guiding principle, such as an approximate flavor symmetry, and then 
the required mass hierarchy and mixing angles would arise naturally (for an early attempt, 
see, for example, Borzumati et al, 1996). Some of the more interesting attempts in this 
direction are the following: 

1. When jj, and B are misaligned, there are contributions related to bilinear couplings at 
both tree and loop level. The ratio between these contributions is roughly ~ 10 -3 . This 
is a little (but certainly not far) below the hierarchy that is appropriate for (Am^/I Am^l) 1 ^ 2 
(see, for example, Grossman and Haber, 1999). 

2. When \i and B are aligned at a high scale, the ratio between the loop contribution from 
trilinear couplings and the tree contribution from RGE-induced misalignment is generally 
of the appropriate hierarchy (see, for example, Nardi, 1997). 

3. If only trilinear couplings are significant, then the ratio between the leading contri- 
butions is likely to be of order ml /3ml or m s/nT-b, which is just the required hierarchy (see, 
for example, Drees et al, 1998). 

We conclude that future neutrino data will provide further guidance into model building 
in the framework of Supersymmetry without i?-parity. 

X. CONCLUSIONS 

Strong evidence for neutrino masses and mixing has been arising from various neutrino 
oscillation experiments in recent years. First, atmospheric neutrinos show deviation from 
the expected ratio between the v^- and z/ e -fluxes. Furthermore, the flux has strong zenith 
angle dependence. The simplest interpretation of these results is that there are — v T 
oscillations with the following mass and mixing parameters: 

1.4 x 10~ 3 eV 2 < Am 2 tm < 6 x 10~ 3 eV 2 , 0.4 < tan 2 # atm < 3.0. (224) 

The ranges quoted here correspond to the results of the global three-neutrino analysis pre- 
sented in Sec. VII [Eq. (166)]. Second, the total rates of solar neutrino fluxes are smaller 
than the theoretical expectations. Furthermore, the suppression is different in the various 
experiments (which are sensitive to different energy ranges). The recently announced mea- 
surements of NC and CC fluxes at SNO provide a 5.3a signal for neutrino flavor transition 
that is independent of the solar model. The simplest interpretation of these results is that 
there are v e — v a oscillations (where v a is some combination of and u T ) with the following 
set of mass and mixing parameters: 

2.4 x 10" 5 eV 2 < Am 2 & < 2.4 x 10~ 4 eV 2 , 0.27 < tan 2 6 & < 0.79 (LMA), (225) 

where again the ranges quoted correspond to the results of the global analysis given in 
Eq. (166). From the global analysis an upper bound on a third mixing angle arises driven 
mainly by the negative results of the reactor experiments in combination with the deduced 
Am^ tm from the atmospheric neutrino data: 
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sin6> rcact or < 0.24 



(226) 



The smallness of this angle guarantees that the results of the three neutrino analysis com- 
bining the atmospheric and solar neutrino data are close to the two separate two neutrino 
analyses. 

The evidence for neutrino masses implies that the SM cannot be a complete picture of 
Nature. In particular, if the SM is only a low energy effective theory, very light neutrino 
masses are expected. The scale at which the SM picture is not valid anymore is inversely 
proportional to the scale of neutrino masses. Specifically 

m v > V / Am2 tm ~ 0.05 eV =}► A NP < 10 15 GeV. (227) 

We learn that there is a scale of NP well below the Planck scale. 

The scale of NP in Eq. (227) is intriguingly close to the scale of coupling unification. 
Indeed, since GUTs with an SO(10) (or larger) gauge group predict neutrino masses, plau- 
sibly at the 0.1 eV scale, the atmospheric neutrino data can be taken as another piece of 
support to the GUT idea. The large mixing angle in (224) also finds a natural and quite 
generic place in GUTs. 

The measured values of the neutrino flavor parameters are useful in testing various ideas 
to explain the flavor puzzle. Quite a few of the simplest models that explain the smallness 
and hierarchy in the quark sector parameters fail to explain the neutrino parameters. The 
neutrino parameters have some features that are quite unique. In particular, the two mixing 
angles atm and 9 Q are large. As concerns the mass hierarchy, Amg/Am^, the situation is 
still ambiguous. If there is strong hierarchy, it is not easy to accommodate it together with 
large mixing angles. If the hierarchy is mild enough to be consistent with just accidental 
suppression, and if, in addition, the third mixing angle (bounded by reactor experiments) is 
not very small, it could well be that there is no hierarchy in the neutrino parameters at all. 
That would call for flavor frameworks that give a different structure for charged and neutral 
fermion parameters. Other possibilities, such as quasi-degeneracy in masses or bi-maximal 
mixing also call for special structure in the neutrino sector that is very different from the 
quark sector. 

The mass scales involved in Eqs. (224) and (225) have implications to many other frame- 
works of NP. In particular, they can help discriminate between various options in the frame- 
work of models with extra dimensions, models of supersymmetry without R parity, left-right 
symmetric models, etc. 

Another hint for neutrino masses comes from the LSND experiment. Here, how- 
ever, the signal is presently observed by a single experiment, and further experimen- 
tal testing is required. The simplest interpretation of the LSND data is that there are 
v e — > oscillations with Aitil SND = 0(1 eV 2 ) and sin 2 26> L snd = 0(0.003). The fact that 
Abil SND 3> Am 2 tm , Arrig means that the three results cannot be explained with oscillations 
among the three active neutrinos alone. 

If the LSND result is confirmed, then more dramatic modifications of the SM will be 
required. The simplest extension, that is, the addition of a light singlet neutrino, is not 
excluded but it does run into difficulties related to the fact that oscillations into purely 
sterile neutrinos fit neither the atmospheric nor the solar neutrino data. 

The good news is that there has been a lot of progress in neutrino physics in recent years. 
Measurements of both atmospheric and solar neutrino fluxes make the case for neutrino 
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masses and mixing more and more convincing. Many theoretical ideas are being excluded, 
while other have at last experimental guidance in constructing them. The other piece of 
good news is that there is a lot of additional experimental information concerning neutrino 
physics to come in the near future. We are guaranteed to learn much more. 

NOTE ADDED 

The original version of this review was completed in January 2002. The year 2002 has 
been a very active and important year in neutrino physics. In April, the SNO collaboration 
have announced the results of their NC measurement, which established the flavour conver- 
sion of solar neutrinos. The review was updated to include the new effects of the NC SNO 
measurement. In December, the KamLAND collaboration released their first data on the 
measurement of the flux of z/ e 's from distant nuclear reactors (Eguchi et al. 2002). Below 
we briefly present the main implications of this result. 

KamLAND results 

KamLAND have measured the flux of P e 's from distant nuclear reactors. In an exposure 
of 162 ton-yr (145.1 days), the ratio of the number of observed inverse /3-decay events to the 
number of events expected without oscillations is 

^KamLAND = 0.611 ± 0.094 . (228) 

for E Pe > 3.4 MeV. This deficit is inconsistent with the expected rate for massless P e 's at 
the 99.95% confidence level. 

KamLAND have also presented the energy dependence of these events in the form of 
the prompt energy (-E pr0 mpt — E Pe + m p — m n ) spectrum with 17 bins of width 0.425 MeV 
for ^prompt > 0.9 MeV. To eliminate the background from geo-neutrinos, only bins with 
-E'prompt > 2.6 MeV are used. The measured spectrum shows a clear deficit, but there is no 
significant signal of energy-dependence of this effect. 

The KamLAND results can be interpreted in terms of V e oscillations with parameters 
shown in Fig.42a (from Bahcall, Gonzalez-Garcia and Pena-Garay, 2002d). Since KamLAND 
is a disappearance experiment and, furthermore, matter effects are small, the experiment 
is not sensitive to the flavour into which v e oscillates. The allowed region has three local 
minima and it is separated into 'islands.' These islands correspond to oscillations with 
wavelengths that are approximately tuned to the average distance between the reactors and 
the detector, 180 km. The local minimum in the lowest mass island (Am 2 = 1.5 x 10 _5 eV 2 ) 
corresponds to an approximate first-maximum in the oscillation probability (minimum in 
the event rate). The overall best-fit point (Am 2 = 7.1 x 10 _5 eV 2 and tan 2 9 = 0.52,1.92) 
corresponds approximately to the second maximum in the oscillation probability. For the 
same Am 2 , maximal mixing is only slightly less favored, A% 2 = 0.4. Thus with the present 
statistical accuracy, KamLAND cannot discriminate between a large-but-not-maximal and 
maximal mixing. 

Because there is no significant evidence of energy distortion, the allowed regions for the 
KamLAND analysis extend to high Am 2 values for which oscillations would be averaged 
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and the event reduction would be energy independent. These solutions, however, are ruled 
out by the non observation of a deficit at shorter baselines, and in particular by the negative 
results of CHOOZ, which rule out solutions with Am 2 > 1 (0.8) x l(r 3 eV 2 at 3a (99% CL). 

Consequences for solar neutrinos 

The most important aspect of Fig. 42a is the demonstration by the KamLAND results 
that anti-neutrinos oscillate with parameters that are consistent with the LMA solar neutrino 
solution. Under the assumption that CPT is satisfied, the anti-neutrino measurements by 
KamLAND apply directly to the neutrino sector and the two sets of data can be combined 
to obtain the globally allowed oscillation parameters (in fact, the KamLAND results already 
show that CPT is satisfied in the neutrino sector to a good approximation). The results 
of such an analysis are shown in Fig. 42b (from Bahcall, Gonzalez-Garcia and Pena-Garay, 
2002d). 

Comparing Fig 18 with Fig 42b, we see that the impact of the KamLAND results is to 
narrow down the allowed parameter space in a very significant way. In particular, the LMA 
region is the only remaining allowed solution. The LOW solution is excluded at 4.8a and 
vacuum solutions are excluded at 4.9a. The once 'favored' SMA solution is now excluded at 
6.1a. Within the LMA region, the main effect of KamLAND is the reduction of the allowed 
range of Am 2 while the impact on the determination of the mixing angle is marginal. 

Furthermore, various alternative solutions to the solar neutrino problem which predict 
that there should be no deficit in the KamLAND experiment (such as neutrino decay or new 
flavor changing interactions) are now excluded at 3.6a. 

Effect on three-neutrino mixing scenarios 

The v e survival probability at KamLAND in the three-neutrino framework is: 

P(&e v e ) — sin 4 #13 + c °s 4 ^13 

where Lj is the distance of reactor % to KamLAND, and j { is the fraction of the total 
neutrino flux which comes from reactor %. Eq. (229) is obtained with two approximations: 
(a) matter effects are neglected at KamLAND-like baselines, and (b) for Am 2 > 1CT 4 eV 2 , 
the KamLAND energy resolution is not good enough to resolve the corresponding oscillation 
wavelength so the oscillation phase associated to Am^ is averaged to 1/2. 

Eq. (229) shows that a small #i 3 does not affect significantly the shape of the measured 
spectrum. On the other hand, the overall normalization is scaled by cos 4 #13 and this factor 
may introduce a non- negligible effect. Within its present accuracy, however, the KamLAND 
experiment cannot provide any further significant constraint on tan 2 13 and the range in 
Eq. (166) still holds (see also Fogli et al. 2002). 

Thus the main effect of the KamLAND results on the global three neutrino analysis of 
solar, atmospheric and reactor data is to modify the allowed range of Am^ in Eq. (166) to 

5.5 x 10~ 5 < Am^/eV 2 < 1.9 x 10~ 4 . (230) 



sin 2 29 12 sin 2 
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E v [GeV] 



, (229) 
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The allowed range of 6u and the corresponding entries of the mixing matrix U are only 
marginally affected by the inclusion of the KamLAND data. 

Effect on four-neutrino mixing scenarios 

In the 2+2 scenarios discussed in Sec. VII. B, the survival probability for z/ e 's in Kam- 
LAND takes the same form as for two neutrino oscillations and depends only on AmJ 2 and 
#12 (as long as matter effects are neglected). Therefore, KamLAND by itself does not pro- 
vide any restriction on the active-sterile admixtures in the oscillating state. Its impact on 
these scenarios is indirect. By reducing the allowed parameter space for AmJ 2 and 9 12 in 
solar neutrino oscillations, it leaves less room for a sterile component in the solar analysis. 
The effect is, however, not very significant, as can be seen in Fig. 43. Comparing with the 
pre-KamLAND figure 41, we see that compatibility with the atmospheric exclusion curve is 
achieved only at the 3.5a (3<r) level for the solar+KamLAND (solar alone) curve. 
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TABLES 



TABLE I. Characteristic values of L and E for various neutrino sources and experiments. 



Experiment 


L (m) 


E (MeV) 


Am 2 (eV 2 ) 


Solar 


10 10 


1 


10 -io 


Atmospheric 


10 4 - 10 7 


10 2 -10 5 


ict 1 - io~ 4 


Reactor 


10 2 - 10 3 


1 


10-2 _ 10-3 


Accelerator 


10 2 


10 3 -10 4 


>0.1 


LBL Accelerator 


10 5 - 10 6 


10 4 


10~ 2 - 10" 3 



TABLE II. Effective potentials for v e (upper sign) and v e (lower sign) in various media. 



medium 


V c 


v N 


e + and e~ 


±V2G F (N e - N s ) 


T^f(A^ e -iVe)(l-4sin 2 9 W ) 


p and p 





±^f(Ap- Ap)(l-4sin 2 M 


n and n 







Neutral (N e = N p ) 


±V2G F N e 





TABLE III. SSM predictions: solar neutrino fluxes and neutrino capture rates for the different 
experiments, with la uncertainties. 



Source 


Flux 


CI 


Ga SK SNO(CC) 


SNO(NC) 

- 2 S - X ) 




(10 10 cm- 2 s- 1 ) 


(SNU) 


(SNU) (10 6 cm-V 1 ) (10 6 cm" 



pp 

pep 
hep 
7 Be 
8 B 

13 N 
15q 
17p 

Total 
Measured 



5.95 ( i.ool8;81 



1-0018:18 



1.40 x 10" 2 (l.OOlg 
9.3 x 10" 7 
4.77 x lO" 1 
5.05 x 10- 4 
5.48 x 10~ 2 
4.80 x 10~ 2 
5.63 x 10- 4 



0.015 
015 



oo +0 - 2 ° N 

uu -0.16 , 

oo +0 - 21 
oo +0 - 25 

uu -0.19 , 

+0.25^ 
-0.25, 



00 



Me 



■d 



0.0 
0.22 
0.04 
1.15 
6.76 
0.09 
0.33 
0.0 



7.6 



+1.3 
-1.1 



69.7 
2.8 
0.1 
34.2 
12.2 
3.4 
5.5 
0.1 

128 



0.0 
0.0 
0.0093 
0.0 
5.05 
0.0 
0.0 
0.0 



0.0 
0.0 
0.0093 
0.0 
5.05 
0.0 
0.0 
0.0 



-9 



5.05: 



-1.01 
-0.81 

2.35 ±0.083 



5.05: 



-1.01 



0.81 

1.76 ±0.11 5.09 ±0.64 



SSM 



2.56 ±0.226 70.8 ±4.3 
0.337 ± 0.065 0.55 ± 0.048 0.465 ± 0.094 0.348 ± 0.073 1.01 ± 0.23 
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TABLE IV. Best-fit global oscillation parameters with all solar neutrino data. This table cor- 
responds to the global solution illustrated in Fig. 18. Sterile solutions do not appear in Fig. 18 
because their Xmm * s l ar g e - 





Solution 


Am 2 [eV 2 ] 


tan 2 6 


/b 


X 2 min (46 d.o.f) 


LMA 


5.0 x 10~ 5 


4.2 x 10" 1 


1.07 


45.5 


LOW 


7.9 x 10~ 8 


6.1 x 10" 1 


0.91 


54.3 


VAC 


4.6 x 10~ 10 


1.8 x 10° 


0.77 


52.0 


SMA 


5.0 x 10~ 6 


1.5 x 10~ 3 


0.89 


62.7 


Just So 2 


5.8 x 10~ 12 


1.0 x 10° 


0.46 


86.3 


Sterile VAC 


4.6 x 10~ 10 


2.3 x 10° 


0.81 


81.6 


Sterile SMA 


3.7 x 10~ 6 


4.7 x 10" 4 


0.55 


89.3 



TABLE V. 90% CL limit on the neutrino oscillation probabilities from the negative searches at 
short baseline experiments. 



Experiment 


Beam 


Channel 


Limit (90%) 


Am 2 lin (eV 2 ) 


Ref. 


CDHSW 


CERN 


V(l -»■ V(i 


P w > 0.95 


0.25 


Dydak et al, 1984 


E776 


BNL 


V(i -»■ 


P efl < 1.5 x 10- 3 


0.075 


Borodovsky et al, 1992 


E734 


BNL 


V(i -»■ 


P e/i < 1.6 x lO- 3 


0.4 


Ahrens et al, 1987 


KARMEN2 


Rutherford 


Vfj, -»■ v e 


P eti < 6.5 x 10~ 4 


0.05 


Wolf et. al, 2001 


E531 


FNAL 


-> V T 


P^ T < 2.5 x lO" 3 


0.9 


Ushida et al, 1986 


CCFR 


FNAL 


-> V T 


P^r < 4 X lO- 3 


1.6 


McFarland et al, 1995 






u e -»■ U T 


P er < 0.10 


20 


Naples et al, 1999 






V(i -»■ v e 


P^e < 9. X lO" 4 


1.6 


Romosan et al, 1997 


Chorus 


CERN 


-> zv r 


P Mr < 3.4 x 10~ 4 


0.6 


Eskut et al, 2001 






U e -> I/ T 


P eT < 2.6 X 10~ 2 


7.5 


Eskut et al, 2001 


Nomad 


CERN 




P^ T < 1.7 X 10~ 4 


0.7 


Astier et al, 2001 






zv e -> i/ r 


P er < 7.5 X lO" 3 


5.9 


Astier et al, 2001 






^ -> 


P^e < 6 X 10" 4 


0.4 


Valuev et al, 2001 



TABLE VI. 90% CL lower bound on P ee and Am 2 sensitivity from searches at reactor experiments. 



Experiment 


L 


Limit (90%) 


A< in (eV 2 ) 


Ref. 


Bugey 


15,40 m 


> 0.91 


lO" 2 


Achkar et al, 1995 


Krasnoyarsk 


57, 230 m 


> 0.93 


7 x lO" 3 


Vidyakin et al, 1994 


Gosgen 


38, 48, 65 m 


> 0.9 


0.02 


Zacek, C, 1986 


CHOOZ 


1 km 


> 0.95 


7 x 10~ 4 


Apollonio et al, 1999 


KamLAND 


150-210 km 


> 0.85 


5 x 10~ 5 


Piepke et al, 2001 
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FIG. 1. The characteristic form of an excluded region from a negative search with fixed L/E 
and of an allowed region from a positive search with varying L/E in the two-neutrino oscillation 
parameter plane. 




FIG. 2. Effective masses acquired in the medium by a system of two massive neutrinos as a function 
of the potential A [see Eq. (60)]. 
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FIG. 3. The mixing angle in matter for a system of two massive neutrinos as a function of the 
potential A for two different values of the mixing angle in vacuum [see Eq. (61)]. 




10 4 10 3 10 2 10 1 1 10 

tcm 2 -tf 

FIG. 4. Isocontours of the survival probablity P ee in the Sun. Also shown is the limit of applica- 
bility of the adiabatic approximation Q = 1 (dashed line). 
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R/ Rsun 

FIG. 5. The solar matter density profile for the BPOO model (full line), and the exponential 
approximation (dashed line). 




4E/Am 2 (eV~ 1 ) 

FIG. 6. The survival probablity for a v e state produced in the center of the sun as a function of 
E/ Am 2 for various values of the mixing angle. 
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FIG. 7. The pp chain in the Sun. 
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FIG. 8. The CNO cycle in the Sun. 
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Neutrino Energy (MeV) 

Neutrino fluxes from the pp chain reactions as a function of the neutrino energy. 




FIG. 10. Production point distribution for the pp chain neutrinos as a function of the distance 
from the solar center. 
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FIG. 11. The electron recoil energy spectrum measured in SK normalized to the SSM prediction, 
and the expectations for the best fit points for the LMA, SMA, LOW and VAC solutions in Fig. 16. 
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FIG. 12. The zenith angle dependence of the solar neutrino flux (statistical error only). The width 
of the night-time bins was chosen to separate solar neutrinos that pass through the Earth's dense 
core (cos# z > 0.84) from those that pass through the mantle (0 < cos^ < 0.84). The horizontal 
line shows the average flux. 
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FIG. 13. Seasonal variation of the solar neutrino flux (statistical errors only). The curve shows 
the expected seasonal variation of the flux induced by the eccentricity of the Earth's orbit. 




4E/Am 2 CeV -1 ) 

FIG. 14. P ee as a function of AE/Am 2 . In panel (a) matter effects are ignored [see Eq. (106)]. 
Panels (b) (large mixing angle) and (c) (small mixing angle) take into account the MSW effect [see 
Eq. (108)]. 
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FIG. 15. Isocontours of the day-night asymmetry at SK. 
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FIG. 16. Allowed oscillation parameters (at 90, 95, 99 and 99.7% CL) from the analysis of the 
total event rates of the Chlorine, Gallium, SK and SNO CC experiments. The best fit point is 
marked with a star. 
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FIG. 17. Excluded oscillation parameters at 95, 99 (shadowed region) and 99.7% CL from the 
analysis of the day-night spectrum data. 
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FIG. 18. Allowed oscillation parameters (at 90, 95, 99 and 99.7% CL) from the global analysis of 
the solar neutrino data. The best fit point (LMA active) is marked with a star. 
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FIG. 19. Event rates as a function of neutrino energy for fully contained events, stopping muons, 
and through-going muons at SuperKamiokande. 
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FIG. 20. The double ratio of to v e events, data divided by theoretical predictions, for various 
underground atmospheric neutrino detectors. 
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FIG. 21. Zenith angle distribution of SuperKamiokande 1289 days data samples. Dots, solid line 
and dashed line correspond to data, MC with no oscillation and MC with best fit oscillation 
parameters, respectively. 
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FIG. 22. The main components of the primary cosmic ray spectrum. 
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FIG. 23. Atmospheric neutrino fluxes as a function of the neutrino energy. 
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FIG. 24. Zenith angular dependence of the atmospheric neutrino fluxes. 
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FIG. 25. The various contributions to the neutrino-nucleon cross section. 
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FIG. 26. The opening angle between the incoming neutrino and the produced charged lepton in 
QE interactions as a function of the charged lepton energy. 
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FIG. 27. as a function of the zenith angle for maximal mixing of with v T (upper panel), v e 
(middle panel) and v s (lower panel). For Am 2 = 5 x 10 -3 eV 2 the curves correspond to E u of 20, 
40, 60 and 80 GeV. The dashed curves are calculated with the approximation of constant average 
densities in the mantle and in the core of the Earth. 
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FIG. 28. The status of the Vn — > v e oscillation solution to the atmospheric neutrino anomaly. 



118 




10 
10 



-2 



-3 



10 



-4 



10 



-2 



10 



-3 



10 



-4 




i i i i i i i i i i i i i i 



J l_l I I I L 



J I I I I I I l_l I L_ 




I I I I I I I I I I I I I I I I I 




J I I l_l I I I I I I L_L 



0,2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 

sin 2 -# 



CO 



c 



c 

o 
o 



en 
c 

o 

CD 



FIG. 29. Allowed regions (at 90, 95 and 99 % CL) from the analysis of various partial data samples 
of atmospheric neutrinos for the oscillation channels Vu — > v T and Vn — * v s . The best fit points are 
marked with a star (see text for details). 
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FIG. 30. Allowed regions (at 90, 95 and 99 % CL) from the analysis of the full data sample of 
atmospheric neutrinos for the oscillation channels — > v T and za, — * v s . The best fit points are 
marked with a star (see text for details). Also shown are the expected sensitivities from LBL 
experiments. 
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FIG. 31. Zenith angular dependence for through-going muon at SuperKamiokande compared with 
the predictions in the case of no-oscillation and for the best fit points of — » v T (Am 2 = 2.6 x 10~ 3 
eV 2 , sin 2 29 = 0.97) and -> v s (Am 2 = 3 x 10~ 3 eV 2 , sin 2 9 = 0.61) oscillations. 
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FIG. 32. Excluded regions at 90% in the v„ — > v T and v e — > v T channels from SBL experiments. 
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FIG. 33. Allowed regions (at 90 and 99 % CL) for v e — > i/„ oscillations from the LSND experiment 
compared with the exclusion regions (at 90% CL) from KARMEN2 and other experiments. The 
90 % CL expected sensitivity curve for MinimBoNE is also shown. 
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FIG. 34. Excluded regions at 90% for v e oscillations from reactors experiments and the expected 
sensitivity from the KamLAND experiment. 
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FIG. 35. Summary of the present pieces of evidence for neutrino masses and mixing as obtained 
from 2-v oscillation analyses. The allowed regions correspond to 90% and 99% CL for 2 d.o.f. We 
also show relevant bounds from laboratory experiments. 
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FIG. 36. Allowed regions (at 90, 95, and 99% CL) in the (Am^, tan 2 9 U ) plane from the global 
analysis of the solar neutrino data in the framework of three-neutrino oscillations for various values 
of sin 2 #13. The global best fit point is marked by the star. 
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FIG. 37. Dependence of A% 2 on sin 2 #13 in the analysis of the atmospheric, solar and CHOOZ 
neutrino data. The dotted horizontal line corresponds to the 3<r limit for a single parameter. 
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FIG. 38. Allowed regions (at 90, 95, 99, and 99.7% CL) in the (Am| 2 , tan 2 6> 23 ) plane from the 
global analysis of the atmospheric neutrino data in the framework of three-neutrino oscillations for 
various values of tan 2 #13. The global best fit point is marked by the star. 
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FIG. 39. The six types of 4-neutrino mass spectra. The different distances between the masses on 
the vertical axes represent the different scales of mass-squared differences required to explain solar, 
atmospheric and LSND data with neutrino oscillations. 




FIG. 40. Upper bounds (at 95 and 99% CL) on A^. e = 4\U ei U ^ in the context of (3+l)-schemes. 
The shaded regions are the regions allowed by LSND at 90 and 99% CL. 
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FIG. 42. (a) Allowed regions (90%, 95%, 99%, and 99.7% CL) for anti-neutrino oscillations as 
implied by the KamLAND result. Results from CHOOZ experiment are also included and rule out 
solutions with Am 2 > 1(0.8) x 10~ 3 eV 2 at 3<r (99% CL) . (b) Allowed oscillation parameters from 
the combined analysis of solar and KamLAND data. 
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